On May 11, 12:48 am, kishore1...@[EMAIL PROTECTED]
wrote:
> Hi,
>
> I am new in IDL language. For reading the CALIPSO satellite data sets
> in that what type of cloud type for that they provided read routine.
> If I use this routine it is printing on IDL log window. How to store
> this information into an string array.
> This is the routine
> vfm_feature_flags,36282
> The output is like this:
> Bit set: 2
> Bit set: 4
> Bit set: 5
> Bit set: 6
> Bit set: 8
> Bit set: 9
> Bit set: 11
> Bit set: 12
> Bit set: 16
> Feature Type : cloud
> Feature Subtype : cirrus (transparent)
> Feature Type QA : high
> Ice/Water Phase : ice
> Ice/Water Phase QA: high
> Cloud/Aerosol/PSC Type QA : not confident
> Horizontal averaging required for detection: 20 km
>
> The above information i want to store in string array.
>
> ;************************routine********************
> pro vfm_feature_flags,val
> ; this routine demonstrates how to read and extract values from a
> feature
> ; classification flag 16-bit integer value in CALIPSO Level 2
> Vertical
> ; Feature Mask files
> ;
> ; INPUT:
> ; val - the feature classification flag value to be decoded
> ;
> ; OUTPUT:
> ; all information is printed into the IDL log window
>
> print, val
>
> feature_type = 0
> feature_type_qa = 0
> ice_water_phase = 0
> ice_water_phase_qa = 0
> feature_subtype = 0
> cloud_aerosol_psc_type_qa = 0
> horizontal_averaging = 0
>
> for i=0,15 do begin
> if ((val and 2L^i) NE 0) then begin
> print,'Bit set: ',i+1
> case i+1 of
> 1 : feature_type = feature_type + 1
> 2 : feature_type = feature_type + 2
> 3 : feature_type = feature_type + 4
> 4 : feature_type_qa = feature_type_qa + 1
> 5 : feature_type_qa = feature_type_qa + 2
> 6 : ice_water_phase = ice_water_phase + 1
> 7 : ice_water_phase = ice_water_phase + 2
> 8 : ice_water_phase_qa = ice_water_phase_qa + 1
> 9 : ice_water_phase_qa = ice_water_phase_qa + 2
> 10 : feature_subtype = feature_subtype + 1
> 11 : feature_subtype = feature_subtype + 2
> 12 : feature_subtype = feature_subtype + 4
> 13 : cloud_aerosol_psc_type_qa = cloud_aerosol_psc_type_qa + 1
> 14 : horizontal_averaging = horizontal_averaging + 1
> 15 : horizontal_averaging = horizontal_averaging + 2
> 16: horizontal_averaging = horizontal_averaging + 4
> else:
> endcase
> endif
> endfor
>
> case feature_type of
> 0 : print,"Feature Type : invalid (bad or missing data)"
> 1 : print,"Feature Type : clear air"
> 2 : begin
> print,"Feature Type : cloud"
> case feature_subtype of
> 0 : print, "Feature Subtype : low overcast, transparent"
> 1 : print, "Feature Subtype : low overcast, opaque"
> 2 : print, "Feature Subtype : transition stratocumulus"
> 3 : print, "Feature Subtype : low, broken cumulus"
> 4 : print, "Feature Subtype : altocumulus (transparent)"
> 5 : print, "Feature Subtype : altostratus (opaque)"
> 6 : print, "Feature Subtype : cirrus (transparent)"
> 7 : print, "Feature Subtype : deep convective (opaque)"
> else : print,"*** error getting Feature Subtype"
> endcase
> end
> 3 : begin
> print,"Feature Type : aerosol"
> case feature_subtype of
> 0 : print, "Feature Subtype : not determined"
> 1 : print, "Feature Subtype : clean marine"
> 2 : print, "Feature Subtype : dust"
> 3 : print, "Feature Subtype : polluted continental"
> 4 : print, "Feature Subtype : clean continental"
> 5 : print, "Feature Subtype : polluted dust"
> 6 : print, "Feature Subtype : smoke"
> 7 : print, "Feature Subtype : other"
> else : print,"*** error getting Feature Subtype"
> endcase
> end
> 4 : begin
> print,"Feature Type : stratospheric feature--PSC or
> stratospheric aerosol"
> case feature_subtype of
> 0 : print, "Feature Subtype : not determined"
> 1 : print, "Feature Subtype : non-depolarizing PSC"
> 2 : print, "Feature Subtype : depolarizing PSC"
> 3 : print, "Feature Subtype : non-depolarizing aerosol"
> 4 : print, "Feature Subtype : depolarizing aerosol"
> 5 : print, "Feature Subtype : spare"
> 6 : print, "Feature Subtype : spare"
> 7 : print, "Feature Subtype : other"
> else : print,"*** error getting Feature Subtype"
> endcase
> end
> 5 : print,"Feature Type : surface"
> 6 : print,"Feature Type : subsurface"
> 7 : print,"Feature Type : no signal (totally attenuated)"
> else : print,"*** error getting Feature Type"
> endcase
>
> case feature_type_qa of
> 0 : print,"Feature Type QA : none"
> 1 : print,"Feature Type QA : low"
> 2 : print,"Feature Type QA : medium"
> 3 : print,"Feature Type QA : high"
> else : print,"*** error getting Feature Type QA"
> endcase
>
> case ice_water_phase of
> 0 : print,"Ice/Water Phase : unknown/not determined"
> 1 : print,"Ice/Water Phase : ice"
> 2 : print,"Ice/Water Phase : water"
> 3 : print,"Ice/Water Phase : mixed phase"
> else : print,"*** error getting Ice/Water Phase"
> endcase
>
> case ice_water_phase_qa of
> 0 : print,"Ice/Water Phase QA: none"
> 1 : print,"Ice/Water Phase QA: low"
> 2 : print,"Ice/Water Phase QA: medium"
> 3 : print,"Ice/Water Phase QA: high"
> else : print,"*** error getting Ice/Water Phase QA"
> endcase
>
> if (cloud_aerosol_psc_type_qa eq 0) then begin
> print,"Cloud/Aerosol/PSC Type QA : not confident"
> endif else begin
> print,"Cloud/Aerosol/PSC Type QA : confident"
> endelse
>
> case horizontal_averaging of
> 0 : print,"Horizontal averaging required for detection: not
> applicable"
> 1 : print,"Horizontal averaging required for detection: 1/3 km"
> 2 : print,"Horizontal averaging required for detection: 1 km"
> 3 : print,"Horizontal averaging required for detection: 5 km"
> 4 : print,"Horizontal averaging required for detection: 20 km"
> 5 : print,"Horizontal averaging required for detection: 80 km"
> else : print,"*** error getting Horizontal averaging"
> endcase
> end
Rather than a string array, I would use an anonymous structure for
this sort of information:
e.g.
Data = {Type:"", Subtype:"", QA:"", Phase:"", PhaseQA:"", TypeQA:""} ;
etc
This creates a set of empty string fields.
then you can read in your strings like this:
data.type = 'cloud'
and retrieve your information like this:
help, data, /struct
IDL> ** Structure <ff8a18>, 6 tags, length=72, data length=72, refs=1:
TYPE STRING 'cloud'
SUBTYPE STRING ''
QA STRING ''
PHASE STRING ''
PHASEQA STRING ''
TYPEQA STRING ''
Reading in your strings could probably be done into a string array
either if you really want to:
Result = STRARR[n]
FOR i = 0, n-1 DO BEGIN
Result[i] = 'Cloud'
ENDFOR ; etc.
Also, you can probably use format codes to convert your bytes/integers
to set bits, something along the lines of:
CloudBits = STRING(CloudByte, FORMAT='(B0)')
NCloudBits = STRLEN(CloudBits)
Which I can imagine would make determining your string contents a bit
easier. I can't remember exactly, and I'm stuck with demo-mode only
today(!), so you'll have to play around with it yourself, unless
someone more knowledgeable jumps in with a fuller explanation.
Good luck!
Chris
|
