Sea Surface Temperature Fronts

Tutorial 2 : Dods to ArcInfo Tutorial
The second part of this tutorial is to walk you through the steps to take the data accessed in Part I and create coverages and grids in ArcInfo. At this point I'm assuming you have installed and used the Perl scripts to access and retrieve your data. I'm also assuming that you have the AML installed and that it is within your AMLPATH.
Directory Stuff
I have the AML set up to run with a particular file structure. The aml will read in an infile which contains pointers to the raw data in GENERATE form, the attribute datain INFO form, and the name of the default coverage. I keep the GENERATE and INFO data files in a raw_data file. The line coverages get placed in a coverage directory, and the grids get placed in a grids directory. These need to be in the same file level, e.g.
|
• |
d:\temp\coverages |
| • |
d:\temp\grids |
|
• |
d:\temp\raw_data |
At some point in the future I will alter the AML to route you to any directory you wish, but right now it is easiest if you navigate to the coverages directory within Arc.
Projection Stuff
Right now the AML creates files in two different projections: Geographic and Lambert Conic Conformal. The Lambert has the following specifics:
| COORDINATE SYSTEM DESCRIPTION |
|
|
|
|
|
|
| Projection |
LAMBERT |
|
|
|
|
| Datum |
NAD83 |
|
|
|
|
| Units |
METERS |
Spheroid |
|
|
|
| GRS1980 |
|
|
|
|
|
| Parameters: |
|
|
|
|
|
| 1st standard parallel |
|
|
36 |
9 |
36.000 |
| 2nd standard parallel |
|
|
43 |
49 |
48.000 |
| central meridian |
|
|
- 70 |
0 |
0.000 |
latitude of projection's
origin |
|
|
40 |
0 |
0.000 |
| false easting (meters) |
|
|
|
|
0.00000 |
| false northing (meters) |
|
|
|
|
0.00000 |
|
|
|
|
|
|
Runtime
If you have all of the above in working order and are in the coverages directory, then you can go ahead and fire up the AML by typing the following from an Arc prompt:
Arc: &run dodslnlp
This will fire up the aml, which will prompt you for some imputs. The first thing the aml asks for is the infile (amlfile.txt). This infile will contain references to the raw generate data, the raw attribute data, and the name of the output coverages:

Now you just sit back (hopefully!), and let the AML process through the data.
The AML will cycle through the dates in the amlfile one at a time, create coverages for each date, convert those coverages to grids, and then create "distance to grids." The ArcInfo messages for one day look like:

Now for the sample list of days you should create 6 coverages (three unique DODS days, but two different projections [see above for projection info]), and 12 grids (a grid of the line, and a distance to grid, both with two projections).
Output
A coverage file should look like this:

The grid file should look like this:

The Distance to Grid should look like this:

Misc
Congratulations. You have successfully made it through the tutorial! If you have questions or comments or a list of suggested improvements, please email me. I'd like to hear how you are using the script. A few changes that will be made include tightening up the aml to make it more flexible around directory structure. Also, I wrote the projection part of the AML to match a custom projection I use here. The AML can be easily modified for any particular projection. This query returns all of the data in the dataset. To subset out areas of specific interest to you, the AML can be altered using SETWINDOW or SETMASK to fit your custom x,y parameters.
|