[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20011204: convert AREA files to HDF files

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: 20011204: convert AREA files to HDF files
  • Date: Tue, 04 Dec 2001 17:30:42 +0000
  • >from: "Miguel A. Martinez" <address@hidden>
  • >keywords: 3C0D0842.AC148617 McIDAS AREA IDL
  • >organization: inm.es
> Sorry if this is off topic for this list but I need to find
> a program or script that will convert McIDAS AREA files
> into HDF format.  Specifically I need to convert GOES-8
> images in AREA format into HDF format for use in IDL/ENVI.
> 
> Does anyone know how to do this?
> 

If you need only a fast way to do this (for example in order to
apply
an IDL command to the image) you can made this procedure that I use for
import the GOES images to PV-WAVE (similar to IDL).

a) Run the McIDAS command AXFORM. With this command you generate one
   bynary or ASCII file with the data.

Example to create bynary files with all the 18-band of GOES
sounder with the radiances:

IMGCOPY GOES_SERVER/IR.0 LOCAL_AREA/MIAREA.1
AXFORM number_area name_file BAND=ALL UNIT=RAD FORMAT=I$

this creates the files: 
$HOME/mcidas/data/name_file.001
$HOME/mcidas/data/name_file.018                 


b) Import to IDL with the read binary command (or ASCII)

Example of procedure for import to PV-WAVE (run
lee_gvar,name_file,mi_gvar to obtain the data):
                
PRO lee_gvar,nombre,gvar


gvar=fltarr(324,189,18)


for i=1,18 do begin 
;  format i3.3 writes 001 002 etc.) 
    fichero=string(i,format='(i3.3)')
    fichero=nombre+'.'+fichero
aux=lonarr(324,189)
        openr,unit,fichero,/get_lun
         readu,unit,aux
        free_lun,unit

; divide by 1000.0 para dar radiancias (si se leen Temperaturas
; modificar y crear otro adaptado (dividir por 100. y restar 273 si 
; se quieren grados centigrados)
        gvar(*,*,i-1)=aux/1000.0
endfor
return
end

-- 
Dr. Miguel A. Martinez Rubio 
Jefe Seccion Satelites. Servicio de Teledeteccion
Instituto Nacional de Meteorologia (I.N.M.). Madrid. SPAIN
e-mail: address@hidden
Tfno: +34 91 5819 662 
Fax:  +34 91 5819 846