Re: [visad] Help script

• To: Boris Vladimir Comi <glez_b@xxxxxxxxxxxxxxxxx>
• Subject: Re: [visad] Help script
• From: Curtis Rueden <ctrueden@xxxxxxxx>
• Date: Wed, 18 Apr 2012 13:42:54 -0500

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/
>
>

• References:
  • [visad] Error Install Visad
    • From: Boris Vladimir Comi
  • Re: [visad] Error Install Visad
    • From: Tom Whittaker
  • Re: [visad] Error Install Visad
    • From: Boris Vladimir Comi
  • Re: [visad] Error Install Visad
    • From: Tom Whittaker
  • [visad] Help script
    • From: Boris Vladimir Comi