| Re: storing in array [message #60233 is a reply to message #60232] |
Mon, 12 May 2008 03:03   |
Spon
Messages: 178
Registered: September 2007
|
Senior Member |
|
|
On May 11, 10:26 pm, kishore1...@gmail.com wrote:
> On May 11, 5:11 am, Spon <christoph.b...@gmail.com> wrote:
>
>
>
>> On May 11, 12:48 am, kishore1...@gmail.com 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
>
> Hi Chris,
> Thanks for your quick reply.
> I am using this routine in my program "vfm_feature_flags,36282"
> This one is changing every time. Just I want to store o/p information
> into a string.
> Your program is simple, but I could not understand clearly. For
> example I want to store "Feature Subtype : cirrus (transparent) " and
> this one is changing every time. Some times cirrus and some times
> convective ..
>
> Could you give me little clearly then I can follow.
>
> Thanking you,
>
> Kishore
------
First, define your output structure to include every possible field
and subfield,
For example, after this:
feature_type = 0
feature_type_qa = 0
ice_water_phase = 0
ice_water_phase_qa = 0
feature_subtype = 0
; etc.
add this:
Output = { feature_type:"", $
feature_type_qa:"", $
ice_water_phase:"", $
ice_water_phase_qa:"", $
feature_subtype:"", $
; etc. , $
}
------
Then, fill your output structure with strings wherever you've used
PRINT.
For example, change this:
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"
To this:
case feature_type of
0 : Output.FeatureType = 'invalid (bad or missing data)'
1 : Output.FeatureType = 'clear air'
2 : begin
Output.FeatureType = 'cloud'
case feature_subtype of
0 : Output.FeatureSubType = 'low overcast, transparent'
1 : Output.FeatureSubType = 'low overcast, opaque'
...
HELP, Output, /STRUCTURE
RETURN, Output
END
This will give you an 'anonymous structure' (look it up in the
documentaton), rather than a string array, as your output. To me it's
a much more managable piece of information, though it takes slightly
longer to get your head around how to juggle with it...
Plus, if you want you can probably concatenate all you structures into
a structure array and you'll have one single output file for all your
string data.
Regards,
Chris
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
