On 24 =CD=C1=D2, 23:11, Jean H <jghas...@[EMAIL PROTECTED]
>
wrote:=
> negra wrote:
> > =9A Hello all!
>
> > Could anybody help me withspatialsubsetting one image by another in
> > idl. It's easy in Envi, but I must process a large number of files.
>
> > Thanks!
> > Gulshat
> > *** Sent From/Enviado desde: =9A =9Ahttp://groups.expo.st***
>
> Hi,
> you can follow these steps in IDL:
>
> get access to ENVI's functions: envi, /restore_base_save_files
> =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A =9A
envi_batch=
_init
> Open your mask image: ENVI_OPEN_FILE
> Get the map information: =9AENVI_GET_MAP_INFO
> Extract the corners coordinates from the map info.
> Open your image tosubset
> get the data of the image tosubset: read_binary
> change the mask coordinates to the image coordinates (geographical
> coords to column/row coordinates): ENVI_CONVERT_FILE_COORDINATES
> Nowsubsetthe image: sub =3D image[x[0]:x[1],y[0]:y[1]]
> Save the data: ENVI_WRITE_ENVI_FILE =9A (you might want to do a
> envi_file_query on the mask image AND on the original image to get the
> proper header values)
>
> Jean
Thanks, I shall try to realize it. But I don't know still how to
> Extract the corners coordinates from the map info.
Gulshat
|
