Læser du BSQ-filformat i IDL?

Bruger du IDL med ENVI eller af sig selv?

Hvis du bruger IDL med ENVI (langt den nemmeste hvis du har ENVI), så fungerer standard ENVI-fillæsningskommandoer. F.eks.

ENVI_SELECT, fid=fid, dims=dims, pos=pos ; Shows a dialog box to allow you to select a file image = ENVI_GET_DATA(fid=fid, dims=dims, pos=pos) ; Gets the data into an array called image 

Hvis du ikke har ENVI, er det lidt sværere. Du kan læse BSQ-data ved hjælp af standard IDL-kommandoer til læsning af filer (se http://www.msi.umn.edu/software/idl/tutorial/idl-rwd.html#Reading%20Binary%20Data , men du skal kender dimensionerne på billedet (prøver, linjer og bånd) for at være i stand til at opdele det fornuftigt. Disse oplysninger gemmes ofte i en .hdr-fil, som ENVI læser automatisk, når du bruger metoden ovenfor.

IDL læser båndsekvens (BSQ) binær filformatstruktur som [Kolonne, række, bånd]

 BSQ, BIL and BIP are multiband encodings where: a 8-bit Grey scale image would encode 1 band of data a 24-bit RGB additive color image would be 3 bands of data a 32-bit CMYK subtractive color image would be 4 bands of data Discrete multispectral images, or contiguous hyperspectral images would have increasing numbers of bands. 

Så som bemærket af @robintw skal du kende dimensionerne og kilden / sensoren på dit billede. Bortset fra BSQ-kodning skal du vide fire ting for at læse dataene ind i en IDL interaktiv matrix. Du har brug for kolonnetælling, rækkeoptælling, antal bånd og datatypen. Naturligvis for ethvert geospatialt arbejde uden for IDL skal du kende den nominelle pixelstørrelse og koordinater på mindst 3 kontrolpunkter til billedregistrering Og yderligere detaljer til s opulate en ESRI Grid World-fil.

For afkodning af dine BSQ-data er IDL-binære datatypevalg

 0 Undefined 1 Byte 2 Integer 3 Longword Integer 4 Floating point 5 Double-precision floating 6 Complex 7 String 8 Structure 9 Double-precision complex 10 Pointer 11 Object reference 12 Unsigned Integer 13 Unsigned Longword Integer 14 64-Bit Integer 15 Unsigned 64-Bit Integer 

Så med minimumsinfoen ved hånden: For at “åbne” dataene opretter du en ny IDL-procedure

 make array of (columns,rows,bands) assign to working_array open the bsq_image for reading nested loops to read-in binary_values for each band for each row for each column read the binary_value of type from the bsq_image write the data to the working_array display working_array ;;ordered as needed for false color aesthetics 

Kommentarer

• @Yudine, skal der være en vedhæftet fil eller nogle detaljer?
• @ Stuart, dette er filoplysningerne ENVI-beskrivelse = {File Resize Result, x resize factor: 1.000000, y resize factor: 1.000000. [Thu Jan 13 18:15:18 2011]} prøver = 338 linjer = 308 bånd = 125 header offset = 0 filtype = ENVI Standard datatype = 2 interleave = bsq sensortype = Ukendt byte rækkefølge = 0 x start = 459 y start = 1735 Antal bånd = 125
• @Yudine, ser ud til at du har alle de ting, du har brug for for at læse billedet i et IDL-array. 338 prøver – > kolonner, 308 linjer – > rækker, 125 bånd. Datatype 2 heltal (sandsynligvis 0-255). Og byteordren er 0, dvs. Little Endian, så direkte læses på en Windows-pc (bliver nødt til at transformere (Swap_If_Big_Endian) hvis en 1. Så lav en ny Pro i IDLDE GUI, og følg pseudokoden. En læsning af Michael Todd ' s ovennævnte note med linket til en klasseøvelse på UIdaho skal hjælpe dig med syntaks for koden. Derefter vil dit eneste spørgsmål være at beslutte hvilke bånd i arrayet der skal vises, 125 til kombinere.

Jeg har ikke erfaring med enten bsq eller idl, men da gdal understøtter bsq , kan du prøve at bruge gdal_translate til at konvertere din fil til en geotiff. IDL ser ud til at understøtte læsning af tif-filer.