Hi Boris,
Perhaps your question would be best directed to the McIDAS-V forums?
http://dcdbs.ssec.wisc.edu/mcidasv/forums/
-Curtis
On Wed, Apr 18, 2012 at 1:04 PM, Boris Vladimir Comi <
glez_b@xxxxxxxxxxxxxxxxx> wrote:
>
> Hello, hope you are fine.
>
> I want to make a script that calculates certain characteristics in a
> satellite image for the day April 3, 2004. The data is downloaded from
> Unidata in the following format: 200404031215.goesi12r and are hourly.
>
> The first thing I want to do:
>
> 1. how to select the Region of Mexico from a script?
>
> Doing some deployments in McIDAS-V and after having researched the types
> of coordinates in McIDAS-V , I note that the values of lines and
> elements which correspond to Region of Mexico are the values determined
> by the "AREA Coordinates" (see 3abril.jpg in attached),they are the values
> associated with the image displayed on screen ?
>
> if that's true, How to convert lines and elements of an image to
> coordinates of latitude and longitude?
>
> There is some method or function to select the domain mentioned?
>
> 2. There is some method or function to display this data on screen of
> McIDAS-V?
>
> 3. I want to display a image with the values of brightness, where the
> range is from 201 to 207 and segment the region corresponding to these
> values. Is there any way to do this?
>
> The script is as follows:
>
> def scm(directory):
> import os;
> fs = os.listdir(directory);
>
> # now fs will be a list of all the files in directory
>
> from edu.wisc.ssec.mcidas import AreaFile;
> for name in fs:
>
> print "Reading in:",name
> af = AreaFile(directory+"/"+name);
> ad = af.getAreaDirectory();
> count = 0;
> data = af.getFloatData();
>
>
> # Mesoscale Convective System detected
>
> for i in xrange(ad.getLines()):
> for j in xrange(ad.getElements()):
> if 552 < i < 900 and 244 < j < 572:
> if (data[0][i][j]) > 206 and (data[0][i][j]) < 208:
> print i, j, data[0][i][j];
>
>
> scm("/home/mcidas/Documents/datos/IR")
>
>
>
> ****Boris Vladimir Comi Gonzalez
> Universidad Nacional Autónoma de México
> Grupo de Tormentas Convecivas
>
>
>
>
>
> _______________________________________________
> visad mailing list
> visad@xxxxxxxxxxxxxxxx
> For list information, to unsubscribe, visit:
> http://www.unidata.ucar.edu/mailing_lists/
>
>
