Čtete formát souboru BSQ v IDL?

Používáte IDL s ENVI nebo samostatně?

Pokud používáte IDL s ENVI (zdaleka nejjednodušší, pokud máte ENVI), pak budou fungovat standardní příkazy pro čtení souborů ENVI. Například:

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 

Pokud nemáte ENVI, je to o něco těžší. Data BSQ můžete číst pomocí standardních příkazů pro čtení binárních souborů IDL (viz http://www.msi.umn.edu/software/idl/tutorial/idl-rwd.html#Reading%20Binary%20Data , ale budete muset znát rozměry obrazu (vzorky, čáry a pruhy), aby je bylo možné rozumně rozdělit. Tyto informace se často ukládají do souboru .hdr, který ENVI čte automaticky, když použijete výše uvedenou metodu.

IDL čte strukturu binárního souboru Band Sequential (BSQ) jako [Column, Row, Band]

 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. 

Takže, jak uvádí @robintw, potřebujete znát rozměry a zdroj / snímač vašeho obrazu. Kromě kódování BSQ musíte znát čtyři věci, abyste mohli načíst data do Interaktivní pole IDL. Potřebujete počet sloupců, počet řádků, počet pásem a datový typ. Samozřejmě pro každou geoprostorovou práci mimo IDL byste potřebovali znát nominální velikost pixelu a souřadnice alespoň 3 kontrolních bodů pro registraci obrázků . A další podrobnosti k str opulujte soubor ESRI Grid World.

Pro dekódování dat BSQ jsou volby binárních datových typů IDL

 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 

Takže s minimem informací po ruce: Chcete-li „otevřít“ data, vytvořte nový postup IDL

 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 

Komentáře

• @ Yudine, měla by tam být nějaká příloha nebo nějaké podrobnosti?
• @ Stuart, toto je detail souboru ENVI description = {Resize Resize Result, x resize factor: 1.000000, y resize factor: 1.000000. [Čt 13. ledna 18:15:18 2011]} ukázky = 338 řádků = 308 pásem = 125 posunutí záhlaví = 0 typ souboru = ENVI Standardní datový typ = 2 prokládané = bsq typ senzoru = Neznámé pořadí bajtů = 0 x start = 459 y start = 1735 Počet pásem = 125
• @Yudine, vypadá to, že máte všechno, co potřebujete k načtení obrázku do pole IDL. 338 vzorků – > sloupce, 308 řádků – > řádky, 125 pásem. Celé číslo datového typu 2 (pravděpodobně 0-255). Pořadí bajtů je 0, tj. Little Endian, takže přímé čtení do počítače se systémem Windows (bude muset transformovat (Swap_If_Big_Endian), pokud je 1. Takže vytvořte nový Pro v grafickém uživatelském rozhraní IDLDE a postupujte podle pseudokódu. Přečtení Michaela Todda ' s výše uvedená poznámka s odkazem na cvičení třídy v UIdaho by vám měla pomoci se syntaxí kódu. Jediným problémem pak bude rozhodování o tom, která pásma v poli se mají zobrazit, 125 kombinovat.

Nemám zkušenosti s bsq ani idl, ale protože gdal podporuje bsq , můžete zkusit použít gdal_translate k převodu souboru na geotiff. Zdá se, že IDL podporuje čtení souborů tif.