#include <windows.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <fcntl.h>
#include <io.h>
#include <sys\stat.h>
#include <winlib/winlib.h>
#include <geolib/cproj.h>
#include <geolib/proj.h>

#include "mapplot.h"
#include "wdb-ii.h"
#include "map.h"


BOOL SetupOtherWorld(VARS *v, HDC hdc)
{
       RECT rect;
       POINT origPt;
       int i;
       double outParm[15];
       long iFlag;
       double lonRad, latRad, x, x2, y;
       double extX, extY;

       if (v->minMax.minLat < -89.5)
               v->minMax.minLat = -89.5;
       if (v->minMax.maxLat > 89.5)
               v->minMax.maxLat = 89.5;

       /* Report errors */
       init(0, -1, NULL, NULL);

       for(i = 0; i < 15; i++)
               outParm[i] = 0.;
/*
       Mercator:

   0, 1        Semi-major/semi-minor axis of ellipsoid
       4               Longitude of the central meridian
       5               Latitude of true scale
   6           False easting in the same units as the semi-major axis
   7           False northing in the same units as the semi-major axis

       Robinson, Hammer:

       0               Sphere
       4               Central meridian
   6           False easting in the same units as the semi-major axis
   7           False northing in the same units as the semi-major axis

       Orthographic;

       0               Sphere
       4               Central longitude
       5               Central latitude
   6           False easting in the same units as the semi-major axis
   7           False northing in the same units as the semi-major axis
*/
       for_init(v->project, 0, outParm, 0, NULL, NULL, &iFlag, v->ForTrans);
       if (iFlag)
               return FALSE;
       inv_init(v->project, 0, outParm, 0, NULL, NULL, &iFlag, v->InvTrans);
       if (iFlag)
               return FALSE;

       switch (v->project) {
               case MERCAT:
                       lonRad = 180. * D2R;
                       latRad = 89.5 * D2R;
                       if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
                               return FALSE;
                       break;

               case ROBIN:
               case HAMMER:
                       lonRad = 180. * D2R;
                       latRad = 0.;
                       if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
                               return FALSE;
                       lonRad = 0.;
                       latRad = 89.5 * D2R;
                       if (v->ForTrans[v->project](lonRad, latRad, &x2, &y))
                               return FALSE;
                       break;

               case ORTHO:
                       lonRad = 90. * D2R;
                       latRad = 0.;
                       if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
                               return FALSE;
                       lonRad = 0. * D2R;
                       latRad = 89.5 * D2R;
                       if (v->ForTrans[v->project](lonRad, latRad, &x2, &y))
                               return FALSE;
                       break;
       }

       extX = (y > x) ? 32767. * (x / y) : 32767.;
       extY = (y > x) ? 32767. : 32767. * (y / x);
       v->scale = 32767. / max(x, y);

       GetClientRect(v->hwndCont, &rect);
       SetMapMode(hdc, MM_ISOTROPIC);
       SetWindowExtEx(hdc, (int)extX, (int)extY, NULL);
       SetViewportExtEx(hdc, (rect.right / 2) * v->zoom,
                       (-rect.bottom / 2) * v->zoom, NULL);

       lonRad = v->origLon * D2R;
       latRad = v->origLat * D2R;
       if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
               return FALSE;

       origPt.x = (LONG)(x * v->scale);
       origPt.y = (LONG)(y * v->scale);
       LPtoDP(hdc, &origPt, 1);
       SetViewportOrgEx(hdc, rect.right / 2 - origPt.x,
                       rect.bottom / 2 - origPt.y, NULL);

       return TRUE;
}


BOOL SetupWorld(VARS *v, HDC hdc)
{
       RECT rect;
       POINT origPt;

       if (v->project != PLATE_CARREE)
               return SetupOtherWorld(v, hdc);

       if (v->minMax.minLat == -89.5)
               v->minMax.minLat = -90.;
       if (v->minMax.maxLat == 89.5)
               v->minMax.maxLat = 90.;

       GetClientRect(v->hwndCont, &rect);
       SetMapMode(hdc, MM_ISOTROPIC);
       SetWindowExtEx(hdc, 18000, 9000, NULL);
       SetViewportExtEx(hdc, (rect.right / 2) * v->zoom,
                       (-rect.bottom / 2) * v->zoom, NULL);
       origPt.x = (LONG)(v->origLon * 100.);
       origPt.y = (LONG)(v->origLat * 100.);
       LPtoDP(hdc, &origPt, 1);
       SetViewportOrgEx(hdc, rect.right / 2 - origPt.x,
                       rect.bottom / 2 - origPt.y, NULL);

       return TRUE;
}


BOOL ReadMap(VARS *v, char *name)
{
       int fd, numRead, toRead;
       DPOINT p;
       MAP_POLY *newPoly;
       char tmp[1024];
       int i, np = 0;
       int maxNp = 0;
       struct stat st;
       MAP_REC *dataBuf;

       sprintf(tmp, "%s\\%s.PNT", v->dataDir, name);
       if ((fd = open(tmp, O_RDONLY | O_BINARY)) == -1) {
               MBPrintf("MapPlot", "Couldn't open data file '%s'", tmp);
               return FALSE;
       }

       stat(tmp, &st);
       dataBuf = (MAP_REC *)malloc(st.st_size);
       read(fd, dataBuf, st.st_size);
       close(fd);

       for (i = 0; i < (st.st_size / sizeof(MAP_REC)); i++) {
               if (dataBuf[i].code > 5) {
                       if (np) {
                               newPoly->nPoints = np;
                               if (! v->mapData)
                                       v->mapData = newPoly;
                               if (v->prevPoly)
                                       v->prevPoly->next = newPoly;
                               v->prevPoly = newPoly;
                       }

                       newPoly = (MAP_POLY *)malloc(sizeof(MAP_POLY));

                       if (dataBuf[i].code < 2001)
                               newPoly->feature = 1;
                       else if (dataBuf[i].code < 4001)
                               newPoly->feature = 2;
                       else if (dataBuf[i].code < 5001)
                               newPoly->feature = 4;
                       else if (dataBuf[i].code < 6001)
                               newPoly->feature = 5;
                       else if (dataBuf[i].code < 7001)
                               newPoly->feature = 6;
                       else
                               newPoly->feature = 7;

                       newPoly->points = (MAP_POINT *)malloc(sizeof(MAP_POINT));
                       newPoly->next = NULL;

                       np = 0;
                       newPoly->points[np].lon = (double)dataBuf[i].lon / 60.;
                       newPoly->points[np].lat = (double)dataBuf[i].lat / 60.;
                       newPoly->points[np].level = 100;
                       np++;
               } else {
                       newPoly->points = (MAP_POINT *)realloc(newPoly->points,
                               sizeof(MAP_POINT) * (np + 1));
                       newPoly->points[np].lon = (double)dataBuf[i].lon / 60.;
                       newPoly->points[np].lat = (double)dataBuf[i].lat / 60.;
                       newPoly->points[np].level = (char)dataBuf[i].code;
                       np++;
               }
       }

       if (np) {
               newPoly->nPoints = np;
               if (! v->mapData)
                       v->mapData = newPoly;
               if (v->prevPoly)
                       v->prevPoly->next = newPoly;
               v->prevPoly = newPoly;
       }

       close(fd);
       free(dataBuf);
       return TRUE;
}


BOOL ReadMaps(VARS *v)
{
       char *names[] = {"PCOAST", "PISLAND", "PLAKE", "RIVER", "PUSA48",
                       NULL};
       char *name;
       int i;

       for (i = 0; name = names[i]; i++) {
               if (! ReadMap(v, name))
                       return FALSE;
       }
       return TRUE;
}


BOOL DrawOtherGrid(VARS *v, HDC hdc)
{
       HFONT hFnt = NULL;
       HBRUSH hBrush;
       HPEN hPen;
       POINT oldPt;
       POINT pts[800];
       RECT rect;
       int i, gridDeg;
       double lonRad, latRad, x, x2, y;
       char tmp[30];
       double deg, deg2;

       hBrush = SelectObject(hdc, CreateSolidBrush(v->colors[CLR_GRID]));
       hPen = SelectObject(hdc, CreatePen(PS_SOLID, 0, v->colors[CLR_GRID]));

       for (deg = v->minMax.minLon; deg <= v->minMax.maxLon;
                       deg += (double)v->gridDeg) {

               for (deg2 = v->minMax.minLat, i = 0; deg2 <= v->minMax.maxLat;
                               deg2++, i++) {
                       lonRad = deg * D2R;
                       latRad = deg2 * D2R;
                       if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
                               goto error;
                       pts[i].x = (int)(x * v->scale);
                       pts[i].y = (int)(y * v->scale);
               }
               Polyline(hdc, pts, i);
       }

       for (deg = v->minMax.minLat; deg <= ((v->minMax.maxLat == 89.5) ? 90. : v->minMax.maxLat);
                       deg += (deg == -89.5) ? (double)v->gridDeg - .5: (double)v->gridDeg) {

               for (deg2 = v->minMax.minLon, i = 0; deg2 <= v->minMax.maxLon;
                               deg2++, i++) {
                       lonRad = deg2 * D2R;
                       latRad = (deg == 90.) ? 89.5 * D2R: deg * D2R;
                       if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
                               goto error;
                       pts[i].x = (int)(x * v->scale);
                       pts[i].y = (int)(y * v->scale);
               }
               Polyline(hdc, pts, i);
       }

       DeleteObject(SelectObject(hdc, hBrush));
       DeleteObject(SelectObject(hdc, hPen));

       if (! (v->shwFeatures & SHW_DEGREES))
               return TRUE;

       GetClientRect(v->hwndCont, (LPRECT) &rect);
       DPtoLP(hdc, (LPPOINT) &rect, 2);

       hFnt = SelectObject(hdc, EzCreateFont(hdc, "Arial", 65, 0, 0, TRUE));
       hBrush = SelectObject(hdc, CreateSolidBrush(0x00000000));
       hPen = SelectObject(hdc, CreatePen(PS_SOLID, 0, 0x00000000));
       SetBkMode(hdc, TRANSPARENT);

       for (deg = v->minMax.minLon; deg <= v->minMax.maxLon; deg += (double)v->gridDeg) {

               if ((deg == v->minMax.minLon) || (deg == v->minMax.maxLon))
                       continue;

               if (v->ForTrans[v->project](deg * D2R, v->minMax.maxLat * D2R, &x, &y))
                       goto error;

               if ((double)(rect.top / v->scale) < y) {
                       if (v->InvTrans[v->project](x, (double)(rect.top / v->scale),
                                       &lonRad, &latRad))
                               goto error;

                       if (v->ForTrans[v->project](deg * D2R, latRad, &x, &y))
                               goto error;
               }

               sprintf(tmp, "%d", (int)deg);
               TextOut(hdc,
                       (int)(x * v->scale), (int)(y * v->scale), tmp, strlen(tmp));
       }

       for (deg = v->minMax.minLat; deg <= v->minMax.maxLat;
                       deg += (deg == -89.5) ? (double)v->gridDeg - .5: (double)v->gridDeg) {

               if ((deg == v->minMax.minLat) || (deg == v->minMax.maxLat))
                       continue;

               if (v->ForTrans[v->project](v->minMax.minLon * D2R, deg * D2R, &x, &y))
                       goto error;

               if ((double)(rect.left / v->scale) > x) {
                       if (v->InvTrans[v->project]((double)rect.left / v->scale, y,
                                       &lonRad, &latRad))
                               goto error;

                       if (v->ForTrans[v->project](lonRad, deg * D2R, &x, &y))
                               goto error;
               }

               sprintf(tmp, "%d", (int)deg);
               TextOut(hdc, (int)(x * v->scale), (int)(y * v->scale),
                               tmp, strlen(tmp));
       }

       DeleteObject(SelectObject(hdc, hBrush));
       DeleteObject(SelectObject(hdc, hPen));
       DeleteObject(SelectObject(hdc, hFnt));
       return TRUE;

error:
       DeleteObject(SelectObject(hdc, hBrush));
       DeleteObject(SelectObject(hdc, hPen));
       if (hFnt)
               DeleteObject(SelectObject(hdc, hFnt));
       return FALSE;
}


BOOL DrawGrid(VARS *v, HDC hdc)
{
       HFONT hFnt;
       HBRUSH hBrush;
       HPEN hPen;
       POINT oldPt;
       RECT rect;
       char tmp[30];
       double deg;

       if (v->project != PLATE_CARREE)
               return DrawOtherGrid(v, hdc);

       hBrush = SelectObject(hdc, CreateSolidBrush(v->colors[CLR_GRID]));
       hPen = SelectObject(hdc, CreatePen(PS_SOLID, 0, v->colors[CLR_GRID]));

       for (deg = v->minMax.minLon; deg <= v->minMax.maxLon; deg += (double)v->gridDeg) {
               MoveToEx(hdc, (int)(deg * 100.), (int)(v->minMax.maxLat * 100.), &oldPt);
               LineTo(hdc, (int)(deg * 100.), (int)(v->minMax.minLat * 100.));
       }

       for (deg = v->minMax.minLat; deg <= v->minMax.maxLat; deg += (double)v->gridDeg) {
               MoveToEx(hdc, (int)(v->minMax.minLon * 100.), (int)(deg * 100.), &oldPt);
               LineTo(hdc, (int)(v->minMax.maxLon * 100.), (int)(deg * 100.));
       }

       DeleteObject(SelectObject(hdc, hBrush));
       DeleteObject(SelectObject(hdc, hPen));

       if (! (v->shwFeatures & SHW_DEGREES))
               return TRUE;

       GetClientRect(v->hwndCont, (LPRECT) &rect);
       DPtoLP(hdc, (LPPOINT) &rect, 2);

       hFnt = SelectObject(hdc, EzCreateFont(hdc, "Arial", 65, 0, 0, TRUE));
       hBrush = SelectObject(hdc, CreateSolidBrush(0x00000000));
       hPen = SelectObject(hdc, CreatePen(PS_SOLID, 0, 0x00000000));
       SetBkMode(hdc, TRANSPARENT);

       for (deg = v->minMax.minLon; deg <= v->minMax.maxLon; deg += (double)v->gridDeg) {
               if ((deg == v->minMax.minLon) || (deg == v->minMax.maxLon))
                       continue;
               sprintf(tmp, "%d", (int)deg);
               TextOut(hdc, (int)(deg * 100.),
                               (rect.top <= (int)(v->minMax.maxLat * 100.)) ? rect.top :
                                               (int)(v->minMax.maxLat * 100.),
                               tmp, strlen(tmp));
       }

       for (deg = v->minMax.minLat; deg <= v->minMax.maxLat; deg += (double)v->gridDeg) {
               if ((deg == v->minMax.minLat) || (deg == v->minMax.maxLat))
                       continue;
               sprintf(tmp, "%d", (int)deg);
               TextOut(hdc,
                               (rect.left >= (int)(v->minMax.minLon * 100.)) ? rect.left :
                                               (int)(v->minMax.minLon * 100.),
                               (int)(deg * 100.), tmp, strlen(tmp));
       }

       DeleteObject(SelectObject(hdc, hBrush));
       DeleteObject(SelectObject(hdc, hPen));
       DeleteObject(SelectObject(hdc, hFnt));
       return TRUE;
}


BOOL DoDrawMap(VARS *v, HDC hdc)
{
       MAP_POLY *mapPoly;
       POINT *pts;
       RECT rect;
       HBRUSH hBrush;
       HPEN hPen = NULL;
       COLORREF color;
       int i, np;
       double lonRad, latRad, x, y;
       double lastLon = 0.;

       ShowStatus(v->hwndStatus, 0, "Processing...");

       GetClientRect(v->hwndCont, &rect);

       hBrush = CreateSolidBrush(v->colors[CLR_BACKGR]);
       FillRect (hdc, &rect, hBrush);
       DeleteObject(hBrush);

       if (! SetupWorld(v, hdc))
               goto error;

       for (mapPoly = v->mapData; mapPoly; mapPoly = mapPoly->next) {
               switch (mapPoly->feature) {
                       case FEAT_COAST:
                               if (! (v->shwFeatures & SHW_COAST))
                                       continue;
                               color = v->colors[CLR_COAST];
                               break;

                       case FEAT_COUNTRY:
                               if (! (v->shwFeatures & SHW_COUNTRY))
                                       continue;
                               color = v->colors[CLR_COUNTRY];
                               break;

                       case FEAT_STATE:
                               if (! (v->shwFeatures & SHW_STATE))
                                       continue;
                               color = v->colors[CLR_STATE];
                               break;

                       case FEAT_ISLAND:
                               if (! (v->shwFeatures & SHW_ISLAND))
                                       continue;
                               color = v->colors[CLR_ISLAND];
                               break;

                       case FEAT_LAKE:
                               if (! (v->shwFeatures & SHW_LAKE))
                                       continue;
                               color = v->colors[CLR_LAKE];
                               break;

                       case FEAT_RIVER:
                               if (! (v->shwFeatures & SHW_RIVER))
                                       continue;
                               color = v->colors[CLR_RIVER];
                               break;

                       default:
                               exit(1);
                               continue;
               }

               pts = (POINT *)malloc(sizeof(POINT) * (int)mapPoly->nPoints);

               hBrush = SelectObject(hdc, CreateSolidBrush(color));
               hPen = SelectObject(hdc, CreatePen(PS_SOLID, 0, color));

               np = 0;
               for (i = 0; i < (int)mapPoly->nPoints; i++) {

                       if (mapPoly->points[i].level < 1)
                               continue;

                       if ((mapPoly->points[i].lon > v->minMax.maxLon) ||
                                       (mapPoly->points[i].lon < v->minMax.minLon) ||
                                       (mapPoly->points[i].lat > v->minMax.maxLat) ||
                                       (mapPoly->points[i].lat < v->minMax.minLat)) {
                               if (np) {
                                       if (np == 1)
                                               SetPixel(hdc, pts[0].x, pts[0].y, color);
                                       else
                                               Polyline(hdc, pts, np);
                                       np = 0;
                               }
                               continue;
                       }

                       if (fabs(lastLon - mapPoly->points[i].lon) > 180.) {
                               if (np) {
                                       if (np == 1)
                                               SetPixel(hdc, pts[0].x, pts[0].y, color);
                                       else
                                               Polyline(hdc, pts, np);
                                       np = 0;
                               }
                       }

                       if (np == 16000) {
                               Polyline(hdc, pts, 16000);
                               np = 0;
                       }

                       if (v->project == PLATE_CARREE) {
                               pts[np].x = (int)(mapPoly->points[i].lon * 100.);
                               pts[np].y = (int)(mapPoly->points[i].lat * 100.);
                       } else {
                               lonRad = mapPoly->points[i].lon * D2R;
                               latRad = mapPoly->points[i].lat * D2R;
                               if (v->ForTrans[v->project](lonRad, latRad, &x, &y))
                                       goto error;
                               pts[np].x = (int)(x * v->scale);
                               pts[np].y = (int)(y * v->scale);
                       }

                       np++;
                       lastLon = mapPoly->points[i].lon;
               }

               if (np) {
                       if (np == 1)
                               SetPixel(hdc, pts[0].x, pts[0].y, color);
                       else
                               Polyline(hdc, pts, np);
               }

               DeleteObject(SelectObject(hdc, hBrush));
               DeleteObject(SelectObject(hdc, hPen));
               hPen = NULL;
               free(pts);
               pts = NULL;
       }

       if (v->shwFeatures & SHW_GRID) {
               if (! DrawGrid(v, hdc))
                       goto error;
       }

       ShowStatus(v->hwndStatus, 0, "Done");
       return TRUE;

error:
       if (pts)
               free(pts);
       if (hPen) {
               DeleteObject(SelectObject(hdc, hBrush));
               DeleteObject(SelectObject(hdc, hPen));
       }
       ShowStatus(v->hwndStatus, 0, "");
       return FALSE;
}


void DrawMap(VARS *v, HDC hdc)
{
       if (v->wdbIIMode)
               DrawWDBIIMaps(v, hdc);
       else
               DoDrawMap(v, hdc);
}

You are viewing proxied material from bitreich.org. The copyright of proxied material belongs to its original authors. Any comments or complaints in relation to proxied material should be directed to the original authors of the content concerned. Please see the disclaimer for more details.