-/****************************************************************************
-* SPLAT: An RF Signal Propagation Loss and Terrain Analysis Tool *
-* Last update: 15-Mar-2007 *
-*****************************************************************************
-* Project started in 1997 by John A. Magliacane, KD2BD *
-*****************************************************************************
-* *
-* Extensively modified by J. D. McDonald in Jan. 2004 to include *
-* the Longley-Rice propagation model using C++ code from NTIA/ITS. *
-* *
-* See: http://flattop.its.bldrdoc.gov/itm.html *
-* *
-*****************************************************************************
-* *
-* This program is free software; you can redistribute it and/or modify it *
-* under the terms of the GNU General Public License as published by the *
-* Free Software Foundation; either version 2 of the License or any later *
-* version. *
-* *
-* This program is distributed in the hope that it will useful, but WITHOUT *
-* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or *
-* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License *
-* for more details. *
-* *
-*****************************************************************************
- g++ -Wall -O3 -s -lm -lbz2 -fomit-frame-pointer itm.cpp splat.cpp -o splat
-*****************************************************************************/
+/****************************************************************************\
+* SPLAT!: An RF Signal Path Loss And Terrain Analysis Tool *
+******************************************************************************
+* Project started in 1997 by John A. Magliacane, KD2BD *
+* Last update: 19-Oct-2007 *
+******************************************************************************
+* Please consult the documentation for a complete list of *
+* individuals who have contributed to this project. *
+******************************************************************************
+* *
+* This program is free software; you can redistribute it and/or modify it *
+* under the terms of the GNU General Public License as published by the *
+* Free Software Foundation; either version 2 of the License or any later *
+* version. *
+* *
+* This program is distributed in the hope that it will useful, but WITHOUT *
+* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or *
+* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License *
+* for more details. *
+* *
+******************************************************************************
+* g++ -Wall -O3 -s -lm -lbz2 -fomit-frame-pointer itm.cpp splat.cpp -o splat *
+\****************************************************************************/
#include <stdio.h>
#include <math.h>
#include <bzlib.h>
#include <unistd.h>
#include "fontdata.h"
-#include "smallfont.h"
#define GAMMA 2.5
-#define MAXSLOTS 9
+#define MAXPAGES 9
#define BZBUFFER 65536
-#if MAXSLOTS==4
+#if MAXPAGES==4
#define ARRAYSIZE 4950
#endif
-#if MAXSLOTS==9
+#if MAXPAGES==9
#define ARRAYSIZE 10870
#endif
-#if MAXSLOTS==16
+#if MAXPAGES==16
#define ARRAYSIZE 19240
#endif
-#if MAXSLOTS==25
+#if MAXPAGES==25
#define ARRAYSIZE 30025
#endif
-char string[255], sdf_path[255], opened=0, *splat_version={"1.2.0b"};
+char string[255], sdf_path[255], opened=0, *splat_version={"1.2.1"};
double TWOPI=6.283185307179586, HALFPI=1.570796326794896,
PI=3.141592653589793, deg2rad=1.74532925199e-02,
EARTHRADIUS=20902230.97, METERS_PER_MILE=1609.344,
METERS_PER_FOOT=0.3048, KM_PER_MILE=1.609344, earthradius,
- max_range=0.0;
+ max_range=0.0, forced_erp=-1.0, fzone_clearance=0.6;
int min_north=90, max_north=-90, min_west=360, max_west=-1,
max_elevation=-32768, min_elevation=32768, bzerror, maxdB=230;
-unsigned char got_elevation_pattern=0, got_azimuth_pattern=0, metric=0;
+unsigned char got_elevation_pattern, got_azimuth_pattern, metric=0;
struct site { double lat;
double lon;
float alt;
char name[50];
+ char filename[255];
} site;
struct path { double lat[ARRAYSIZE];
int min_el;
short data[1200][1200];
unsigned char mask[1200][1200];
- } dem[MAXSLOTS];
+ unsigned char signal[1200][1200];
+ } dem[MAXPAGES];
struct LR { double eps_dielect;
double sgm_conductivity;
double frq_mhz;
double conf;
double rel;
+ double erp;
int radio_climate;
int pol;
float antenna_pattern[361][1001];
} LR;
+struct region { unsigned char color[32][3];
+ int level[32];
+ int levels;
+ } region;
+
double elev_l[ARRAYSIZE+10];
void point_to_point(double elev[], double tht_m, double rht_m,
return (string);
}
+int PutMask(double lat, double lon, int value)
+{
+ /* Lines, text, markings, and coverage areas are stored in a
+ mask that is combined with topology data when topographic
+ maps are generated by SPLAT!. This function sets and resets
+ bits in the mask based on the latitude and longitude of the
+ area pointed to. */
+
+ int x, y, indx;
+ char found;
+
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<=(double)dem[indx].max_north && lon>=(double)dem[indx].min_west && lon<=(double)dem[indx].max_west)
+ found=1;
+ else
+ indx++;
+
+ if (found)
+ {
+ x=(int)(1199.0*(lat-floor(lat)));
+ y=(int)(1199.0*(lon-floor(lon)));
+
+ dem[indx].mask[x][y]=value;
+
+ return (dem[indx].mask[x][y]);
+ }
+
+ else
+ return -1;
+}
+
int OrMask(double lat, double lon, int value)
{
/* Lines, text, markings, and coverage areas are stored in a
the mask based on the latitude and longitude of the area
pointed to. */
- int x, y, indx, minlat, minlon;
+ int x, y, indx;
char found;
- minlat=(int)floor(lat);
- minlon=(int)floor(lon);
-
- for (indx=0, found=0; indx<MAXSLOTS && found==0;)
- if (minlat==dem[indx].min_north && minlon==dem[indx].min_west)
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<=(double)dem[indx].max_north && lon>=(double)dem[indx].min_west && lon<=(double)dem[indx].max_west)
found=1;
else
indx++;
return (OrMask(lat,lon,0));
}
+int PutSignal(double lat, double lon, unsigned char signal)
+{
+ /* This function writes a signal level (0-255)
+ at the specified location for later recall. */
+
+ int x, y, indx;
+ char found;
+
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<=(double)dem[indx].max_north && lon>=(double)dem[indx].min_west && lon<=(double)dem[indx].max_west)
+ found=1;
+ else
+ indx++;
+
+ if (found)
+ {
+ x=(int)(1199.0*(lat-floor(lat)));
+ y=(int)(1199.0*(lon-floor(lon)));
+
+ dem[indx].signal[x][y]=signal;
+
+ return (dem[indx].signal[x][y]);
+ }
+
+ else
+ return 0;
+}
+
+unsigned char GetSignal(double lat, double lon)
+{
+ /* This function reads the signal level (0-255) at the
+ specified location that was previously written by the
+ complimentary PutSignal() function. */
+
+ int x, y, indx;
+ char found;
+
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<=(double)dem[indx].max_north && lon>=(double)dem[indx].min_west && lon<=(double)dem[indx].max_west)
+ found=1;
+ else
+ indx++;
+
+ if (found)
+ {
+ x=(int)(1199.0*(lat-floor(lat)));
+ y=(int)(1199.0*(lon-floor(lon)));
+
+ return (dem[indx].signal[x][y]);
+ }
+
+ else
+ return 0;
+}
+
double GetElevation(struct site location)
{
/* This function returns the elevation (in feet) of any location
Function returns -5000.0 for locations not found in memory. */
char found;
- int x, y, indx, minlat, minlon;
+ int x, y, indx;
double elevation;
elevation=-5000.0;
- minlat=(int)floor(location.lat);
- minlon=(int)floor(location.lon);
-
x=(int)(1199.0*(location.lat-floor(location.lat)));
y=(int)(1199.0*(location.lon-floor(location.lon)));
- for (indx=0, found=0; indx<MAXSLOTS && found==0; indx++)
+ for (indx=0, found=0; indx<MAXPAGES && found==0; indx++)
{
- if (minlat==dem[indx].min_north && minlon==dem[indx].min_west)
+ if (location.lat>=(double)dem[indx].min_north && location.lat<=(double)dem[indx].max_north && location.lon>=(double)dem[indx].min_west && location.lon<=(double)dem[indx].max_west)
{
-
elevation=3.28084*dem[indx].data[x][y];
found=1;
}
not found in memory. */
char found;
- int x, y, indx, minlat, minlon;
-
- minlat=(int)floor(lat);
- minlon=(int)floor(lon);
+ int x, y, indx;
x=(int)(1199.0*(lat-floor(lat)));
y=(int)(1199.0*(lon-floor(lon)));
- for (indx=0, found=0; indx<MAXSLOTS && found==0; indx++)
+ for (indx=0, found=0; indx<MAXPAGES && found==0; indx++)
{
- if (minlat==dem[indx].min_north && minlon==dem[indx].min_west)
+ if (lat>=(double)dem[indx].min_north && lat<=(double)dem[indx].max_north && lon>=(double)dem[indx].min_west && lon<=(double)dem[indx].max_west)
{
-
dem[indx].data[x][y]+=(short)rint(height);
found=1;
}
increment=total_distance/path_length; /* Miles per sample */
- for (distance=0, c=0; distance<=total_distance; distance+=increment)
+ for (distance=0, c=0; (distance<=total_distance && c<ARRAYSIZE); distance+=increment, c++)
{
beta=distance/3959.0;
lat2=asin(sin(lat1)*cos(beta)+cos(azimuth)*sin(beta)*cos(lat1));
lat2=lat2/deg2rad;
lon2=lon2/deg2rad;
- if (c<ARRAYSIZE)
- {
- path.lat[c]=lat2;
- path.lon[c]=lon2;
- tempsite.lat=lat2;
- tempsite.lon=lon2;
- path.elevation[c]=GetElevation(tempsite);
- path.distance[c]=distance;
- c++;
- }
+ path.lat[c]=lat2;
+ path.lon[c]=lon2;
+ tempsite.lat=lat2;
+ tempsite.lon=lon2;
+ path.elevation[c]=GetElevation(tempsite);
+ path.distance[c]=distance;
}
/* Make sure exact destination point is recorded at path.length-1 */
/* Anything else returns a 0.0 */
- if (bearing>360.0 || bearing<-90.0)
+ if (bearing>360.0 || bearing<-360.0)
bearing=0.0;
return bearing;
tempsite.lon=361.0;
tempsite.alt=0.0;
tempsite.name[0]=0;
+ tempsite.filename[0]=0;
fd=fopen(qthfile,"r");
fgets(string,49,fd);
tempsite.lon=ReadBearing(string);
+ if (tempsite.lon<0.0)
+ tempsite.lon+=360.0;
+
/* Antenna Height */
fgets(string,49,fd);
fclose(fd);
must be converted to feet before exiting. */
for (x=0; string[x]!='M' && string[x]!='m' && string[x]!=0 && x<48; x++);
+
if (string[x]=='M' || string[x]=='m')
{
string[x]=0;
string[x]=0;
sscanf(string,"%f",&tempsite.alt);
}
+
+ for (x=0; x<254 && qthfile[x]!=0; x++)
+ tempsite.filename[x]=qthfile[x];
+
+ tempsite.filename[x]=0;
}
return tempsite;
float az, xx, elevation, amplitude, rotation, valid1, valid2,
delta, azimuth[361], azimuth_pattern[361], el_pattern[10001],
elevation_pattern[361][1001], slant_angle[361], tilt,
- mechanical_tilt, tilt_azimuth, tilt_increment, sum;
+ mechanical_tilt=0.0, tilt_azimuth, tilt_increment, sum;
FILE *fd=NULL;
unsigned char read_count[10001];
rotation=0.0;
+ got_azimuth_pattern=0;
+ got_elevation_pattern=0;
+
/* Load .az antenna pattern file */
fd=fopen(azfile,"r");
dem[] structure. */
int x, y, data, indx, minlat, minlon, maxlat, maxlon;
- char found, free_slot=0, line[20], sdf_file[255],
+ char found, free_page=0, line[20], sdf_file[255],
path_plus_name[255];
FILE *fd;
/* Is it already in memory? */
- for (indx=0, found=0; indx<MAXSLOTS && found==0; indx++)
+ for (indx=0, found=0; indx<MAXPAGES && found==0; indx++)
{
if (minlat==dem[indx].min_north && minlon==dem[indx].min_west && maxlat==dem[indx].max_north && maxlon==dem[indx].max_west)
found=1;
if (found==0)
{
- for (indx=0, free_slot=0; indx<MAXSLOTS && free_slot==0; indx++)
+ for (indx=0, free_page=0; indx<MAXPAGES && free_page==0; indx++)
if (dem[indx].max_north==-90)
- free_slot=1;
+ free_page=1;
}
indx--;
- if (free_slot && found==0 && indx>=0 && indx<MAXSLOTS)
+ if (free_page && found==0 && indx>=0 && indx<MAXPAGES)
{
/* Search for SDF file in current working directory first */
if (fd!=NULL)
{
- fprintf(stdout,"Loading \"%s\" into slot %d...",path_plus_name,indx+1);
+ fprintf(stdout,"Loading \"%s\" into page %d...",path_plus_name,indx+1);
fflush(stdout);
fgets(line,19,fd);
for (y=0; y<1200; y++)
{
fgets(line,19,fd);
- sscanf(line,"%d",&data);
+ data=atoi(line);
dem[indx].data[x][y]=data;
+ dem[indx].signal[x][y]=0;
+ dem[indx].mask[x][y]=0;
if (data>dem[indx].max_el)
dem[indx].max_el=data;
stored in the first available dem[] structure. */
int x, y, data, indx, minlat, minlon, maxlat, maxlon;
- char found, free_slot=0, sdf_file[255], path_plus_name[255];
+ char found, free_page=0, sdf_file[255], path_plus_name[255], *string;
FILE *fd;
BZFILE *bzfd;
/* Is it already in memory? */
- for (indx=0, found=0; indx<MAXSLOTS && found==0; indx++)
+ for (indx=0, found=0; indx<MAXPAGES && found==0; indx++)
{
if (minlat==dem[indx].min_north && minlon==dem[indx].min_west && maxlat==dem[indx].max_north && maxlon==dem[indx].max_west)
found=1;
if (found==0)
{
- for (indx=0, free_slot=0; indx<MAXSLOTS && free_slot==0; indx++)
+ for (indx=0, free_page=0; indx<MAXPAGES && free_page==0; indx++)
if (dem[indx].max_north==-90)
- free_slot=1;
+ free_page=1;
}
indx--;
- if (free_slot && found==0 && indx>=0 && indx<MAXSLOTS)
+ if (free_page && found==0 && indx>=0 && indx<MAXPAGES)
{
/* Search for SDF file in current working directory first */
if (fd!=NULL && bzerror==BZ_OK)
{
- fprintf(stdout,"Loading \"%s\" into slot %d...",path_plus_name,indx+1);
+ fprintf(stdout,"Loading \"%s\" into page %d...",path_plus_name,indx+1);
fflush(stdout);
sscanf(BZfgets(bzfd,255),"%d",&dem[indx].max_west);
for (x=0; x<1200; x++)
for (y=0; y<1200; y++)
{
- sscanf(BZfgets(bzfd,20),"%d",&data);
+ string=BZfgets(bzfd,20);
+ data=atoi(string);
dem[indx].data[x][y]=data;
+ dem[indx].signal[x][y]=0;
+ dem[indx].mask[x][y]=0;
if (data>dem[indx].max_el)
dem[indx].max_el=data;
requested must be entirely over water. */
int x, y, indx, minlat, minlon, maxlat, maxlon;
- char found, free_slot=0;
+ char found, free_page=0;
int return_value=-1;
/* Try to load an uncompressed SDF first. */
/* Is it already in memory? */
- for (indx=0, found=0; indx<MAXSLOTS && found==0; indx++)
+ for (indx=0, found=0; indx<MAXPAGES && found==0; indx++)
{
if (minlat==dem[indx].min_north && minlon==dem[indx].min_west && maxlat==dem[indx].max_north && maxlon==dem[indx].max_west)
found=1;
if (found==0)
{
- for (indx=0, free_slot=0; indx<MAXSLOTS && free_slot==0; indx++)
+ for (indx=0, free_page=0; indx<MAXPAGES && free_page==0; indx++)
if (dem[indx].max_north==-90)
- free_slot=1;
+ free_page=1;
}
indx--;
- if (free_slot && found==0 && indx>=0 && indx<MAXSLOTS)
+ if (free_page && found==0 && indx>=0 && indx<MAXPAGES)
{
- fprintf(stdout,"Region \"%s\" assumed as sea-level into slot %d...",name,indx+1);
+ fprintf(stdout,"Region \"%s\" assumed as sea-level into page %d...",name,indx+1);
fflush(stdout);
dem[indx].max_west=maxlon;
for (y=0; y<1200; y++)
{
dem[indx].data[x][y]=0;
+ dem[indx].signal[x][y]=0;
+ dem[indx].mask[x][y]=0;
if (dem[indx].min_el>0)
dem[indx].min_el=0;
{
fgets(input,78,fd);
- fprintf(stdout,"Reading \"%s\"... ",filename);
+ fprintf(stdout,"\nReading \"%s\"... ",filename);
fflush(stdout);
while (fd!=NULL && feof(fd)==0)
city_site.lon=ReadBearing(str[2]);
city_site.alt=0.0;
+ if (city_site.lon<0.0)
+ city_site.lon+=360.0;
+
PlaceMarker(city_site);
fgets(input,78,fd);
}
fclose(fd);
- fprintf(stdout,"Done!\n");
+ fprintf(stdout,"Done!");
fflush(stdout);
}
if (pointer!=NULL)
*pointer=0;
- fprintf(stdout,"Reading \"%s\"... ",filename);
+ fprintf(stdout,"\nReading \"%s\"... ",filename);
fflush(stdout);
while (feof(fd1)==0)
latitude=ReadBearing(str[0]);
longitude=ReadBearing(str[1]);
+ if (longitude<0.0)
+ longitude+=360.0;
+
/* Remove <CR> and/or <LF> from antenna height string */
for (i=0; str[2][i]!=13 && str[2][i]!=10 && str[2][i]!=0; i++);
fclose(fd2);
close(fd);
- fprintf(stdout,"Done!\n");
+ fprintf(stdout,"Done!");
fflush(stdout);
fd1=fopen(tempname,"r");
}
else
- fprintf(stderr,"*** ERROR: \"%s\": not found!\n",filename);
+ fprintf(stderr,"\n*** ERROR: \"%s\": not found!",filename);
}
void LoadBoundaries(char *filename)
{
fgets(string,78,fd);
- fprintf(stdout,"Reading \"%s\"... ",filename);
+ fprintf(stdout,"\nReading \"%s\"... ",filename);
fflush(stdout);
do
fclose(fd);
- fprintf(stdout,"Done!\n");
+ fprintf(stdout,"Done!");
fflush(stdout);
}
else
- fprintf(stderr,"*** ERROR: \"%s\": not found!\n",filename);
+ fprintf(stderr,"\n*** ERROR: \"%s\": not found!",filename);
}
-void ReadLRParm(char *txsite_filename)
+char ReadLRParm(struct site txsite, char forced_read)
{
/* This function reads Longley-Rice parameter data for the
transmitter site. The file name is the same as the txsite,
except the filename extension is .lrp. If the needed file
is not found, then the file "splat.lrp" is read from the
- current working directory. Failure to load this file will
- result in the default parameters hard coded into this
- function to be used and written to "splat.lrp". */
+ current working directory. Failure to load this file under
+ a forced_read condition will result in the default parameters
+ hard coded into this function to be used and written to
+ "splat.lrp". */
double din;
- char filename[255], string[80], *pointer=NULL;
+ char filename[255], string[80], *pointer=NULL, return_value=0;
int iin, ok=0, x;
FILE *fd=NULL, *outfile=NULL;
- /* Default parameters in case things go bad */
+ /* Default parameters */
- LR.eps_dielect=15.0;
- LR.sgm_conductivity=0.005;
- LR.eno_ns_surfref=301.0;
- LR.frq_mhz=300.0;
- LR.radio_climate=5;
+ LR.eps_dielect=0.0;
+ LR.sgm_conductivity=0.0;
+ LR.eno_ns_surfref=0.0;
+ LR.frq_mhz=0.0;
+ LR.radio_climate=0;
LR.pol=0;
- LR.conf=0.50;
- LR.rel=0.50;
+ LR.conf=0.0;
+ LR.rel=0.0;
+ LR.erp=0.0;
- /* Modify txsite filename to one with a .lrp extension. */
+ /* Generate .lrp filename from txsite filename. */
- strncpy(filename,txsite_filename,255);
-
- for (x=0; filename[x]!='.' && filename[x]!=0 && filename[x]!='\n' && x<249; x++);
+ for (x=0; txsite.filename[x]!='.' && txsite.filename[x]!=0 && x<250; x++)
+ filename[x]=txsite.filename[x];
filename[x]='.';
filename[x+1]='l';
ok=sscanf(string,"%lf", &din);
}
- fclose(fd);
-
if (ok)
{
LR.rel=din;
- LoadPAT(filename);
+ din=0.0;
+ return_value=1;
+
+ if (fgets(string,80,fd)!=NULL)
+ {
+ pointer=strchr(string,':');
+
+ if (pointer!=NULL)
+ *pointer=0;
+
+ if (sscanf(string,"%lf", &din))
+ LR.erp=din;
+ }
}
+
+ fclose(fd);
+
+ if (forced_erp!=-1.0)
+ LR.erp=forced_erp;
+
+ if (ok)
+ LoadPAT(filename);
}
- if (fd==NULL)
+ if (fd==NULL && forced_read)
{
- /* Create a "splat.lrp" file since one
- could not be successfully loaded. */
+ /* Assign some default parameters
+ for use in this run. */
+
+ LR.eps_dielect=15.0;
+ LR.sgm_conductivity=0.005;
+ LR.eno_ns_surfref=301.0;
+ LR.frq_mhz=300.0;
+ LR.radio_climate=5;
+ LR.pol=0;
+ LR.conf=0.50;
+ LR.rel=0.50;
+ LR.erp=0.0;
+
+ /* Write them to a "splat.lrp" file. */
outfile=fopen("splat.lrp","w");
fprintf(outfile,"%d\t; Radio Climate\n",LR.radio_climate);
fprintf(outfile,"%d\t; Polarization (0 = Horizontal, 1 = Vertical)\n", LR.pol);
fprintf(outfile,"%.2f\t; Fraction of situations\n",LR.conf);
- fprintf(outfile, "%.2f\t; Fraction of time\n",LR.rel);
+ fprintf(outfile,"%.2f\t; Fraction of time\n",LR.rel);
+ fprintf(outfile,"%.2f\t; Transmitter Effective Radiated Power in Watts (optional)\n",LR.erp);
fprintf(outfile,"\nPlease consult SPLAT! documentation for the meaning and use of this data.\n");
fclose(outfile);
- fprintf(stderr,"\n%c*** There were problems reading your \"%s\" file! ***\nA \"splat.lrp\" file was written to your directory with default data.\n",7,filename);
+ return_value=1;
+
+ fprintf(stderr,"\n\n%c*** There were problems reading your \"%s\" file! ***\nA \"splat.lrp\" file was written to your directory with default data.\n",7,filename);
}
- if (fd==NULL || ok==0)
+ else if (forced_read==0)
+ return_value=0;
+
+ if (forced_read && (fd==NULL || ok==0))
+ {
+ LR.eps_dielect=15.0;
+ LR.sgm_conductivity=0.005;
+ LR.eno_ns_surfref=301.0;
+ LR.frq_mhz=300.0;
+ LR.radio_climate=5;
+ LR.pol=0;
+ LR.conf=0.50;
+ LR.rel=0.50;
+ LR.erp=0.0;
+
fprintf(stderr,"Longley-Rice default parameters have been assumed for this analysis.\n");
+
+ return_value=1;
+ }
+
+ return (return_value);
}
void PlotPath(struct site source, struct site destination, char mask_value)
}
}
-void PlotLRPath(struct site source, struct site destination, FILE *fd)
+void PlotLRPath(struct site source, struct site destination, unsigned char mask_value, FILE *fd)
{
/* This function plots the RF path loss between source and
destination points based on the Longley-Rice propagation
model, taking into account antenna pattern data, if available. */
+ int x, y, ifs, ofs, errnum;
char block=0, strmode[100];
- int x, y, errnum;
double loss, azimuth, pattern=0.0,
- source_alt, dest_alt, source_alt2, dest_alt2,
- cos_xmtr_angle, cos_test_angle=0.0, test_alt,
- elevation, distance=0.0, four_thirds_earth;
+ xmtr_alt, dest_alt, xmtr_alt2, dest_alt2,
+ cos_rcvr_angle, cos_test_angle=0.0, test_alt,
+ elevation=0.0, distance=0.0, four_thirds_earth,
+ field_strength=0.0;
struct site temp;
ReadPath(source,destination);
/* Process this point only if it
has not already been processed. */
- if (GetMask(path.lat[y],path.lon[y])==0)
+ if ((GetMask(path.lat[y],path.lon[y])&248)!=(mask_value<<3))
{
distance=5280.0*path.distance[y];
- source_alt=four_thirds_earth+source.alt+path.elevation[0];
+ xmtr_alt=four_thirds_earth+source.alt+path.elevation[0];
dest_alt=four_thirds_earth+destination.alt+path.elevation[y];
dest_alt2=dest_alt*dest_alt;
- source_alt2=source_alt*source_alt;
+ xmtr_alt2=xmtr_alt*xmtr_alt;
/* Calculate the cosine of the elevation of
the receiver as seen by the transmitter. */
- cos_xmtr_angle=((source_alt2)+(distance*distance)-(dest_alt2))/(2.0*source_alt*distance);
+ cos_rcvr_angle=((xmtr_alt2)+(distance*distance)-(dest_alt2))/(2.0*xmtr_alt*distance);
+
+ if (cos_rcvr_angle>1.0)
+ cos_rcvr_angle=1.0;
+
+ if (cos_rcvr_angle<-1.0)
+ cos_rcvr_angle=-1.0;
if (got_elevation_pattern || fd!=NULL)
{
- /* If no antenna elevation pattern is available, and
- no output file is designated, the following code
- that determines the elevation angle to the first
- obstruction along the path is bypassed. */
+ /* Determine the elevation angle to the first obstruction
+ along the path IF elevation pattern data is available
+ or an output (.plo) file has been designated. */
for (x=2, block=0; (x<y && block==0); x++)
{
distance=5280.0*path.distance[x];
+
test_alt=four_thirds_earth+path.elevation[x];
/* Calculate the cosine of the elevation
angle of the terrain (test point)
as seen by the transmitter. */
- cos_test_angle=((source_alt2)+(distance*distance)-(test_alt*test_alt))/(2.0*source_alt*distance);
+ cos_test_angle=((xmtr_alt2)+(distance*distance)-(test_alt*test_alt))/(2.0*xmtr_alt*distance);
+
+ if (cos_test_angle>1.0)
+ cos_test_angle=1.0;
+
+ if (cos_test_angle<-1.0)
+ cos_test_angle=-1.0;
/* Compare these two angles to determine if
an obstruction exists. Since we're comparing
the cosines of these angles rather than
the angles themselves, the sense of the
following "if" statement is reversed from
- what it would be if the angles themselves
+ what it would be if the angles themselves
were compared. */
- if (cos_xmtr_angle>cos_test_angle)
+ if (cos_rcvr_angle>cos_test_angle)
block=1;
}
- /* At this point, we have the elevation angle
- to the first obstruction (if it exists). */
+ if (block)
+ elevation=((acos(cos_test_angle))/deg2rad)-90.0;
+ else
+ elevation=((acos(cos_rcvr_angle))/deg2rad)-90.0;
}
/* Determine attenuation for each point along the
LR.radio_climate, LR.pol, LR.conf, LR.rel, loss,
strmode, errnum);
- if (block)
- elevation=((acos(cos_test_angle))/deg2rad)-90.0;
-
- else
- elevation=((acos(cos_xmtr_angle))/deg2rad)-90.0;
-
temp.lat=path.lat[y];
temp.lon=path.lon[y];
azimuth=(Azimuth(source,temp));
- if (fd!=NULL)
- {
- /* Write path loss data to output file */
+ /* Write path loss data to output file */
- fprintf(fd,"%.7f, %.7f, %.3f, %.3f, %.2f\n",path.lat[y], path.lon[y], azimuth, elevation, loss);
- }
+ if (fd!=NULL)
+ fprintf(fd,"%.7f, %.7f, %.3f, %.3f, %.2f",path.lat[y], path.lon[y], azimuth, elevation, loss);
/* Integrate the antenna's radiation
pattern into the overall path loss. */
{
pattern=20.0*log10(pattern);
loss-=pattern;
+
+ if (fd!=NULL && (got_elevation_pattern || got_azimuth_pattern))
+ fprintf(fd,", %.2f",loss);
}
}
- if (loss>225.0)
- loss=225.0;
+ if (LR.erp!=0.0)
+ {
+ field_strength=(137.26+(20.0*log10(LR.frq_mhz))-loss)+(10.0*log10(LR.erp/1000.0));
+
+ ifs=100+(int)rint(field_strength);
- if (loss<75.0)
- loss=75.0;
+ if (ifs<0)
+ ifs=0;
- loss-=75.0;
- loss/=10.0;
- loss+=1.0;
+ if (ifs>255)
+ ifs=255;
+
+ ofs=GetSignal(path.lat[y],path.lon[y]);
+
+ if (ofs>ifs)
+ ifs=ofs;
+
+ PutSignal(path.lat[y],path.lon[y],(unsigned char)ifs);
- OrMask(path.lat[y],path.lon[y],((unsigned char)(loss))<<3);
- }
+ if (fd!=NULL)
+ fprintf(fd,", %.3f",field_strength);
+ }
+
+ else
+ {
+ if (loss>255)
+ ifs=255;
+ else
+ ifs=(int)rint(loss);
+
+ ofs=GetSignal(path.lat[y],path.lon[y]);
+
+ if (ofs<ifs && ofs!=0)
+ ifs=ofs;
+
+ PutSignal(path.lat[y],path.lon[y],(unsigned char)ifs);
+ }
+
+ if (fd!=NULL)
+ {
+ if (block)
+ fprintf(fd," *");
- else if (GetMask(path.lat[y],path.lon[y])==0 && path.distance[y]>max_range)
- OrMask(path.lat[y],path.lon[y],1);
+ fprintf(fd,"\n");
+ }
+
+ /* Mark this point as being analyzed */
+
+ PutMask(path.lat[y],path.lon[y],(GetMask(path.lat[y],path.lon[y])&7)+mask_value<<3);
+ }
}
}
is later invoked. */
float lat, lon, one_pixel;
- static unsigned char mask_value;
+ static unsigned char mask_value=1;
int z, count;
struct site edge;
unsigned char symbol[4], x;
- /* Initialize mask_value */
-
- if (mask_value!=8 && mask_value!=16 && mask_value!=32)
- mask_value=1;
-
one_pixel=1.0/1200.0;
symbol[0]='.';
count=0;
- fprintf(stdout,"\nComputing line-of-sight coverage of %s with an RX antenna\nat %.2f %s AGL:\n\n 0%c to 25%c ",source.name,metric?altitude*METERS_PER_FOOT:altitude,metric?"meters":"feet",37,37);
+ fprintf(stdout,"\n\nComputing line-of-sight coverage of \"%s\" with an RX antenna\nat %.2f %s AGL...\n\n 0%c to 25%c ",source.name,metric?altitude*METERS_PER_FOOT:altitude,metric?"meters":"feet",37,37);
fflush(stdout);
/* 18.75=1200 pixels/degree divided by 64 loops
topographic map based on a receiver located at
the specified altitude (in feet AGL). Results
are stored in memory, and written out in the form
- of a topographic map when the WritePPMLR() function
- is later invoked. */
+ of a topographic map when the WritePPMLR() or
+ WritePPMSS() functions are later invoked. */
int z, count;
struct site edge;
float lat, lon, one_pixel;
- unsigned char symbol[4], x;
+ unsigned char x, symbol[4];
+ static unsigned char mask_value=1;
FILE *fd=NULL;
one_pixel=1.0/1200.0;
count=0;
- fprintf(stdout,"\nComputing Longley-Rice coverage of %s ", source.name);
+ fprintf(stdout,"\n\nComputing Longley-Rice contours of \"%s\" ", source.name);
- fprintf(stdout,"out to a radius\nof %.2f %s with an RX antenna at %.2f %s AGL:\n\n 0%c to 25%c ",metric?max_range*KM_PER_MILE:max_range,metric?"kilometers":"miles",metric?altitude*METERS_PER_FOOT:altitude,metric?"meters":"feet",37,37);
+ fprintf(stdout,"out to a radius\nof %.2f %s with an RX antenna at %.2f %s AGL...\n\n 0%c to 25%c ",metric?max_range*KM_PER_MILE:max_range,metric?"kilometers":"miles",metric?altitude*METERS_PER_FOOT:altitude,metric?"meters":"feet",37,37);
fflush(stdout);
if (plo_filename[0]!=0)
edge.lon=lon;
edge.alt=altitude;
- PlotLRPath(source,edge,fd);
+ PlotLRPath(source,edge,mask_value,fd);
count++;
if (count==z)
edge.lon=min_west;
edge.alt=altitude;
- PlotLRPath(source,edge,fd);
+ PlotLRPath(source,edge,mask_value,fd);
count++;
if (count==z)
edge.lon=lon;
edge.alt=altitude;
- PlotLRPath(source,edge,fd);
+ PlotLRPath(source,edge,mask_value,fd);
count++;
if (count==z)
edge.lon=max_west;
edge.alt=altitude;
- PlotLRPath(source,edge,fd);
+ PlotLRPath(source,edge,mask_value,fd);
count++;
if (count==z)
fprintf(stdout,"\nDone!\n");
fflush(stdout);
+
+ if (mask_value<30)
+ mask_value++;
}
-void WritePPM(char *filename, unsigned char geo)
+void LoadSignalColors(struct site xmtr)
{
- /* This function generates a topographic map in Portable Pix Map
- (PPM) format based on logarithmically scaled topology data,
- as well as the content of flags held in the mask[][] array.
- The image created is rotated counter-clockwise 90 degrees
- from its representation in dem[][] so that north points
- up and east points right in the image generated. */
-
- char mapfile[255], geofile[255];
- unsigned char found, mask;
- unsigned width, height, output;
- int indx, x, y, x0=0, y0=0, minlat, minlon;
- float lat, lon, one_pixel, conversion, one_over_gamma;
- FILE *fd;
+ int x, y, ok, val[4];
+ char filename[255], string[80], *pointer=NULL;
+ FILE *fd=NULL;
- one_pixel=1.0/1200.0;
- one_over_gamma=1.0/GAMMA;
- conversion=255.0/pow((double)(max_elevation-min_elevation),one_over_gamma);
+ for (x=0; xmtr.filename[x]!='.' && xmtr.filename[x]!=0 && x<250; x++)
+ filename[x]=xmtr.filename[x];
- width=(unsigned)(1200*ReduceAngle(max_west-min_west));
- height=(unsigned)(1200*ReduceAngle(max_north-min_north));
+ filename[x]='.';
+ filename[x+1]='s';
+ filename[x+2]='c';
+ filename[x+3]='f';
+ filename[x+4]=0;
- if (filename[0]==0)
- strncpy(mapfile, "map.ppm\0",8);
- else
- {
- for (x=0; filename[x]!='.' && filename[x]!=0 && x<250; x++)
- {
- mapfile[x]=filename[x];
- geofile[x]=filename[x];
- }
+ /* Default values */
+
+ region.level[0]=128;
+ region.color[0][0]=255;
+ region.color[0][1]=0;
+ region.color[0][2]=0;
+
+ region.level[1]=118;
+ region.color[1][0]=255;
+ region.color[1][1]=165;
+ region.color[1][2]=0;
+
+ region.level[2]=108;
+ region.color[2][0]=255;
+ region.color[2][1]=206;
+ region.color[2][2]=0;
+
+ region.level[3]=98;
+ region.color[3][0]=255;
+ region.color[3][1]=255;
+ region.color[3][2]=0;
+
+ region.level[4]=88;
+ region.color[4][0]=184;
+ region.color[4][1]=255;
+ region.color[4][2]=0;
+
+ region.level[5]=78;
+ region.color[5][0]=0;
+ region.color[5][1]=255;
+ region.color[5][2]=0;
+
+ region.level[6]=68;
+ region.color[6][0]=0;
+ region.color[6][1]=208;
+ region.color[6][2]=0;
+
+ region.level[7]=58;
+ region.color[7][0]=0;
+ region.color[7][1]=196;
+ region.color[7][2]=196;
+
+ region.level[8]=48;
+ region.color[8][0]=0;
+ region.color[8][1]=148;
+ region.color[8][2]=255;
+
+ region.level[9]=38;
+ region.color[9][0]=80;
+ region.color[9][1]=80;
+ region.color[9][2]=255;
+
+ region.level[10]=28;
+ region.color[10][0]=0;
+ region.color[10][1]=38;
+ region.color[10][2]=255;
+
+ region.level[11]=18;
+ region.color[11][0]=142;
+ region.color[11][1]=63;
+ region.color[11][2]=255;
+
+ region.level[12]=8;
+ region.color[12][0]=140;
+ region.color[12][1]=0;
+ region.color[12][2]=128;
+
+ region.levels=13;
+
+ fd=fopen("splat.scf","r");
- mapfile[x]='.';
- geofile[x]='.';
- mapfile[x+1]='p';
- geofile[x+1]='g';
- mapfile[x+2]='p';
- geofile[x+2]='e';
- mapfile[x+3]='m';
- geofile[x+3]='o';
- mapfile[x+4]=0;
- geofile[x+4]=0;
- }
+ if (fd==NULL)
+ fd=fopen(filename,"r");
- if (geo)
+ if (fd==NULL)
{
- fd=fopen(geofile,"wb");
+ fd=fopen(filename,"w");
- fprintf(fd,"FILENAME\t%s\n",mapfile);
- fprintf(fd,"#\t\tX\tY\tLong\t\tLat\n");
- fprintf(fd,"TIEPOINT\t0\t0\t%d.000\t\t%d.000\n",(max_west<180?-max_west:360-max_west),max_north);
- fprintf(fd,"TIEPOINT\t%u\t%u\t%d.000\t\t%d.000\n",width-1,height-1,(min_west<180?-min_west:360-min_west),min_north);
- fprintf(fd,"IMAGESIZE\t%u\t%u\n",width,height);
- fprintf(fd,"#\n# Auto Generated by SPLAT! v%s\n#\n",splat_version);
+ fprintf(fd,"; SPLAT! Auto-generated Signal Color Definition (\"%s\") File\n",filename);
+ fprintf(fd,";\n; Format for the parameters held in this file is as follows:\n;\n");
+ fprintf(fd,"; dBuV/m: red, green, blue\n;\n");
+ fprintf(fd,"; ...where \"dBuV/m\" is the signal strength (in dBuV/m) and\n");
+ fprintf(fd,"; \"red\", \"green\", and \"blue\" are the corresponding RGB color\n");
+ fprintf(fd,"; definitions ranging from 0 to 255 for the region specified.\n");
+ fprintf(fd,";\n; The following parameters may be edited and/or expanded\n");
+ fprintf(fd,"; for future runs of SPLAT! A total of 32 contour regions\n");
+ fprintf(fd,"; may be defined in this file.\n;\n;\n");
+
+ for (x=0; x<region.levels; x++)
+ fprintf(fd,"%3d: %3d, %3d, %3d\n",region.level[x], region.color[x][0], region.color[x][1], region.color[x][2]);
fclose(fd);
}
- fd=fopen(mapfile,"wb");
-
- fprintf(fd,"P6\n%u %u\n255\n",width,height);
-
- fprintf(stdout,"\nWriting \"%s\" (%ux%u pixmap image)... ",mapfile,width,height);
- fflush(stdout);
-
- for (y=0, lat=((double)max_north)-one_pixel; y<(int)height; y++, lat-=one_pixel)
+ else
{
- minlat=(int)floor(lat);
+ x=0;
+ fgets(string,80,fd);
- for (x=0, lon=((double)max_west)-one_pixel; x<(int)width; x++, lon-=one_pixel)
+ while (x<32 && feof(fd)==0)
{
- if (lon<0.0)
- lon+=360.0;
+ pointer=strchr(string,';');
- minlon=(int)floor(lon);
+ if (pointer!=NULL)
+ *pointer=0;
- for (indx=0, found=0; indx<MAXSLOTS && found==0;)
- if (minlat==dem[indx].min_north && minlon==dem[indx].min_west)
- found=1;
- else
- indx++;
+ ok=sscanf(string,"%d: %d, %d, %d", &val[0], &val[1], &val[2], &val[3]);
- if (found)
+ if (ok==4)
{
- x0=(int)(1199.0*(lat-floor(lat)));
- y0=(int)(1199.0*(lon-floor(lon)));
+ for (y=0; y<4; y++)
+ {
+ if (val[y]>255)
+ val[y]=255;
- mask=dem[indx].mask[x0][y0];
+ if (val[y]<0)
+ val[y]=0;
+ }
+
+ region.level[x]=val[0];
+ region.color[x][0]=val[1];
+ region.color[x][1]=val[2];
+ region.color[x][2]=val[3];
+ x++;
+ }
- if (mask&2)
- /* Text Labels: Red */
- fprintf(fd,"%c%c%c",255,0,0);
+ fgets(string,80,fd);
+ }
- else if (mask&4)
- /* County Boundaries: Light Cyan */
- fprintf(fd,"%c%c%c",128,128,255);
+ fclose(fd);
+ region.levels=x;
+ }
+}
- else switch (mask&57)
- {
- case 1:
- /* TX1: Green */
- fprintf(fd,"%c%c%c",0,255,0);
- break;
+void LoadLossColors(struct site xmtr)
+{
+ int x, y, ok, val[4];
+ char filename[255], string[80], *pointer=NULL;
+ FILE *fd=NULL;
- case 8:
- /* TX2: Cyan */
- fprintf(fd,"%c%c%c",0,255,255);
- break;
+ for (x=0; xmtr.filename[x]!='.' && xmtr.filename[x]!=0 && x<250; x++)
+ filename[x]=xmtr.filename[x];
- case 9:
- /* TX1 + TX2: Yellow */
- fprintf(fd,"%c%c%c",255,255,0);
- break;
+ filename[x]='.';
+ filename[x+1]='l';
+ filename[x+2]='c';
+ filename[x+3]='f';
+ filename[x+4]=0;
- case 16:
- /* TX3: Medium Violet */
- fprintf(fd,"%c%c%c",147,112,219);
- break;
+ /* Default values */
+
+ region.level[0]=80;
+ region.color[0][0]=255;
+ region.color[0][1]=0;
+ region.color[0][2]=0;
+
+ region.level[1]=90;
+ region.color[1][0]=255;
+ region.color[1][1]=128;
+ region.color[1][2]=0;
+
+ region.level[2]=100;
+ region.color[2][0]=255;
+ region.color[2][1]=165;
+ region.color[2][2]=0;
+
+ region.level[3]=110;
+ region.color[3][0]=255;
+ region.color[3][1]=206;
+ region.color[3][2]=0;
+
+ region.level[4]=120;
+ region.color[4][0]=255;
+ region.color[4][1]=255;
+ region.color[4][2]=0;
+
+ region.level[5]=130;
+ region.color[5][0]=184;
+ region.color[5][1]=255;
+ region.color[5][2]=0;
+
+ region.level[6]=140;
+ region.color[6][0]=0;
+ region.color[6][1]=255;
+ region.color[6][2]=0;
+
+ region.level[7]=150;
+ region.color[7][0]=0;
+ region.color[7][1]=208;
+ region.color[7][2]=0;
+
+ region.level[8]=160;
+ region.color[8][0]=0;
+ region.color[8][1]=196;
+ region.color[8][2]=196;
+
+ region.level[9]=170;
+ region.color[9][0]=0;
+ region.color[9][1]=148;
+ region.color[9][2]=255;
+
+ region.level[10]=180;
+ region.color[10][0]=80;
+ region.color[10][1]=80;
+ region.color[10][2]=255;
+
+ region.level[11]=190;
+ region.color[11][0]=0;
+ region.color[11][1]=38;
+ region.color[11][2]=255;
+
+ region.level[12]=200;
+ region.color[12][0]=142;
+ region.color[12][1]=63;
+ region.color[12][2]=255;
+
+ region.level[13]=210;
+ region.color[13][0]=196;
+ region.color[13][1]=54;
+ region.color[13][2]=255;
+
+ region.level[14]=220;
+ region.color[14][0]=255;
+ region.color[14][1]=0;
+ region.color[14][2]=255;
+
+ region.level[15]=230;
+ region.color[15][0]=255;
+ region.color[15][1]=194;
+ region.color[15][2]=204;
+
+ region.levels=16;
+
+ fd=fopen("splat.lcf","r");
+
+ if (fd==NULL)
+ fd=fopen(filename,"r");
+
+ if (fd==NULL)
+ {
+ fd=fopen(filename,"w");
+
+ fprintf(fd,"; SPLAT! Auto-generated Path-Loss Color Definition (\"%s\") File\n",filename);
+ fprintf(fd,";\n; Format for the parameters held in this file is as follows:\n;\n");
+ fprintf(fd,"; dB: red, green, blue\n;\n");
+ fprintf(fd,"; ...where \"dB\" is the path loss (in dB) and\n");
+ fprintf(fd,"; \"red\", \"green\", and \"blue\" are the corresponding RGB color\n");
+ fprintf(fd,"; definitions ranging from 0 to 255 for the region specified.\n");
+ fprintf(fd,";\n; The following parameters may be edited and/or expanded\n");
+ fprintf(fd,"; for future runs of SPLAT! A total of 32 contour regions\n");
+ fprintf(fd,"; may be defined in this file.\n;\n;\n");
+
+ for (x=0; x<region.levels; x++)
+ fprintf(fd,"%3d: %3d, %3d, %3d\n",region.level[x], region.color[x][0], region.color[x][1], region.color[x][2]);
+
+ fclose(fd);
+ }
+
+ else
+ {
+ x=0;
+ fgets(string,80,fd);
+
+ while (x<32 && feof(fd)==0)
+ {
+ pointer=strchr(string,';');
+
+ if (pointer!=NULL)
+ *pointer=0;
+
+ ok=sscanf(string,"%d: %d, %d, %d", &val[0], &val[1], &val[2], &val[3]);
+
+ if (ok==4)
+ {
+ for (y=0; y<4; y++)
+ {
+ if (val[y]>255)
+ val[y]=255;
+
+ if (val[y]<0)
+ val[y]=0;
+ }
+
+ region.level[x]=val[0];
+ region.color[x][0]=val[1];
+ region.color[x][1]=val[2];
+ region.color[x][2]=val[3];
+ x++;
+ }
+
+ fgets(string,80,fd);
+ }
+
+ fclose(fd);
+ region.levels=x;
+ }
+}
+
+void WritePPM(char *filename, unsigned char geo, unsigned char kml, unsigned char ngs)
+{
+ /* This function generates a topographic map in Portable Pix Map
+ (PPM) format based on logarithmically scaled topology data,
+ as well as the content of flags held in the mask[][] array.
+ The image created is rotated counter-clockwise 90 degrees
+ from its representation in dem[][] so that north points
+ up and east points right in the image generated. */
+
+ char mapfile[255], geofile[255], kmlfile[255];
+ unsigned char found, mask;
+ unsigned width, height, terrain;
+ int indx, x, y, x0=0, y0=0;
+ double lat, lon, one_pixel, conversion, one_over_gamma; /* USED to be float... */
+ FILE *fd;
+
+ one_pixel=1.0/1200.0;
+ one_over_gamma=1.0/GAMMA;
+ conversion=255.0/pow((double)(max_elevation-min_elevation),one_over_gamma);
+
+ width=(unsigned)(1200*ReduceAngle(max_west-min_west));
+ height=(unsigned)(1200*ReduceAngle(max_north-min_north));
+
+ if (filename[0]==0)
+ strncpy(filename, "map.ppm\0",8);
+
+ for (x=0; filename[x]!='.' && filename[x]!=0 && x<250; x++)
+ {
+ mapfile[x]=filename[x];
+ geofile[x]=filename[x];
+ kmlfile[x]=filename[x];
+ }
+
+ mapfile[x]='.';
+ geofile[x]='.';
+ kmlfile[x]='.';
+ mapfile[x+1]='p';
+ geofile[x+1]='g';
+ kmlfile[x+1]='k';
+ mapfile[x+2]='p';
+ geofile[x+2]='e';
+ kmlfile[x+2]='m';
+ mapfile[x+3]='m';
+ geofile[x+3]='o';
+ kmlfile[x+3]='l';
+ mapfile[x+4]=0;
+ geofile[x+4]=0;
+ kmlfile[x+4]=0;
+
+ if (kml==0 && geo)
+ {
+ fd=fopen(geofile,"wb");
+
+ fprintf(fd,"FILENAME\t%s\n",mapfile);
+ fprintf(fd,"#\t\tX\tY\tLong\t\tLat\n");
+ fprintf(fd,"TIEPOINT\t0\t0\t%d.000\t\t%d.000\n",(max_west<180?-max_west:360-max_west),max_north);
+ fprintf(fd,"TIEPOINT\t%u\t%u\t%d.000\t\t%d.000\n",width-1,height-1,(min_west<180?-min_west:360-min_west),min_north);
+ fprintf(fd,"IMAGESIZE\t%u\t%u\n",width,height);
+ fprintf(fd,"#\n# Auto Generated by SPLAT! v%s\n#\n",splat_version);
+
+ fclose(fd);
+ }
+
+ if (kml && geo==0)
+ {
+ fd=fopen(kmlfile,"wb");
+
+ fprintf(fd,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
+ fprintf(fd,"<kml xmlns=\"http://earth.google.com/kml/2.1\">\n");
+ fprintf(fd," <Folder>\n");
+ fprintf(fd," <name>SPLAT!</name>\n");
+ fprintf(fd," <description>Line-of-Sight Overlay</description>\n");
+ fprintf(fd," <GroundOverlay>\n");
+ fprintf(fd," <name>SPLAT! Line-of-Sight Overlay</name>\n");
+ fprintf(fd," <description>SPLAT! Coverage</description>\n");
+ fprintf(fd," <Icon>\n");
+ fprintf(fd," <href>%s</href>\n",mapfile);
+ fprintf(fd," </Icon>\n");
+ fprintf(fd," <opacity>128</opacity>\n");
+ fprintf(fd," <LatLonBox>\n");
+ fprintf(fd," <north>%.5f</north>\n",(double)max_north-one_pixel);
+ fprintf(fd," <south>%.5f</south>\n",(double)min_north);
+ fprintf(fd," <east>%.5f</east>\n",((double)min_west<180.0?(double)-min_west:360.0-(double)min_west));
+ fprintf(fd," <west>%.5f</west>\n",(((double)max_west-one_pixel)<180.0?-((double)max_west-one_pixel):(360.0-(double)max_west-one_pixel)));
+ fprintf(fd," <rotation>0.0</rotation>\n");
+ fprintf(fd," </LatLonBox>\n");
+ fprintf(fd," </GroundOverlay>\n");
+ fprintf(fd," </Folder>\n");
+ fprintf(fd,"</kml>\n");
+
+ fclose(fd);
+ }
+
+ fd=fopen(mapfile,"wb");
+
+ fprintf(fd,"P6\n%u %u\n255\n",width,height);
+ fprintf(stdout,"\nWriting \"%s\" (%ux%u pixmap image)... ",mapfile,width,height);
+ fflush(stdout);
+
+ for (y=0, lat=(double)max_north-one_pixel; y<(int)height; y++, lat=(double)max_north-(one_pixel*(double)y))
+ {
+ for (x=0, lon=(double)max_west-one_pixel; x<(int)width; x++, lon=(double)max_west-(one_pixel*(double)x))
+ {
+ if (lon<0.0)
+ lon+=360.0;
+
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<(double)dem[indx].max_north && LonDiff(lon,(double)dem[indx].min_west)>=0.0 && LonDiff(lon,(double)dem[indx].max_west)<0.0)
+ found=1;
+ else
+ indx++;
+
+ if (found)
+ {
+ x0=(int)(1199.0*(lat-floor(lat)));
+ y0=(int)(1199.0*(lon-floor(lon)));
+
+ mask=dem[indx].mask[x0][y0];
+
+ if (mask&2)
+ /* Text Labels: Red */
+ fprintf(fd,"%c%c%c",255,0,0);
+
+ else if (mask&4)
+ /* County Boundaries: Light Cyan */
+ fprintf(fd,"%c%c%c",128,128,255);
+
+ else switch (mask&57)
+ {
+ case 1:
+ /* TX1: Green */
+ fprintf(fd,"%c%c%c",0,255,0);
+ break;
+
+ case 8:
+ /* TX2: Cyan */
+ fprintf(fd,"%c%c%c",0,255,255);
+ break;
+
+ case 9:
+ /* TX1 + TX2: Yellow */
+ fprintf(fd,"%c%c%c",255,255,0);
+ break;
+
+ case 16:
+ /* TX3: Medium Violet */
+ fprintf(fd,"%c%c%c",147,112,219);
+ break;
case 17:
/* TX1 + TX3: Pink */
break;
default:
- /* Water: Medium Blue */
- if (dem[indx].data[x0][y0]==0)
- fprintf(fd,"%c%c%c",0,0,170);
+ if (ngs) /* No terrain */
+ fprintf(fd,"%c%c%c",255,255,255);
else
{
- /* Elevation: Greyscale */
- output=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
- fprintf(fd,"%c%c%c",output,output,output);
+ /* Sea-level: Medium Blue */
+ if (dem[indx].data[x0][y0]==0)
+ fprintf(fd,"%c%c%c",0,0,170);
+ else
+ {
+ /* Elevation: Greyscale */
+ terrain=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
+ fprintf(fd,"%c%c%c",terrain,terrain,terrain);
+ }
}
}
}
fflush(stdout);
}
-void WritePPMLR(char *filename, unsigned char geo)
+void WritePPMLR(char *filename, unsigned char geo, unsigned char kml, unsigned char ngs, struct site *xmtr, unsigned char txsites)
{
/* This function generates a topographic map in Portable Pix Map
(PPM) format based on the content of flags held in the mask[][]
90 degrees from its representation in dem[][] so that north
points up and east points right in the image generated. */
- char mapfile[255], geofile[255];
- unsigned width, height, output;
+ char mapfile[255], geofile[255], kmlfile[255], color=0;
+ unsigned width, height, red, green, blue, terrain=0;
unsigned char found, mask, cityorcounty;
- int indx, x, y, t, t2, x0, y0, minlat, minlon, loss;
- float lat, lon, one_pixel, conversion, one_over_gamma;
+ int indx, x, y, z, colorwidth, x0, y0, loss, level,
+ hundreds, tens, units, match;
+ double lat, lon, one_pixel, conversion, one_over_gamma;
FILE *fd;
one_pixel=1.0/1200.0;
width=(unsigned)(1200*ReduceAngle(max_west-min_west));
height=(unsigned)(1200*ReduceAngle(max_north-min_north));
- if (filename[0]==0)
- strncpy(mapfile, "map.ppm\0",8);
- else
- {
- for (x=0; filename[x]!='.' && filename[x]!=0 && x<250; x++)
- {
- mapfile[x]=filename[x];
- geofile[x]=filename[x];
- }
+ LoadLossColors(xmtr[0]);
- mapfile[x]='.';
- geofile[x]='.';
- mapfile[x+1]='p';
- geofile[x+1]='g';
- mapfile[x+2]='p';
- geofile[x+2]='e';
- mapfile[x+3]='m';
- geofile[x+3]='o';
- mapfile[x+4]=0;
- geofile[x+4]=0;
- }
+ if (filename[0]==0)
+ strncpy(filename, xmtr[0].filename,254);
- if (geo)
+ for (x=0; filename[x]!='.' && filename[x]!=0 && x<250; x++)
+ {
+ mapfile[x]=filename[x];
+ geofile[x]=filename[x];
+ kmlfile[x]=filename[x];
+ }
+
+ mapfile[x]='.';
+ geofile[x]='.';
+ kmlfile[x]='.';
+ mapfile[x+1]='p';
+ geofile[x+1]='g';
+ kmlfile[x+1]='k';
+ mapfile[x+2]='p';
+ geofile[x+2]='e';
+ kmlfile[x+2]='m';
+ mapfile[x+3]='m';
+ geofile[x+3]='o';
+ kmlfile[x+3]='l';
+ mapfile[x+4]=0;
+ geofile[x+4]=0;
+ kmlfile[x+4]=0;
+
+ if (kml==0 && geo)
{
fd=fopen(geofile,"wb");
fclose(fd);
}
- fd=fopen(mapfile,"wb");
+ if (kml && geo==0)
+ {
+ fd=fopen(kmlfile,"wb");
+
+ fprintf(fd,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
+ fprintf(fd,"<kml xmlns=\"http://earth.google.com/kml/2.1\">\n");
+ fprintf(fd,"<!-- Generated by SPLAT! Version %s -->\n",splat_version);
+ fprintf(fd," <Folder>\n");
+ fprintf(fd," <name>SPLAT!</name>\n");
+ fprintf(fd," <description>%s Transmitter Path Loss Overlay</description>\n",xmtr[0].name);
+ fprintf(fd," <GroundOverlay>\n");
+ fprintf(fd," <name>SPLAT! Path Loss Overlay</name>\n");
+ fprintf(fd," <description>SPLAT! Coverage</description>\n");
+ fprintf(fd," <Icon>\n");
+ fprintf(fd," <href>%s</href>\n",mapfile);
+ fprintf(fd," </Icon>\n");
+ fprintf(fd," <opacity>128</opacity>\n");
+ fprintf(fd," <LatLonBox>\n");
+ fprintf(fd," <north>%.5f</north>\n",(double)max_north-one_pixel);
+ fprintf(fd," <south>%.5f</south>\n",(double)min_north);
+ fprintf(fd," <east>%.5f</east>\n",((double)min_west<180.0?(double)-min_west:360.0-(double)min_west));
+ fprintf(fd," <west>%.5f</west>\n",(((double)max_west-one_pixel)<180.0?-((double)max_west-one_pixel):(360.0-(double)max_west-one_pixel)));
+ fprintf(fd," <rotation>0.0</rotation>\n");
+ fprintf(fd," </LatLonBox>\n");
+ fprintf(fd," </GroundOverlay>\n");
+
+ for (x=0; x<txsites; x++)
+ {
+ fprintf(fd," <Placemark>\n");
+ fprintf(fd," <name>%s</name>\n",xmtr[x].name);
+ fprintf(fd," <visibility>1</visibility>\n");
+ fprintf(fd," <Style>\n");
+ fprintf(fd," <IconStyle>\n");
+ fprintf(fd," <Icon>\n");
+ fprintf(fd," <href>root://icons/palette-5.png</href>\n");
+ fprintf(fd," <x>224</x>\n");
+ fprintf(fd," <y>224</y>\n");
+ fprintf(fd," <w>32</w>\n");
+ fprintf(fd," <h>32</h>\n");
+ fprintf(fd," </Icon>\n");
+ fprintf(fd," </IconStyle>\n");
+ fprintf(fd," </Style>\n");
+ fprintf(fd," <Point>\n");
+ fprintf(fd," <extrude>1</extrude>\n");
+ fprintf(fd," <altitudeMode>relativeToGround</altitudeMode>\n");
+ fprintf(fd," <coordinates>%f,%f,%f</coordinates>\n",(xmtr[x].lon<180.0?-xmtr[x].lon:360.0-xmtr[x].lon), xmtr[x].lat, xmtr[x].alt);
+ fprintf(fd," </Point>\n");
+ fprintf(fd," </Placemark>\n");
+ }
+
+ fprintf(fd," </Folder>\n");
+ fprintf(fd,"</kml>\n");
+
+ fclose(fd);
+ }
- fprintf(fd,"P6\n%u %u\n255\n",width,height+30);
+ fd=fopen(mapfile,"wb");
- fprintf(stdout,"\nWriting \"%s\" (%ux%u pixmap image)... ",mapfile,width,height+30);
+ fprintf(fd,"P6\n%u %u\n255\n",width,(kml?height:height+30));
+ fprintf(stdout,"\nWriting \"%s\" (%ux%u pixmap image)... ",mapfile,width,(kml?height:height+30));
fflush(stdout);
- for (y=0, lat=((double)max_north)-one_pixel; y<(int)height; y++, lat-=one_pixel)
+ for (y=0, lat=(double)max_north-one_pixel; y<(int)height; y++, lat=(double)max_north-(one_pixel*(double)y))
{
- minlat=(int)floor(lat);
-
- for (x=0, lon=((double)max_west)-one_pixel; x<(int)width; x++, lon-=one_pixel)
+ for (x=0, lon=(double)max_west-one_pixel; x<(int)width; x++, lon=(double)max_west-(one_pixel*(double)x))
{
if (lon<0.0)
lon+=360.0;
- minlon=(int)floor(lon);
-
- for (indx=0, found=0; indx<MAXSLOTS && found==0;)
- if (minlat==dem[indx].min_north && minlon==dem[indx].min_west)
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<(double)dem[indx].max_north && LonDiff(lon,(double)dem[indx].min_west)>=0.0 && LonDiff(lon,(double)dem[indx].max_west)<0.0)
found=1;
else
indx++;
+
if (found)
{
x0=(int)(1199.0*(lat-floor(lat)));
y0=(int)(1199.0*(lon-floor(lon)));
mask=dem[indx].mask[x0][y0];
- loss=70+(10*(int)((mask&248)>>3));
+ loss=(dem[indx].signal[x0][y0]);
cityorcounty=0;
- if (mask&2)
+ match=255;
+
+ red=0;
+ green=0;
+ blue=0;
+
+ if (loss<=region.level[0])
+ match=0;
+ else
+ {
+ for (z=1; (z<region.levels && match==255); z++)
+ {
+ if (loss>=region.level[z-1] && loss<region.level[z])
+ match=z;
+ }
+ }
+
+ if (match<region.levels)
{
- /* Text Labels - Black or Red */
+ red=region.color[match][0];
+ green=region.color[match][1];
+ blue=region.color[match][2];
- if ((mask&120) && (loss<=90))
- fprintf(fd,"%c%c%c",0,0,0);
- else
- fprintf(fd,"%c%c%c",255,0,0);
+ color=1;
+ }
- cityorcounty=1;
+ if ((mask&2) && (kml==0))
+ {
+ /* Text Labels: Red or otherwise */
+
+ if (red>=180 && green<=75 && blue<=75 && loss!=0)
+ fprintf(fd,"%c%c%c",255^red,255^green,255^blue);
+ else
+ fprintf(fd,"%c%c%c",255,0,0);
+
+ cityorcounty=1;
}
- else if (mask&4)
+ else if ((mask&4) && (kml==0))
{
/* County Boundaries: Black */
if (cityorcounty==0)
{
- if (loss>maxdB)
-
- { /* Display land or sea elevation */
-
- if (dem[indx].data[x0][y0]==0)
- fprintf(fd,"%c%c%c",0,0,170);
+ if (loss>maxdB || loss==0)
+ {
+ if (ngs) /* No terrain */
+ fprintf(fd,"%c%c%c",255,255,255);
else
{
- /* Elevation: Greyscale */
- output=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
- fprintf(fd,"%c%c%c",output,output,output);
+ /* Display land or sea elevation */
+
+ if (dem[indx].data[x0][y0]==0)
+ fprintf(fd,"%c%c%c",0,0,170);
+ else
+ {
+ terrain=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
+ fprintf(fd,"%c%c%c",terrain,terrain,terrain);
+ }
}
}
- else switch (loss)
+ else
{
- /* Plot signal loss in color */
-
- case 80:
- fprintf(fd,"%c%c%c",255,0,0);
- break;
-
- case 90:
- fprintf(fd,"%c%c%c",255,128,0);
- break;
-
- case 100:
- fprintf(fd,"%c%c%c",255,165,0);
- break;
-
- case 110:
- fprintf(fd,"%c%c%c",255,206,0);
- break;
-
- case 120:
- fprintf(fd,"%c%c%c",255,255,0);
- break;
-
- case 130:
- fprintf(fd,"%c%c%c",184,255,0);
- break;
-
- case 140:
- fprintf(fd,"%c%c%c",0,255,0);
- break;
-
- case 150:
- fprintf(fd,"%c%c%c",0,208,0);
- break;
-
- case 160:
- fprintf(fd,"%c%c%c",0,196,196);
- break;
-
- case 170:
- fprintf(fd,"%c%c%c",0,148,255);
- break;
-
- case 180:
- fprintf(fd,"%c%c%c",80,80,255);
- break;
-
- case 190:
- fprintf(fd,"%c%c%c",0,38,255);
- break;
-
- case 200:
- fprintf(fd,"%c%c%c",142,63,255);
- break;
-
- case 210:
- fprintf(fd,"%c%c%c",196,54,255);
- break;
-
- case 220:
- fprintf(fd,"%c%c%c",255,0,255);
- break;
-
- case 230:
- fprintf(fd,"%c%c%c",255,194,204);
- break;
+ /* Plot path loss in color */
- default:
+ if (red!=0 || green!=0 || blue!=0)
+ fprintf(fd,"%c%c%c",red,green,blue);
- if (dem[indx].data[x0][y0]==0)
- fprintf(fd,"%c%c%c",0,0,170);
- else
+ else /* terrain / sea-level */
{
- /* Elevation: Greyscale */
- output=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
- fprintf(fd,"%c%c%c",output,output,output);
+ if (dem[indx].data[x0][y0]==0)
+ fprintf(fd,"%c%c%c",0,0,170);
+ else
+ {
+ /* Elevation: Greyscale */
+ terrain=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
+ fprintf(fd,"%c%c%c",terrain,terrain,terrain);
+ }
}
}
}
}
}
- /* Display legend along bottom of image */
+ if (kml==0 && color)
+ {
+ /* Display legend along bottom of image */
- x0=width/16;
+ colorwidth=(int)rint((float)width/(float)region.levels);
- for (y0=0; y0<30; y0++)
- {
- for (indx=0; indx<16; indx++)
+ for (y0=0; y0<30; y0++)
{
- for (x=0; x<x0; x++)
+ for (x0=0; x0<(int)width; x0++)
{
- t=indx;
- t2=indx+8;
+ indx=x0/colorwidth;
+ x=x0%colorwidth;
+ level=region.level[indx];
+
+ hundreds=level/100;
+
+ if (hundreds>0)
+ level-=(hundreds*100);
+
+ tens=level/10;
+
+ if (tens>0)
+ level-=(tens*10);
+
+ units=level;
- if (y0>=10 && y0<=18)
+ if (y0>=8 && y0<=23)
{
- if (t2>9)
+ if (hundreds>0)
{
- if (x>=11 && x<=17)
- if (smallfont[t2/10][y0-10][x-11])
- t=255;
+ if (x>=11 && x<=18)
+ if (fontdata[16*(hundreds+'0')+(y0-8)]&(128>>(x-11)))
+ indx=255;
}
- if (x>=19 && x<=25)
- if (smallfont[t2%10][y0-10][x-19])
- t=255;
+ if (tens>0 || hundreds>0)
+ {
+ if (x>=19 && x<=26)
+ if (fontdata[16*(tens+'0')+(y0-8)]&(128>>(x-19)))
+ indx=255;
+ }
- if (x>=27 && x<=33)
- if (smallfont[0][y0-10][x-27])
- t=255;
+ if (x>=27 && x<=34)
+ if (fontdata[16*(units+'0')+(y0-8)]&(128>>(x-27)))
+ indx=255;
+
+ if (x>=42 && x<=49)
+ if (fontdata[16*('d')+(y0-8)]&(128>>(x-42)))
+ indx=255;
+
+ if (x>=50 && x<=57)
+ if (fontdata[16*('B')+(y0-8)]&(128>>(x-50)))
+ indx=255;
}
- switch (t)
+ if (indx>region.levels)
+ fprintf(fd,"%c%c%c",0,0,0);
+ else
{
- case 0:
- fprintf(fd,"%c%c%c",255,0,0);
- break;
+ red=region.color[indx][0];
+ green=region.color[indx][1];
+ blue=region.color[indx][2];
- case 1:
- fprintf(fd,"%c%c%c",255,128,0);
- break;
+ fprintf(fd,"%c%c%c",red,green,blue);
+ }
+ }
+ }
+ }
- case 2:
- fprintf(fd,"%c%c%c",255,165,0);
- break;
+ fclose(fd);
+ fprintf(stdout,"Done!\n");
+ fflush(stdout);
+}
- case 3:
- fprintf(fd,"%c%c%c",255,206,0);
- break;
+void WritePPMSS(char *filename, unsigned char geo, unsigned char kml, unsigned char ngs, struct site *xmtr, unsigned char txsites)
+{
+ /* This function generates a topographic map in Portable Pix Map
+ (PPM) format based on the signal strength values held in the
+ signal[][] array. The image created is rotated counter-clockwise
+ 90 degrees from its representation in dem[][] so that north
+ points up and east points right in the image generated. */
- case 4:
- fprintf(fd,"%c%c%c",255,255,0);
- break;
+ char mapfile[255], geofile[255], kmlfile[255], color=0;
+ unsigned width, height, terrain, red, green, blue;
+ unsigned char found, mask, cityorcounty;
+ int indx, x, y, z=1, x0, y0, signal, level, hundreds,
+ tens, units, match, colorwidth;
+ double lat, lon, one_pixel, conversion, one_over_gamma;
+ FILE *fd;
- case 5:
- fprintf(fd,"%c%c%c",184,255,0);
- break;
+ one_pixel=1.0/1200.0;
+ one_over_gamma=1.0/GAMMA;
+ conversion=255.0/pow((double)(max_elevation-min_elevation),one_over_gamma);
- case 6:
- fprintf(fd,"%c%c%c",0,255,0);
- break;
+ width=(unsigned)(1200*ReduceAngle(max_west-min_west));
+ height=(unsigned)(1200*ReduceAngle(max_north-min_north));
- case 7:
- fprintf(fd,"%c%c%c",0,208,0);
- break;
+ LoadSignalColors(xmtr[0]);
- case 8:
- fprintf(fd,"%c%c%c",0,196,196);
- break;
+ if (filename[0]==0)
+ strncpy(filename, xmtr[0].filename,254);
- case 9:
- fprintf(fd,"%c%c%c",0,148,255);
- break;
+ for (x=0; filename[x]!='.' && filename[x]!=0 && x<250; x++)
+ {
+ mapfile[x]=filename[x];
+ geofile[x]=filename[x];
+ kmlfile[x]=filename[x];
+ }
+
+ mapfile[x]='.';
+ geofile[x]='.';
+ kmlfile[x]='.';
+ mapfile[x+1]='p';
+ geofile[x+1]='g';
+ kmlfile[x+1]='k';
+ mapfile[x+2]='p';
+ geofile[x+2]='e';
+ kmlfile[x+2]='m';
+ mapfile[x+3]='m';
+ geofile[x+3]='o';
+ kmlfile[x+3]='l';
+ mapfile[x+4]=0;
+ geofile[x+4]=0;
+ kmlfile[x+4]=0;
+
+ if (geo && kml==0)
+ {
+ fd=fopen(geofile,"wb");
- case 10:
- fprintf(fd,"%c%c%c",80,80,255);
- break;
+ fprintf(fd,"FILENAME\t%s\n",mapfile);
+ fprintf(fd,"#\t\tX\tY\tLong\t\tLat\n");
+ fprintf(fd,"TIEPOINT\t0\t0\t%d.000\t\t%d.000\n",(max_west<180?-max_west:360-max_west),max_north);
+ fprintf(fd,"TIEPOINT\t%u\t%u\t%d.000\t\t%.3f\n",width-1,height+29,(min_west<180?-min_west:360-min_west),(double)(min_north-0.025));
+ fprintf(fd,"IMAGESIZE\t%u\t%u\n",width,height+30);
+ fprintf(fd,"#\n# Auto Generated by SPLAT! v%s\n#\n",splat_version);
- case 11:
- fprintf(fd,"%c%c%c",0,38,255);
- break;
+ fclose(fd);
+ }
- case 12:
- fprintf(fd,"%c%c%c",142,63,255);
- break;
+ if (kml && geo==0)
+ {
+ fd=fopen(kmlfile,"wb");
+
+ fprintf(fd,"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");
+ fprintf(fd,"<kml xmlns=\"http://earth.google.com/kml/2.1\">\n");
+ fprintf(fd,"<!-- Generated by SPLAT! Version %s -->\n",splat_version);
+ fprintf(fd," <Folder>\n");
+ fprintf(fd," <name>SPLAT!</name>\n");
+ fprintf(fd," <description>%s Transmitter Coverage Overlay</description>\n",xmtr[0].name);
+ fprintf(fd," <GroundOverlay>\n");
+ fprintf(fd," <name>SPLAT! Signal Strength Overlay</name>\n");
+ fprintf(fd," <description>SPLAT! Coverage</description>\n");
+ fprintf(fd," <Icon>\n");
+ fprintf(fd," <href>%s</href>\n",mapfile);
+ fprintf(fd," </Icon>\n");
+ fprintf(fd," <opacity>128</opacity>\n");
+ fprintf(fd," <LatLonBox>\n");
+ fprintf(fd," <north>%.5f</north>\n",(double)max_north-one_pixel);
+ fprintf(fd," <south>%.5f</south>\n",(double)min_north);
+ fprintf(fd," <east>%.5f</east>\n",((double)min_west<180.0?(double)-min_west:360.0-(double)min_west));
+ fprintf(fd," <west>%.5f</west>\n",(((double)max_west-one_pixel)<180.0?-((double)max_west-one_pixel):(360.0-(double)max_west-one_pixel)));
+ fprintf(fd," <rotation>0.0</rotation>\n");
+ fprintf(fd," </LatLonBox>\n");
+ fprintf(fd," </GroundOverlay>\n");
- case 13:
- fprintf(fd,"%c%c%c",196,54,255);
- break;
+ for (x=0; x<txsites; x++)
+ {
+ fprintf(fd," <Placemark>\n");
+ fprintf(fd," <name>%s</name>\n",xmtr[x].name);
+ fprintf(fd," <visibility>1</visibility>\n");
+ fprintf(fd," <Style>\n");
+ fprintf(fd," <IconStyle>\n");
+ fprintf(fd," <Icon>\n");
+ fprintf(fd," <href>root://icons/palette-5.png</href>\n");
+ fprintf(fd," <x>224</x>\n");
+ fprintf(fd," <y>224</y>\n");
+ fprintf(fd," <w>32</w>\n");
+ fprintf(fd," <h>32</h>\n");
+ fprintf(fd," </Icon>\n");
+ fprintf(fd," </IconStyle>\n");
+ fprintf(fd," </Style>\n");
+ fprintf(fd," <Point>\n");
+ fprintf(fd," <extrude>1</extrude>\n");
+ fprintf(fd," <altitudeMode>relativeToGround</altitudeMode>\n");
+ fprintf(fd," <coordinates>%f,%f,%f</coordinates>\n",(xmtr[x].lon<180.0?-xmtr[x].lon:360.0-xmtr[x].lon), xmtr[x].lat, xmtr[x].alt);
+ fprintf(fd," </Point>\n");
+ fprintf(fd," </Placemark>\n");
+ }
- case 14:
- fprintf(fd,"%c%c%c",255,0,255);
- break;
+ fprintf(fd," </Folder>\n");
+ fprintf(fd,"</kml>\n");
+
+ fclose(fd);
+ }
+
+ fd=fopen(mapfile,"wb");
+
+ fprintf(fd,"P6\n%u %u\n255\n",width,(kml?height:height+30));
+ fprintf(stdout,"\nWriting \"%s\" (%ux%u pixmap image)... ",mapfile,width,(kml?height:height+30));
+ fflush(stdout);
+
+ for (y=0, lat=(double)max_north-one_pixel; y<(int)height; y++, lat=(double)max_north-(one_pixel*(double)y))
+ {
+ for (x=0, lon=(double)max_west-one_pixel; x<(int)width; x++, lon=(double)max_west-(one_pixel*(double)x))
+ {
+ if (lon<0.0)
+ lon+=360.0;
+
+ for (indx=0, found=0; indx<MAXPAGES && found==0;)
+ if (lat>=(double)dem[indx].min_north && lat<(double)dem[indx].max_north && LonDiff(lon,(double)dem[indx].min_west)>=0.0 && LonDiff(lon,(double)dem[indx].max_west)<0.0)
+ found=1;
+ else
+ indx++;
+
+ if (found)
+ {
+ x0=(int)(1199.0*(lat-floor(lat)));
+ y0=(int)(1199.0*(lon-floor(lon)));
+
+ mask=dem[indx].mask[x0][y0];
+ signal=(dem[indx].signal[x0][y0])-100;
+ cityorcounty=0;
+
+ match=255;
+
+ red=0;
+ green=0;
+ blue=0;
+
+ if (signal>=region.level[0])
+ match=0;
+ else
+ {
+ for (z=1; (z<region.levels && match==255); z++)
+ {
+ if (signal<region.level[z-1] && signal>=region.level[z])
+ match=z;
+ }
+ }
+
+ if (match<region.levels)
+ {
+ red=region.color[match][0];
+ green=region.color[match][1];
+ blue=region.color[match][2];
+
+ color=1;
+ }
+
+ if ((mask&2) && (kml==0))
+ {
+ /* Text Labels: Red or otherwise */
+
+ if (red>=180 && green<=75 && blue<=75)
+ fprintf(fd,"%c%c%c",255^red,255^green,255^blue);
+ else
+ fprintf(fd,"%c%c%c",255,0,0);
+
+ cityorcounty=1;
+ }
+
+ else if ((mask&4) && (kml==0))
+ {
+ /* County Boundaries: Black */
- case 255:
- /* Black */
fprintf(fd,"%c%c%c",0,0,0);
- break;
- default:
- fprintf(fd,"%c%c%c",255,194,204);
+ cityorcounty=1;
+ }
+
+ if (cityorcounty==0)
+ {
+ if (dem[indx].signal[x0][y0]==0)
+ {
+ if (ngs)
+ fprintf(fd,"%c%c%c",255,255,255);
+ else
+ {
+ /* Display land or sea elevation */
+
+ if (dem[indx].data[x0][y0]==0)
+ fprintf(fd,"%c%c%c",0,0,170);
+ else
+ {
+ terrain=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
+ fprintf(fd,"%c%c%c",terrain,terrain,terrain);
+ }
+ }
+ }
+
+ else
+ {
+ /* Plot field strength regions in color */
+
+ if (red!=0 || green!=0 || blue!=0)
+ fprintf(fd,"%c%c%c",red,green,blue);
+
+ else /* terrain / sea-level */
+ {
+ if (ngs)
+ fprintf(fd,"%c%c%c",255,255,255);
+ else
+ {
+ if (dem[indx].data[x0][y0]==0)
+ fprintf(fd,"%c%c%c",0,0,170);
+ else
+ {
+ /* Elevation: Greyscale */
+ terrain=(unsigned)(0.5+pow((double)(dem[indx].data[x0][y0]-min_elevation),one_over_gamma)*conversion);
+ fprintf(fd,"%c%c%c",terrain,terrain,terrain);
+ }
+ }
+ }
+ }
+ }
+ }
+
+ else
+ {
+ /* We should never get here, but if */
+ /* we do, display the region as black */
+
+ fprintf(fd,"%c%c%c",0,0,0);
+ }
+ }
+ }
+
+ if (kml==0 && color)
+ {
+ /* Display legend along bottom of image */
+
+ colorwidth=(int)rint((float)width/(float)region.levels);
+
+ for (y0=0; y0<30; y0++)
+ {
+ for (x0=0; x0<(int)width; x0++)
+ {
+ indx=x0/colorwidth;
+ x=x0%colorwidth;
+ level=region.level[indx];
+
+ hundreds=level/100;
+
+ if (hundreds>0)
+ level-=(hundreds*100);
+
+ tens=level/10;
+
+ if (tens>0)
+ level-=(tens*10);
+
+ units=level;
+
+ if (y0>=8 && y0<=23)
+ {
+ if (hundreds>0)
+ {
+ if (x>=5 && x<=12)
+ if (fontdata[16*(hundreds+'0')+(y0-8)]&(128>>(x-5)))
+ indx=255;
+ }
+
+ if (tens>0 || hundreds>0)
+ {
+ if (x>=13 && x<=20)
+ if (fontdata[16*(tens+'0')+(y0-8)]&(128>>(x-13)))
+ indx=255;
+ }
+
+ if (x>=21 && x<=28)
+ if (fontdata[16*(units+'0')+(y0-8)]&(128>>(x-21)))
+ indx=255;
+
+ if (x>=36 && x<=43)
+ if (fontdata[16*('d')+(y0-8)]&(128>>(x-36)))
+ indx=255;
+
+ if (x>=44 && x<=51)
+ if (fontdata[16*('B')+(y0-8)]&(128>>(x-44)))
+ indx=255;
+
+ if (x>=52 && x<=59)
+ if (fontdata[16*('u')+(y0-8)]&(128>>(x-52)))
+ indx=255;
+
+ if (x>=60 && x<=67)
+ if (fontdata[16*('V')+(y0-8)]&(128>>(x-60)))
+ indx=255;
+
+ if (x>=68 && x<=75)
+ if (fontdata[16*('/')+(y0-8)]&(128>>(x-68)))
+ indx=255;
+
+ if (x>=76 && x<=83)
+ if (fontdata[16*('m')+(y0-8)]&(128>>(x-76)))
+ indx=255;
+ }
+
+ if (indx>region.levels)
+ fprintf(fd,"%c%c%c",0,0,0);
+ else
+ {
+ red=region.color[indx][0];
+ green=region.color[indx][1];
+ blue=region.color[indx][2];
+
+ fprintf(fd,"%c%c%c",red,green,blue);
}
}
}
{
if (metric)
fprintf(fd,"%f\t%f\n",KM_PER_MILE*path.distance[x],METERS_PER_FOOT*path.elevation[x]);
-
else
fprintf(fd,"%f\t%f\n",path.distance[x],path.elevation[x]);
}
else if (strncmp(ext,"ps",2)==0)
strncpy(term,"postscript enhanced color\0",26);
- fprintf(stdout,"Writing \"%s.%s\"...",filename,ext);
- fflush(stdout);
-
fd=fopen("splat.gp","w");
fprintf(fd,"set grid\n");
fprintf(fd,"set autoscale\n");
{
unlink("splat.gp");
unlink("profile.gp");
- fprintf(stdout," Done!\n");
+
+ fprintf(stdout,"\nTerrain plot written to: \"%s.%s\"",filename,ext);
fflush(stdout);
}
else if (strncmp(ext,"ps",2)==0)
strncpy(term,"postscript enhanced color\0",26);
- fprintf(stdout,"Writing \"%s.%s\"...",filename,ext);
- fflush(stdout);
-
fd=fopen("splat.gp","w");
fprintf(fd,"set grid\n");
unlink("profile.gp");
unlink("reference.gp");
- fprintf(stdout," Done!\n");
+ fprintf(stdout,"\nElevation plot written to: \"%s.%s\"",filename,ext);
fflush(stdout);
}
fd4=fopen("fresnel_pt_6.gp", "wb");
}
- for (x=0; x<path.length; x++)
+ for (x=0; x<path.length-1; x++)
{
remote.lat=path.lat[x];
remote.lon=path.lon[x];
if (f)
{
d1=5280.0*path.distance[x];
- f_zone=-1*sqrt(lambda*d1*(d-d1)/d);
- fpt6_zone=0.6*f_zone;
+ f_zone=-1.0*sqrt(lambda*d1*(d-d1)/d);
+ fpt6_zone=f_zone*fzone_clearance;
}
if (n)
else if (strncmp(ext,"ps",2)==0)
strncpy(term,"postscript enhanced color\0",26);
- fprintf(stdout,"Writing \"%s.%s\"...",filename,ext);
- fflush(stdout);
-
fd=fopen("splat.gp","w");
dheight=maxheight-minheight;
fprintf(fd,"set output \"%s.%s\"\n",filename,ext);
if (f)
- fprintf(fd,"plot \"profile.gp\" title \"Point-to-Point Profile\" with lines, \"reference.gp\" title \"Line of Sight Path\" with lines, \"curvature.gp\" axes x1y2 title \"Earth's Curvature Contour\" with lines, \"fresnel.gp\" axes x1y1 title \"First Fresnel Zone (%.3f MHz)\" with lines, \"fresnel_pt_6.gp\" title \"60%% of First Fresnel Zone\" with lines\n",f);
+ fprintf(fd,"plot \"profile.gp\" title \"Point-to-Point Profile\" with lines, \"reference.gp\" title \"Line of Sight Path\" with lines, \"curvature.gp\" axes x1y2 title \"Earth's Curvature Contour\" with lines, \"fresnel.gp\" axes x1y1 title \"First Fresnel Zone (%.3f MHz)\" with lines, \"fresnel_pt_6.gp\" title \"%.0f%% of First Fresnel Zone\" with lines\n",f,fzone_clearance*100.0);
else
fprintf(fd,"plot \"profile.gp\" title \"Point-to-Point Profile\" with lines, \"reference.gp\" title \"Line Of Sight Path\" with lines, \"curvature.gp\" axes x1y2 title \"Earth's Curvature Contour\" with lines\n");
unlink("fresnel_pt_6.gp");
}
- fprintf(stdout," Done!\n");
+ fprintf(stdout,"\nHeight plot written to: \"%s.%s\"",filename,ext);
fflush(stdout);
}
fprintf(stderr,"\n*** ERROR: Error occurred invoking gnuplot!\n");
}
-void GraphLongley(struct site source, struct site destination, char *name)
+void ObstructionAnalysis(struct site xmtr, struct site rcvr, double f, FILE *outfile)
{
- /* This function invokes gnuplot to generate an appropriate
- output file indicating the Longley-Rice model loss between
- the source and destination locations. "filename" is the
- name assigned to the output file generated by gnuplot.
- The filename extension is used to set gnuplot's terminal
- setting and output file type. If no extension is found,
- .png is assumed. */
-
- int x, y, z, errnum, errflag=0;
+ /* Perform an obstruction analysis along the
+ path between receiver and transmitter. */
+
+ int x;
+ struct site site_x;
+ double h_r, h_t, h_x, h_r_orig, cos_tx_angle, cos_test_angle,
+ cos_tx_angle_f1, cos_tx_angle_fpt6, d_tx, d_x,
+ h_r_f1, h_r_fpt6, h_f, h_los, lambda=0.0;
+ char string[255], string_fpt6[255], string_f1[255];
+
+ ReadPath(xmtr,rcvr);
+ h_r=GetElevation(rcvr)+rcvr.alt+earthradius;
+ h_r_f1=h_r;
+ h_r_fpt6=h_r;
+ h_r_orig=h_r;
+ h_t=GetElevation(xmtr)+xmtr.alt+earthradius;
+ d_tx=5280.0*Distance(rcvr,xmtr);
+ cos_tx_angle=((h_r*h_r)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r*d_tx);
+ cos_tx_angle_f1=cos_tx_angle;
+ cos_tx_angle_fpt6=cos_tx_angle;
+
+ if (f)
+ lambda=9.8425e8/(f*1e6);
+
+ /* At each point along the path calculate the cosine
+ of a sort of "inverse elevation angle" at the receiver.
+ From the antenna, 0 deg. looks at the ground, and 90 deg.
+ is parallel to the ground.
+
+ Start at the receiver. If this is the lowest antenna,
+ then terrain obstructions will be nearest to it. (Plus,
+ that's the way SPLAT!'s original los() did it.)
+
+ Calculate cosines only. That's sufficient to compare
+ angles and it saves the extra computational burden of
+ acos(). However, note the inverted comparison: if
+ acos(A) > acos(B), then B > A. */
+
+ for (x=path.length-1; x>0; x--)
+ {
+ site_x.lat=path.lat[x];
+ site_x.lon=path.lon[x];
+ site_x.alt=0.0;
+
+ h_x=GetElevation(site_x)+earthradius;
+ d_x=5280.0*Distance(rcvr,site_x);
+
+ /* Deal with the LOS path first. */
+
+ cos_test_angle=((h_r*h_r)+(d_x*d_x)-(h_x*h_x))/(2.0*h_r*d_x);
+
+ if (cos_tx_angle>cos_test_angle)
+ {
+ if (h_r==h_r_orig)
+ fprintf(outfile,"Between %s and %s, SPLAT! detected obstructions at:\n\n",rcvr.name,xmtr.name);
+
+ if (site_x.lat>=0.0)
+ {
+ if (metric)
+ fprintf(outfile,"\t%.4f N, %.4f W, %5.2f kilometers, %6.2f meters AMSL\n",site_x.lat, site_x.lon, KM_PER_MILE*(d_x/5280.0), METERS_PER_FOOT*(h_x-earthradius));
+ else
+ fprintf(outfile,"\t%.4f N, %.4f W, %5.2f miles, %6.2f feet AMSL\n",site_x.lat, site_x.lon, d_x/5280.0, h_x-earthradius);
+ }
+
+ else
+ {
+ if (metric)
+ fprintf(outfile,"\t%.4f S, %.4f W, %5.2f kilometers, %6.2f meters AMSL\n",-site_x.lat, site_x.lon, KM_PER_MILE*(d_x/5280.0), METERS_PER_FOOT*(h_x-earthradius));
+ else
+
+ fprintf(outfile,"\t%.4f S, %.4f W, %5.2f miles, %6.2f feet AMSL\n",-site_x.lat, site_x.lon, d_x/5280.0, h_x-earthradius);
+ }
+ }
+
+ while (cos_tx_angle>cos_test_angle)
+ {
+ h_r+=1;
+ cos_test_angle=((h_r*h_r)+(d_x*d_x)-(h_x*h_x))/(2.0*h_r*d_x);
+ cos_tx_angle=((h_r*h_r)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r*d_tx);
+ }
+
+ if (f)
+ {
+ /* Now clear the first Fresnel zone... */
+
+ cos_tx_angle_f1=((h_r_f1*h_r_f1)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_f1*d_tx);
+ h_los=sqrt(h_r_f1*h_r_f1+d_x*d_x-2*h_r_f1*d_x*cos_tx_angle_f1);
+ h_f=h_los-sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
+
+ while (h_f<h_x)
+ {
+ h_r_f1+=1;
+ cos_tx_angle_f1=((h_r_f1*h_r_f1)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_f1*d_tx);
+ h_los=sqrt(h_r_f1*h_r_f1+d_x*d_x-2*h_r_f1*d_x*cos_tx_angle_f1);
+ h_f=h_los-sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
+ }
+
+ /* and clear the 60% F1 zone. */
+
+ cos_tx_angle_fpt6=((h_r_fpt6*h_r_fpt6)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_fpt6*d_tx);
+ h_los=sqrt(h_r_fpt6*h_r_fpt6+d_x*d_x-2*h_r_fpt6*d_x*cos_tx_angle_fpt6);
+ h_f=h_los-fzone_clearance*sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
+
+ while (h_f<h_x)
+ {
+ h_r_fpt6+=1;
+ cos_tx_angle_fpt6=((h_r_fpt6*h_r_fpt6)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_fpt6*d_tx);
+ h_los=sqrt(h_r_fpt6*h_r_fpt6+d_x*d_x-2*h_r_fpt6*d_x*cos_tx_angle_fpt6);
+ h_f=h_los-fzone_clearance*sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
+ }
+ }
+ }
+
+ if (h_r>h_r_orig)
+ {
+ if (metric)
+ sprintf(string,"\nAntenna at %s must be raised to at least %.2f meters AGL\nto clear all obstructions detected by SPLAT!\n",rcvr.name, METERS_PER_FOOT*(h_r-GetElevation(rcvr)-earthradius));
+ else
+ sprintf(string,"\nAntenna at %s must be raised to at least %.2f feet AGL\nto clear all obstructions detected by SPLAT!\n",rcvr.name, h_r-GetElevation(rcvr)-earthradius);
+ }
+
+ else
+ sprintf(string,"\nNo obstructions to LOS path due to terrain were detected by SPLAT!\n");
+
+ if (f)
+ {
+ if (h_r_fpt6>h_r_orig)
+ {
+ if (metric)
+ sprintf(string_fpt6,"\nAntenna at %s must be raised to at least %.2f meters AGL\nto clear %.0f%c of the first Fresnel zone.\n",rcvr.name, METERS_PER_FOOT*(h_r_fpt6-GetElevation(rcvr)-earthradius),fzone_clearance*100.0,37);
+
+ else
+ sprintf(string_fpt6,"\nAntenna at %s must be raised to at least %.2f feet AGL\nto clear %.0f%c of the first Fresnel zone.\n",rcvr.name, h_r_fpt6-GetElevation(rcvr)-earthradius,fzone_clearance*100.0,37);
+ }
+
+ else
+ sprintf(string_fpt6,"\n%.0f%c of the first Fresnel zone is clear.\n",fzone_clearance*100.0,37);
+
+ if (h_r_f1>h_r_orig)
+ {
+ if (metric)
+ sprintf(string_f1,"\nAntenna at %s must be raised to at least %.2f meters AGL\nto clear the first Fresnel zone.\n",rcvr.name, METERS_PER_FOOT*(h_r_f1-GetElevation(rcvr)-earthradius));
+
+ else
+ sprintf(string_f1,"\nAntenna at %s must be raised to at least %.2f feet AGL\nto clear the first Fresnel zone.\n",rcvr.name, h_r_f1-GetElevation(rcvr)-earthradius);
+
+ }
+
+ else
+ sprintf(string_f1,"\nThe first Fresnel zone is clear.\n\n");
+ }
+
+ fprintf(outfile,"%s",string);
+
+ if (f)
+ {
+ fprintf(outfile,"%s",string_f1);
+ fprintf(outfile,"%s",string_fpt6);
+ }
+}
+
+void PathReport(struct site source, struct site destination, char *name, char graph_it)
+{
+ /* This function writes a SPLAT! Path Report (name.txt) to
+ the filesystem. If (graph_it == 1), then gnuplot is invoked
+ to generate an appropriate output file indicating the Longley-Rice
+ model loss between the source and destination locations.
+ "filename" is the name assigned to the output file generated
+ by gnuplot. The filename extension is used to set gnuplot's
+ terminal setting and output file type. If no extension is
+ found, .png is assumed. */
+
+ int x, y, z, errnum;
char filename[255], term[30], ext[15], strmode[100],
report_name[80], block=0;
double maxloss=-100000.0, minloss=100000.0, loss, haavt,
angle1, angle2, azimuth, pattern=1.0, patterndB=0.0,
total_loss=0.0, cos_xmtr_angle, cos_test_angle=0.0,
source_alt, test_alt, dest_alt, source_alt2, dest_alt2,
- distance, elevation, four_thirds_earth;
+ distance, elevation, four_thirds_earth, field_strength,
+ free_space_loss=0.0, voltage;
FILE *fd=NULL, *fd2=NULL;
- sprintf(report_name,"%s-to-%s.lro",source.name,destination.name);
+ sprintf(report_name,"%s-to-%s.txt",source.name,destination.name);
four_thirds_earth=EARTHRADIUS*(4.0/3.0);
fd2=fopen(report_name,"w");
- fprintf(fd2,"\n\t--==[ SPLAT! v%s Longley-Rice Model Path Loss Report ]==--\n\n",splat_version);
- fprintf(fd2,"Analysis of RF path conditions between %s and %s:\n",source.name, destination.name);
- fprintf(fd2,"\n-------------------------------------------------------------------------\n\n");
+ fprintf(fd2,"\n\t\t--==[ SPLAT! v%s Path Analysis ]==--\n\n",splat_version);
+ fprintf(fd2,"-------------------------------------------------------------------------\n\n");
fprintf(fd2,"Transmitter site: %s\n",source.name);
if (source.lat>=0.0)
pattern=(double)LR.antenna_pattern[(int)rint(azimuth)][x];
patterndB=20.0*log10(pattern);
-
- fprintf(fd2,"Antenna pattern between %s and %s: %.3f (%.2f dB)\n",source.name, destination.name, pattern, patterndB);
-
}
if (metric)
else
fprintf(fd2,"Depression angle to %s: %+.4f degrees\n",destination.name,angle1);
- if (angle1!=angle2)
+ if ((angle2-angle1)>0.0001)
{
if (angle2<0.0)
fprintf(fd2,"Depression");
else
fprintf(fd2,"Depression angle to %s: %+.4f degrees\n",source.name,angle1);
- if (angle1!=angle2)
+ if ((angle2-angle1)>0.0001)
{
if (angle2<0.0)
fprintf(fd2,"Depression");
fprintf(fd2,"\n-------------------------------------------------------------------------\n\n");
- fprintf(fd2,"Longley-Rice path calculation parameters used in this analysis:\n\n");
- fprintf(fd2,"Earth's Dielectric Constant: %.3lf\n",LR.eps_dielect);
- fprintf(fd2,"Earth's Conductivity: %.3lf\n",LR.sgm_conductivity);
- fprintf(fd2,"Atmospheric Bending Constant (N): %.3lf\n",LR.eno_ns_surfref);
- fprintf(fd2,"Frequency: %.3lf (MHz)\n",LR.frq_mhz);
- fprintf(fd2,"Radio Climate: %d (",LR.radio_climate);
-
- switch (LR.radio_climate)
+ if (LR.frq_mhz>0.0)
{
- case 1:
- fprintf(fd2,"Equatorial");
- break;
+ fprintf(fd2,"Longley-Rice path calculation parameters used in this analysis:\n\n");
+ fprintf(fd2,"Earth's Dielectric Constant: %.3lf\n",LR.eps_dielect);
+ fprintf(fd2,"Earth's Conductivity: %.3lf Siemens/meter\n",LR.sgm_conductivity);
+ fprintf(fd2,"Atmospheric Bending Constant (N-units): %.3lf ppm\n",LR.eno_ns_surfref);
+ fprintf(fd2,"Frequency: %.3lf MHz\n",LR.frq_mhz);
+ fprintf(fd2,"Radio Climate: %d (",LR.radio_climate);
- case 2:
- fprintf(fd2,"Continental Subtropical");
- break;
+ switch (LR.radio_climate)
+ {
+ case 1:
+ fprintf(fd2,"Equatorial");
+ break;
- case 3:
- fprintf(fd2,"Maritime Subtropical");
- break;
+ case 2:
+ fprintf(fd2,"Continental Subtropical");
+ break;
- case 4:
- fprintf(fd2,"Desert");
- break;
+ case 3:
+ fprintf(fd2,"Maritime Subtropical");
+ break;
- case 5:
- fprintf(fd2,"Continental Temperate");
- break;
+ case 4:
+ fprintf(fd2,"Desert");
+ break;
- case 6:
- fprintf(fd2,"Martitime Temperate, Over Land");
- break;
+ case 5:
+ fprintf(fd2,"Continental Temperate");
+ break;
- case 7:
- fprintf(fd2,"Maritime Temperate, Over Sea");
- break;
+ case 6:
+ fprintf(fd2,"Martitime Temperate, Over Land");
+ break;
- default:
- fprintf(fd2,"Unknown");
- }
-
- fprintf(fd2,")\nPolarization: %d (",LR.pol);
-
- if (LR.pol==0)
- fprintf(fd2,"Horizontal");
-
- if (LR.pol==1)
- fprintf(fd2,"Vertical");
-
- fprintf(fd2,")\nFraction of Situations: %.1lf%c\n",LR.conf*100.0,37);
- fprintf(fd2,"Fraction of Time: %.1lf%c\n",LR.rel*100.0,37);
- fprintf(fd2,"\n-------------------------------------------------------------------------\n\n");
-
- fprintf(fd2,"Analysis Results:\n\n");
-
- ReadPath(source, destination); /* source=TX, destination=RX */
-
- /* Copy elevations along path into the elev_l[] array. */
-
- for (x=0; x<path.length; x++)
- elev_l[x+2]=path.elevation[x]*METERS_PER_FOOT;
-
- if (metric)
- fprintf(fd2,"Distance (km)");
- else
- fprintf(fd2,"Distance (mi)");
-
- fprintf(fd2,"\tLoss (dB)\tErrnum\tComment\n\n");
-
- fd=fopen("profile.gp","w");
-
- azimuth=rint(Azimuth(source,destination));
-
- for (y=2; y<(path.length-1); y++) /* path.length-1 avoids LR error */
- {
- distance=5280.0*path.distance[y];
- source_alt=four_thirds_earth+source.alt+path.elevation[0];
- dest_alt=four_thirds_earth+destination.alt+path.elevation[y];
- dest_alt2=dest_alt*dest_alt;
- source_alt2=source_alt*source_alt;
-
- /* Calculate the cosine of the elevation of
- the receiver as seen by the transmitter. */
-
- cos_xmtr_angle=((source_alt2)+(distance*distance)-(dest_alt2))/(2.0*source_alt*distance);
-
- if (got_elevation_pattern)
- {
- /* If an antenna elevation pattern is available, the
- following code determines the elevation angle to
- the first obstruction along the path. */
-
- for (x=2, block=0; x<y && block==0; x++)
- {
- distance=5280.0*(path.distance[y]-path.distance[x]);
- test_alt=four_thirds_earth+path.elevation[x];
-
- /* Calculate the cosine of the elevation
- angle of the terrain (test point)
- as seen by the transmitter. */
-
- cos_test_angle=((source_alt2)+(distance*distance)-(test_alt*test_alt))/(2.0*source_alt*distance);
-
- /* Compare these two angles to determine if
- an obstruction exists. Since we're comparing
- the cosines of these angles rather than
- the angles themselves, the sense of the
- following "if" statement is reversed from
- what it would be if the angles themselves
- were compared. */
-
- if (cos_xmtr_angle>cos_test_angle)
- block=1;
- }
+ case 7:
+ fprintf(fd2,"Maritime Temperate, Over Sea");
+ break;
- /* At this point, we have the elevation angle
- to the first obstruction (if it exists). */
+ default:
+ fprintf(fd2,"Unknown");
}
- /* Determine path loss for each point along the
- path using Longley-Rice's point_to_point mode
- starting at x=2 (number_of_points = 1), the
- shortest distance terrain can play a role in
- path loss. */
-
- elev_l[0]=y-1; /* (number of points - 1) */
-
- /* Distance between elevation samples */
- elev_l[1]=METERS_PER_MILE*(path.distance[y]-path.distance[y-1]);
+ fprintf(fd2,")\nPolarization: %d (",LR.pol);
- point_to_point(elev_l, source.alt*METERS_PER_FOOT,
- destination.alt*METERS_PER_FOOT, LR.eps_dielect,
- LR.sgm_conductivity, LR.eno_ns_surfref, LR.frq_mhz,
- LR.radio_climate, LR.pol, LR.conf, LR.rel, loss,
- strmode, errnum);
+ if (LR.pol==0)
+ fprintf(fd2,"Horizontal");
- if (block)
- elevation=((acos(cos_test_angle))/deg2rad)-90.0;
- else
- elevation=((acos(cos_xmtr_angle))/deg2rad)-90.0;
-
- /* Integrate the antenna's radiation
- pattern into the overall path loss. */
+ if (LR.pol==1)
+ fprintf(fd2,"Vertical");
- x=(int)rint(10.0*(10.0-elevation));
+ fprintf(fd2,")\nFraction of Situations: %.1lf%c\n",LR.conf*100.0,37);
+ fprintf(fd2,"Fraction of Time: %.1lf%c\n",LR.rel*100.0,37);
- if (x>=0 && x<=1000)
+ if (LR.erp!=0.0)
{
- pattern=(double)LR.antenna_pattern[(int)azimuth][x];
+ fprintf(fd2,"Transmitter ERP: ");
- if (pattern!=0.0)
- patterndB=20.0*log10(pattern);
- }
+ if (LR.erp<1.0)
+ fprintf(fd2,"%.1lf milliwatts\n",1000.0*LR.erp);
- else
- patterndB=0.0;
+ if (LR.erp>=1.0 && LR.erp<10.0)
+ fprintf(fd2,"%.1lf Watts\n",LR.erp);
- total_loss=loss-patterndB;
+ if (LR.erp>=10.0 && LR.erp<10.0e3)
+ fprintf(fd2,"%.0lf Watts\n",LR.erp);
- if (metric)
- {
- fprintf(fd,"%f\t%f\n",KM_PER_MILE*(path.distance[path.length-1]-path.distance[y]),total_loss);
- fprintf(fd2,"%7.2f\t\t%7.2f\t\t %d\t%s\n",KM_PER_MILE*path.distance[y],total_loss, errnum, strmode);
+ if (LR.erp>=10.0e3)
+ fprintf(fd2,"%.3lf kilowatts\n",LR.erp/1.0e3);
}
- else
- {
- fprintf(fd,"%f\t%f\n",path.distance[path.length-1]-path.distance[y],total_loss);
- fprintf(fd2,"%7.2f\t\t%7.2f\t\t %d\t%s\n",path.distance[y],total_loss, errnum, strmode);
- }
+ fprintf(fd2,"\n-------------------------------------------------------------------------\n\n");
- errflag|=errnum;
+ fprintf(fd2,"Summary for the link between %s and %s:\n\n",source.name, destination.name);
- if (total_loss>maxloss)
- maxloss=total_loss;
+ if (patterndB!=0.0)
+ fprintf(fd2,"%s antenna pattern towards %s: %.3f (%.2f dB)\n", source.name, destination.name, pattern, patterndB);
- if (total_loss<minloss)
- minloss=total_loss;
- }
-
- fclose(fd);
-
- fprintf(fd2,"\nLongley-Rice path loss between %s and %s is %.2f dB.\n",source.name, destination.name, loss);
-
- if (patterndB!=0.0)
- fprintf(fd2,"Total path loss including TX antenna pattern is %.2f dB.\n",total_loss);
+ ReadPath(source, destination); /* source=TX, destination=RX */
- if (errflag)
- {
- fprintf(fd2,"\nNotes on \"errnum\":\n\n");
- fprintf(fd2," 0: No error. :-)\n");
- fprintf(fd2," 1: Warning! Some parameters are nearly out of range.\n");
- fprintf(fd2," Results should be used with caution.\n");
- fprintf(fd2," 2: Note: Default parameters have been substituted for impossible ones.\n");
- fprintf(fd2," 3: Warning! A combination of parameters is out of range.\n");
- fprintf(fd2," Results are probably invalid.\n");
- fprintf(fd2," Other: Warning! Some parameters are out of range.\n");
- fprintf(fd2," Results are probably invalid.\n\nEnd of Report\n");
- }
+ /* Copy elevations along path into the elev_l[] array. */
- fprintf(stdout,"Longley-Rice path loss between %s and %s is %.2f dB.\n",source.name, destination.name, loss);
+ for (x=0; x<path.length; x++)
+ elev_l[x+2]=path.elevation[x]*METERS_PER_FOOT;
- if (patterndB!=0.0)
- fprintf(stdout,"Total path loss including TX antenna pattern is %.2f dB.\n",total_loss);
+ fd=fopen("profile.gp","w");
- fprintf(stdout,"Path Loss Report written to: \"%s\"\n",report_name);
- fflush(stdout);
-
- fclose(fd2);
+ azimuth=rint(Azimuth(source,destination));
- if (name[0]==0)
- {
- /* Default filename and output file type */
-
- strncpy(filename,"loss\0",5);
- strncpy(term,"png\0",4);
- strncpy(ext,"png\0",4);
- }
-
- else
- {
- /* Grab extension and terminal type from "name" */
-
- for (x=0; name[x]!='.' && name[x]!=0 && x<254; x++)
- filename[x]=name[x];
-
- if (name[x]=='.')
+ for (y=2; y<(path.length-1); y++) /* path.length-1 avoids LR error */
{
- for (y=0, z=x, x++; name[x]!=0 && x<254 && y<14; x++, y++)
- {
- term[y]=tolower(name[x]);
- ext[y]=term[y];
- }
+ distance=5280.0*path.distance[y];
+ source_alt=four_thirds_earth+source.alt+path.elevation[0];
+ dest_alt=four_thirds_earth+destination.alt+path.elevation[y];
+ dest_alt2=dest_alt*dest_alt;
+ source_alt2=source_alt*source_alt;
- ext[y]=0;
- term[y]=0;
- filename[z]=0;
- }
+ /* Calculate the cosine of the elevation of
+ the receiver as seen by the transmitter. */
- else
- { /* No extension -- Default is png */
+ cos_xmtr_angle=((source_alt2)+(distance*distance)-(dest_alt2))/(2.0*source_alt*distance);
- filename[x]=0;
- strncpy(term,"png\0",4);
- strncpy(ext,"png\0",4);
- }
- }
+ if (got_elevation_pattern)
+ {
+ /* If an antenna elevation pattern is available, the
+ following code determines the elevation angle to
+ the first obstruction along the path. */
- /* Either .ps or .postscript may be used
- as an extension for postscript output. */
+ for (x=2, block=0; x<y && block==0; x++)
+ {
+ distance=5280.0*(path.distance[y]-path.distance[x]);
+ test_alt=four_thirds_earth+path.elevation[x];
- if (strncmp(term,"postscript",10)==0)
- strncpy(ext,"ps\0",3);
+ /* Calculate the cosine of the elevation
+ angle of the terrain (test point)
+ as seen by the transmitter. */
- else if (strncmp(ext,"ps",2)==0)
- strncpy(term,"postscript enhanced color\0",26);
+ cos_test_angle=((source_alt2)+(distance*distance)-(test_alt*test_alt))/(2.0*source_alt*distance);
- fprintf(stdout,"Writing \"%s.%s\"...",filename,ext);
- fflush(stdout);
+ /* Compare these two angles to determine if
+ an obstruction exists. Since we're comparing
+ the cosines of these angles rather than
+ the angles themselves, the sense of the
+ following "if" statement is reversed from
+ what it would be if the angles themselves
+ were compared. */
- fd=fopen("splat.gp","w");
+ if (cos_xmtr_angle>cos_test_angle)
+ block=1;
+ }
- fprintf(fd,"set grid\n");
- fprintf(fd,"set yrange [%2.3f to %2.3f]\n", minloss, maxloss);
- fprintf(fd,"set encoding iso_8859_1\n");
- fprintf(fd,"set term %s\n",term);
- fprintf(fd,"set title \"SPLAT! Loss Profile Along Path Between %s and %s (%.2f%c azimuth)\"\n",destination.name, source.name, Azimuth(destination,source),176);
+ /* At this point, we have the elevation angle
+ to the first obstruction (if it exists). */
+ }
- if (metric)
- fprintf(fd,"set xlabel \"Distance Between %s and %s (%.2f kilometers)\"\n",destination.name,source.name,KM_PER_MILE*Distance(destination,source));
- else
- fprintf(fd,"set xlabel \"Distance Between %s and %s (%.2f miles)\"\n",destination.name,source.name,Distance(destination,source));
+ /* Determine path loss for each point along the
+ path using Longley-Rice's point_to_point mode
+ starting at x=2 (number_of_points = 1), the
+ shortest distance terrain can play a role in
+ path loss. */
- if (got_azimuth_pattern || got_elevation_pattern)
- fprintf(fd,"set ylabel \"Total Path Loss (including TX antenna pattern) (dB)");
- else
- fprintf(fd,"set ylabel \"Longley-Rice Path Loss (dB)");
+ elev_l[0]=y-1; /* (number of points - 1) */
- fprintf(fd,"\"\nset output \"%s.%s\"\n",filename,ext);
- fprintf(fd,"plot \"profile.gp\" title \"Path Loss\" with lines\n");
+ /* Distance between elevation samples */
+ elev_l[1]=METERS_PER_MILE*(path.distance[y]-path.distance[y-1]);
- fclose(fd);
-
- x=system("gnuplot splat.gp");
+ point_to_point(elev_l, source.alt*METERS_PER_FOOT,
+ destination.alt*METERS_PER_FOOT, LR.eps_dielect,
+ LR.sgm_conductivity, LR.eno_ns_surfref, LR.frq_mhz,
+ LR.radio_climate, LR.pol, LR.conf, LR.rel, loss,
+ strmode, errnum);
- if (x!=-1)
- {
- unlink("splat.gp");
- unlink("profile.gp");
- unlink("reference.gp");
+ if (block)
+ elevation=((acos(cos_test_angle))/deg2rad)-90.0;
+ else
+ elevation=((acos(cos_xmtr_angle))/deg2rad)-90.0;
- fprintf(stdout," Done!\n");
- fflush(stdout);
- }
+ /* Integrate the antenna's radiation
+ pattern into the overall path loss. */
- else
- fprintf(stderr,"\n*** ERROR: Error occurred invoking gnuplot!\n");
-}
+ x=(int)rint(10.0*(10.0-elevation));
-void ObstructionReport(struct site xmtr, struct site rcvr, char report, double f)
-{
- int x;
- struct site site_x;
- double h_r, h_t, h_x, h_r_orig, cos_tx_angle, cos_test_angle,
- cos_tx_angle_f1, cos_tx_angle_fpt6, haavt, d_tx, d_x,
- h_r_f1, h_r_fpt6, h_f, h_los, lambda=0.0, azimuth,
- pattern, patterndB, distance, angle1, angle2;
- char report_name[80], string[255], string_fpt6[255],
- string_f1[255];
- FILE *fd;
+ if (x>=0 && x<=1000)
+ {
+ pattern=(double)LR.antenna_pattern[(int)azimuth][x];
- sprintf(report_name,"%s-to-%s.txt",xmtr.name,rcvr.name);
+ if (pattern!=0.0)
+ patterndB=20.0*log10(pattern);
+ }
- for (x=0; report_name[x]!=0; x++)
- if (report_name[x]==32 || report_name[x]==17 || report_name[x]==92 || report_name[x]==42 || report_name[x]==47)
- report_name[x]='_';
+ else
+ patterndB=0.0;
- fd=fopen(report_name,"w");
+ total_loss=loss-patterndB;
- fprintf(fd,"\n\t\t--==[ SPLAT! v%s Obstruction Report ]==--\n\n",splat_version);
- fprintf(fd,"Analysis of great circle path between %s and %s:\n",xmtr.name, rcvr.name);
- fprintf(fd,"\n-------------------------------------------------------------------------\n\n");
- fprintf(fd,"Transmitter site: %s\n",xmtr.name);
+ if (metric)
+ fprintf(fd,"%f\t%f\n",KM_PER_MILE*(path.distance[path.length-1]-path.distance[y]),total_loss);
+ else
+ fprintf(fd,"%f\t%f\n",path.distance[path.length-1]-path.distance[y],total_loss);
- if (xmtr.lat>=0.0)
- {
- fprintf(fd,"Site location: %.4f North / %.4f West",xmtr.lat, xmtr.lon);
- fprintf(fd, " (%s N / ", dec2dms(xmtr.lat));
- }
+ if (total_loss>maxloss)
+ maxloss=total_loss;
- else
- {
- fprintf(fd,"Site location: %.4f South / %.4f West",-xmtr.lat, xmtr.lon);
- fprintf(fd, " (%s S / ", dec2dms(xmtr.lat));
- }
+ if (total_loss<minloss)
+ minloss=total_loss;
+ }
- fprintf(fd, "%s W)\n", dec2dms(xmtr.lon));
+ fclose(fd);
- if (metric)
- {
- fprintf(fd,"Ground elevation: %.2f meters AMSL\n",METERS_PER_FOOT*GetElevation(xmtr));
- fprintf(fd,"Antenna height: %.2f meters AGL / %.2f meters AMSL\n",METERS_PER_FOOT*xmtr.alt, METERS_PER_FOOT*(xmtr.alt+GetElevation(xmtr)));
- }
+ distance=Distance(source,destination);
- else
- {
- fprintf(fd,"Ground elevation: %.2f feet AMSL\n",GetElevation(xmtr));
- fprintf(fd,"Antenna height: %.2f feet AGL / %.2f feet AMSL\n",xmtr.alt, xmtr.alt+GetElevation(xmtr));
- }
- haavt=haat(xmtr);
+ if (distance!=0.0)
+ {
+ free_space_loss=36.6+(20.0*log10(LR.frq_mhz))+(20.0*log10(distance));
- if (haavt>-4999.0)
- {
- if (metric)
- fprintf(fd,"Antenna height above average terrain: %.2f meters\n",METERS_PER_FOOT*haavt);
- else
- fprintf(fd,"Antenna height above average terrain: %.2f feet\n",haavt);
- }
+ fprintf(fd2,"Free space path loss: %.2f dB\n",free_space_loss);
+ }
- pattern=1.0;
- patterndB=0.0;
- distance=Distance(xmtr,rcvr);
- azimuth=Azimuth(xmtr,rcvr);
- angle1=ElevationAngle(xmtr,rcvr);
- angle2=ElevationAngle2(xmtr,rcvr,earthradius);
+ fprintf(fd2,"Longley-Rice path loss: %.2f dB\n",loss);
- if (got_azimuth_pattern || got_elevation_pattern)
- {
- x=(int)rint(10.0*(10.0-angle2));
+ if (free_space_loss!=0.0)
+ fprintf(fd2,"Attenuation due to effects of terrain: %.2f dB\n",loss-free_space_loss);
- if (x>=0 && x<=1000)
- pattern=(double)LR.antenna_pattern[(int)rint(azimuth)][x];
+ if (patterndB!=0.0)
+ fprintf(fd2,"Total path loss including %s antenna pattern: %.2f dB\n",source.name,total_loss);
- if (pattern!=1.0)
+ if (LR.erp!=0.0)
{
- fprintf(fd,"Antenna pattern toward %s: %.3f",rcvr.name,pattern);
- patterndB=20.0*log10(pattern);
- fprintf(fd," (%.2f dB)\n",patterndB);
+ field_strength=(137.26+(20.0*log10(LR.frq_mhz))-total_loss)+(10.0*log10(LR.erp/1000.0));
+ fprintf(fd2,"Field strength at %s: %.2f dBuV/meter\n", destination.name,field_strength);
+ voltage=(pow(10.0,field_strength/20.0)*39.52726907)/LR.frq_mhz;
+ fprintf(fd2,"Voltage produced by a terminated 50 ohm 0 dBd gain antenna: %.2f uV\n",voltage);
+ voltage=(pow(10.0,field_strength/20.0)*48.41082007)/LR.frq_mhz;
+ fprintf(fd2,"Voltage produced by a terminated 75 ohm 0 dBd gain antenna: %.2f uV\n",voltage);
}
- }
-
- if (metric)
- fprintf(fd,"Distance to %s: %.2f kilometers\n",rcvr.name,KM_PER_MILE*distance);
-
- else
- fprintf(fd,"Distance to %s: %.2f miles\n",rcvr.name,distance);
-
- fprintf(fd,"Azimuth to %s: %.2f degrees\n",rcvr.name,azimuth);
-
- if (angle1>=0.0)
- fprintf(fd,"Elevation angle to %s: %+.4f degrees\n",rcvr.name,angle1);
-
- else
- fprintf(fd,"Depression angle to %s: %+.4f degrees\n",rcvr.name,angle1);
-
- if (angle1!=angle2)
- {
- if (angle2<0.0)
- fprintf(fd,"Depression");
- else
- fprintf(fd,"Elevation");
-
- fprintf(fd," angle to the first obstruction: %+.4f degrees\n",angle2);
- }
-
- fprintf(fd,"\n-------------------------------------------------------------------------\n\n");
-
- /* Receiver */
-
- fprintf(fd,"Receiver site: %s\n",rcvr.name);
-
- if (rcvr.lat>=0.0)
- {
- fprintf(fd,"Site location: %.4f North / %.4f West",rcvr.lat, rcvr.lon);
- fprintf(fd, " (%s N / ", dec2dms(rcvr.lat));
- }
-
- else
- {
- fprintf(fd,"Site location: %.4f South / %.4f West",-rcvr.lat, rcvr.lon);
- fprintf(fd, " (%s S / ", dec2dms(rcvr.lat));
- }
-
- fprintf(fd, "%s W)\n", dec2dms(rcvr.lon));
-
- if (metric)
- {
- fprintf(fd,"Ground elevation: %.2f meters AMSL\n",METERS_PER_FOOT*GetElevation(rcvr));
- fprintf(fd,"Antenna height: %.2f meters AGL / %.2f meters AMSL\n",METERS_PER_FOOT*rcvr.alt, METERS_PER_FOOT*(rcvr.alt+GetElevation(rcvr)));
- }
-
- else
- {
- fprintf(fd,"Ground elevation: %.2f feet AMSL\n",GetElevation(rcvr));
- fprintf(fd,"Antenna height: %.2f feet AGL / %.2f feet AMSL\n",rcvr.alt, rcvr.alt+GetElevation(rcvr));
- }
- haavt=haat(rcvr);
+ fprintf(fd2,"Mode of propagation: %s\n",strmode);
- if (haavt>-4999.0)
- {
- if (metric)
- fprintf(fd,"Antenna height above average terrain: %.2f meters\n",METERS_PER_FOOT*haavt);
- else
- fprintf(fd,"Antenna height above average terrain: %.2f feet\n",haavt);
+ fprintf(fd2,"\n-------------------------------------------------------------------------\n\n");
}
- azimuth=Azimuth(rcvr,xmtr);
- angle1=ElevationAngle(rcvr,xmtr);
- angle2=ElevationAngle2(rcvr,xmtr,earthradius);
-
- if (metric)
- fprintf(fd,"Distance to %s: %.2f kilometers\n",xmtr.name,KM_PER_MILE*distance);
- else
- fprintf(fd,"Distance to %s: %.2f miles\n",xmtr.name,distance);
-
- fprintf(fd,"Azimuth to %s: %.2f degrees\n",xmtr.name,azimuth);
-
- if (angle1>=0.0)
- fprintf(fd,"Elevation to %s: %+.4f degrees\n",xmtr.name,angle1);
-
- else
- fprintf(fd,"Depression angle to %s: %+.4f degrees\n",xmtr.name,angle1);
-
- if (angle1!=angle2)
- {
- if (angle2<0.0)
- fprintf(fd,"Depression");
- else
- fprintf(fd,"Elevation");
+ fprintf(stdout,"\nPath Loss Report written to: \"%s\"\n",report_name);
+ fflush(stdout);
- fprintf(fd," angle to the first obstruction: %+.4f degrees\n",angle2);
+ ObstructionAnalysis(source, destination, LR.frq_mhz, fd2);
- }
+ fclose(fd2);
- fprintf(fd,"\n-------------------------------------------------------------------------\n\n");
+ /* Skip plotting the graph if ONLY a path-loss report is needed. */
- if (report=='y')
+ if (graph_it)
{
- /* Generate profile of the terrain. Create the path
- from transmitter to receiver because that's the
- way the original los() function did it, and going
- the other way can yield slightly different results. */
-
- ReadPath(xmtr,rcvr);
- h_r=GetElevation(rcvr)+rcvr.alt+earthradius;
- h_r_f1=h_r;
- h_r_fpt6=h_r;
- h_r_orig=h_r;
- h_t=GetElevation(xmtr)+xmtr.alt+earthradius;
- d_tx=5280.0*Distance(rcvr,xmtr);
- cos_tx_angle=((h_r*h_r)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r*d_tx);
- cos_tx_angle_f1=cos_tx_angle;
- cos_tx_angle_fpt6=cos_tx_angle;
-
- if (f)
- lambda=9.8425e8/(f*1e6);
-
- /* At each point along the path calculate the cosine
- of a sort of "inverse elevation angle" at the receiver.
- From the antenna, 0 deg. looks at the ground, and 90 deg.
- is parallel to the ground.
-
- Start at the receiver. If this is the lowest antenna,
- then terrain obstructions will be nearest to it. (Plus,
- that's the way the original los() did it.)
-
- Calculate cosines only. That's sufficient to compare
- angles and it saves the extra computational burden of
- acos(). However, note the inverted comparison: if
- acos(A) > acos(B), then B > A. */
-
- for (x=path.length-1; x>0; x--)
+ if (name[0]==0)
{
- site_x.lat=path.lat[x];
- site_x.lon=path.lon[x];
- site_x.alt=0.0;
-
- h_x=GetElevation(site_x)+earthradius;
- d_x=5280.0*Distance(rcvr,site_x);
-
- /* Deal with the LOS path first. */
-
- cos_test_angle=((h_r*h_r)+(d_x*d_x)-(h_x*h_x))/(2.0*h_r*d_x);
-
- if (cos_tx_angle>cos_test_angle)
- {
- if (h_r==h_r_orig)
- fprintf(fd,"SPLAT! detected obstructions at:\n\n");
-
- if (site_x.lat>=0.0)
- {
- if (metric)
- fprintf(fd,"\t%.4f N, %.4f W, %5.2f kilometers, %6.2f meters AMSL\n",site_x.lat, site_x.lon, KM_PER_MILE*(d_x/5280.0), METERS_PER_FOOT*(h_x-earthradius));
- else
- fprintf(fd,"\t%.4f N, %.4f W, %5.2f miles, %6.2f feet AMSL\n",site_x.lat, site_x.lon, d_x/5280.0, h_x-earthradius);
- }
-
- else
- {
- if (metric)
- fprintf(fd,"\t%.4f S, %.4f W, %5.2f kilometers, %6.2f meters AMSL\n",-site_x.lat, site_x.lon, KM_PER_MILE*(d_x/5280.0), METERS_PER_FOOT*(h_x-earthradius));
- else
-
- fprintf(fd,"\t%.4f S, %.4f W, %5.2f miles, %6.2f feet AMSL\n",-site_x.lat, site_x.lon, d_x/5280.0, h_x-earthradius);
- }
- }
+ /* Default filename and output file type */
- while (cos_tx_angle>cos_test_angle)
- {
- h_r+=1;
- cos_test_angle=((h_r*h_r)+(d_x*d_x)-(h_x*h_x))/(2.0*h_r*d_x);
- cos_tx_angle=((h_r*h_r)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r*d_tx);
- }
-
- if (f)
- {
- /* Now clear the first Fresnel zone, but don't
- clutter the obstruction report. */
-
- cos_tx_angle_f1=((h_r_f1*h_r_f1)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_f1*d_tx);
- h_los=sqrt(h_r_f1*h_r_f1+d_x*d_x-2*h_r_f1*d_x*cos_tx_angle_f1);
- h_f=h_los-sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
-
- while (h_f<h_x)
- {
- h_r_f1+=1;
- cos_tx_angle_f1=((h_r_f1*h_r_f1)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_f1*d_tx);
- h_los=sqrt(h_r_f1*h_r_f1+d_x*d_x-2*h_r_f1*d_x*cos_tx_angle_f1);
- h_f=h_los-sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
- }
-
- /* And clear the 60% F1 zone. */
-
- cos_tx_angle_fpt6=((h_r_fpt6*h_r_fpt6)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_fpt6*d_tx);
- h_los=sqrt(h_r_fpt6*h_r_fpt6+d_x*d_x-2*h_r_fpt6*d_x*cos_tx_angle_fpt6);
- h_f=h_los-0.6*sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
-
- while (h_f<h_x)
- {
- h_r_fpt6+=1;
- cos_tx_angle_fpt6=((h_r_fpt6*h_r_fpt6)+(d_tx*d_tx)-(h_t*h_t))/(2.0*h_r_fpt6*d_tx);
- h_los=sqrt(h_r_fpt6*h_r_fpt6+d_x*d_x-2*h_r_fpt6*d_x*cos_tx_angle_fpt6);
- h_f=h_los-0.6*sqrt(lambda*d_x*(d_tx-d_x)/d_tx);
- }
- }
- }
-
- if (h_r>h_r_orig)
- {
- if (metric)
- sprintf(string,"\nAntenna at %s must be raised to at least %.2f meters AGL\nto clear all obstructions detected by SPLAT!\n",rcvr.name, METERS_PER_FOOT*(h_r-GetElevation(rcvr)-earthradius));
- else
- sprintf(string,"\nAntenna at %s must be raised to at least %.2f feet AGL\nto clear all obstructions detected by SPLAT!\n",rcvr.name, h_r-GetElevation(rcvr)-earthradius);
+ strncpy(filename,"loss\0",5);
+ strncpy(term,"png\0",4);
+ strncpy(ext,"png\0",4);
}
else
- sprintf(string,"\nNo obstructions to LOS path due to terrain were detected by SPLAT!\n");
-
- if (f)
{
- if (h_r_fpt6>h_r_orig)
- {
- if (metric)
- sprintf(string_fpt6,"\nAntenna at %s must be raised to at least %.2f meters AGL\nto clear 60%c of the first Fresnel zone.\n",rcvr.name, METERS_PER_FOOT*(h_r_fpt6-GetElevation(rcvr)-earthradius),37);
-
- else
- sprintf(string_fpt6,"\nAntenna at %s must be raised to at least %.2f feet AGL\nto clear 60%c of the first Fresnel zone.\n",rcvr.name, h_r_fpt6-GetElevation(rcvr)-earthradius,37);
- }
+ /* Grab extension and terminal type from "name" */
- else
- sprintf(string_fpt6,"\n60%c of the first Fresnel zone is clear.\n",37);
-
- if (h_r_f1>h_r_orig)
- {
- if (metric)
- sprintf(string_f1,"\nAntenna at %s must be raised to at least %.2f meters AGL\nto clear the first Fresnel zone.\n",rcvr.name, METERS_PER_FOOT*(h_r_f1-GetElevation(rcvr)-earthradius));
+ for (x=0; name[x]!='.' && name[x]!=0 && x<254; x++)
+ filename[x]=name[x];
- else
- sprintf(string_f1,"\nAntenna at %s must be raised to at least %.2f feet AGL\nto clear the first Fresnel zone.\n",rcvr.name, h_r_f1-GetElevation(rcvr)-earthradius);
+ if (name[x]=='.')
+ {
+ for (y=0, z=x, x++; name[x]!=0 && x<254 && y<14; x++, y++)
+ {
+ term[y]=tolower(name[x]);
+ ext[y]=term[y];
+ }
+ ext[y]=0;
+ term[y]=0;
+ filename[z]=0;
}
else
- sprintf(string_f1,"\nThe first Fresnel zone is clear.\n\n");
+ { /* No extension -- Default is png */
+
+ filename[x]=0;
+ strncpy(term,"png\0",4);
+ strncpy(ext,"png\0",4);
+ }
}
- }
-
- fprintf(fd,"%s",string);
- if (f)
- {
- fprintf(fd,"%s",string_f1);
- fprintf(fd,"%s",string_fpt6);
- }
+ /* Either .ps or .postscript may be used
+ as an extension for postscript output. */
- fclose(fd);
+ if (strncmp(term,"postscript",10)==0)
+ strncpy(ext,"ps\0",3);
- /* Display report summary on terminal */
+ else if (strncmp(ext,"ps",2)==0)
+ strncpy(term,"postscript enhanced color\0",26);
- /* Line-of-sight status */
+ fd=fopen("splat.gp","w");
- fprintf(stdout,"%s",string);
+ fprintf(fd,"set grid\n");
+ fprintf(fd,"set yrange [%2.3f to %2.3f]\n", minloss, maxloss);
+ fprintf(fd,"set encoding iso_8859_1\n");
+ fprintf(fd,"set term %s\n",term);
+ fprintf(fd,"set title \"SPLAT! Loss Profile Along Path Between %s and %s (%.2f%c azimuth)\"\n",destination.name, source.name, Azimuth(destination,source),176);
- if (f)
- {
- /* Fresnel zone status */
+ if (metric)
+ fprintf(fd,"set xlabel \"Distance Between %s and %s (%.2f kilometers)\"\n",destination.name,source.name,KM_PER_MILE*Distance(destination,source));
+ else
+ fprintf(fd,"set xlabel \"Distance Between %s and %s (%.2f miles)\"\n",destination.name,source.name,Distance(destination,source));
- fprintf(stdout,"%s",string_f1);
- fprintf(stdout,"%s",string_fpt6);
- }
+ if (got_azimuth_pattern || got_elevation_pattern)
+ fprintf(fd,"set ylabel \"Total Path Loss (including TX antenna pattern) (dB)");
+ else
+ fprintf(fd,"set ylabel \"Longley-Rice Path Loss (dB)");
- fprintf(stdout, "\nObstruction report written to: \"%s\"\n",report_name);
+ fprintf(fd,"\"\nset output \"%s.%s\"\n",filename,ext);
+ fprintf(fd,"plot \"profile.gp\" title \"Path Loss\" with lines\n");
- fflush(stdout);
+ fclose(fd);
+
+ x=system("gnuplot splat.gp");
+
+ if (x!=-1)
+ {
+ unlink("splat.gp");
+ unlink("profile.gp");
+ unlink("reference.gp");
+
+ fprintf(stdout,"Path loss plot written to: \"%s.%s\"\n",filename,ext);
+ fflush(stdout);
+ }
+
+ else
+ fprintf(stderr,"\n*** ERROR: Error occurred invoking gnuplot!\n");
+ }
+
+ if (x!=-1)
+ unlink("profile.gp");
}
void SiteReport(struct site xmtr)
fprintf(fd,"\n---------------------------------------------------------------------------\n\n");
fclose(fd);
- fprintf(stdout,"\nSite analysis report written to: \"%s\"\n",report_name);
+ fprintf(stdout,"\nSite analysis report written to: \"%s\"",report_name);
}
void LoadTopoData(int max_lon, int min_lon, int max_lat, int min_lat)
/* This function loads the SDF files required
to cover the limits of the region specified. */
- int x, y, width, ymin, ymax;
+ int x, y, width, ymin, ymax;
width=ReduceAngle(max_lon-min_lon);
int LoadPLI(char *filename)
{
+ /* This function reads a SPLAT! path-loss output
+ file (-plo) for analysis and/or map generation. */
+
int error=0, max_west, min_west, max_north, min_north;
char string[80], *pointer=NULL;
double latitude=0.0, longitude=0.0, azimuth=0.0, elevation=0.0,
fprintf(stdout,"\nReading \"%s\"... ",filename);
fflush(stdout);
- fscanf(fd,"%lf, %lf, %lf, %lf, %lf",&latitude, &longitude, &azimuth, &elevation, &loss);
+ fgets(string,78,fd);
+ sscanf(string,"%lf, %lf, %lf, %lf, %lf",&latitude, &longitude, &azimuth, &elevation, &loss);
while (feof(fd)==0)
{
- if (loss>225.0)
- loss=225.0;
-
- if (loss<75.0)
- loss=75.0;
-
- loss-=75.0;
- loss/=10.0;
- loss+=1.0;
+ if (loss>255.0)
+ loss=255.0;
if (loss<=(double)maxdB)
- OrMask(latitude,longitude,((unsigned char)(loss))<<3);
+ PutSignal(latitude,longitude,((unsigned char)round(loss)));
- fscanf(fd,"%lf, %lf, %lf, %lf, %lf",&latitude, &longitude, &azimuth, &elevation, &loss);
+ fgets(string,78,fd);
+ sscanf(string,"%lf, %lf, %lf, %lf, %lf",&latitude, &longitude, &azimuth, &elevation, &loss);
}
fclose(fd);
fclose(fd);
- fprintf(stdout, "KML file written to: \"%s\"\n",report_name);
+ fprintf(stdout, "\nKML file written to: \"%s\"",report_name);
fflush(stdout);
}
rxlat, rxlon, txlat, txlon, west_min, west_max,
north_min, north_max;
- unsigned char coverage=0, LRmap=0, ext[20], terrain_plot=0,
- elevation_plot=0, height_plot=0,
+ unsigned char coverage=0, LRmap=0, terrain_plot=0,
+ elevation_plot=0, height_plot=0, map=0, nf=0,
longley_plot=0, cities=0, bfs=0, txsites=0,
- count, report='y', norm=0, topomap=0, geo=0,
- kml=0;
+ norm=0, topomap=0, geo=0, kml=0, pt2pt_mode=0,
+ area_mode=0, max_txsites, ngs=0, nolospath=0,
+ nositereports=0;
char mapfile[255], header[80], city_file[5][255],
elevation_file[255], height_file[255],
longley_file[255], terrain_file[255],
string[255], rxfile[255], *env=NULL,
- txfile[255], map=0, boundary_file[5][255],
+ txfile[255], boundary_file[5][255],
udt_file[255], rxsite=0, plo_filename[255],
- pli_filename[255], nf=0;
+ pli_filename[255], ext[20];
double altitude=0.0, altitudeLR=0.0, tx_range=0.0,
rx_range=0.0, deg_range=0.0, deg_limit,
deg_range_lon, er_mult, freq=0.0;
- struct site tx_site[4], rx_site;
+ struct site tx_site[32], rx_site;
FILE *fd;
if (argc==1)
{
fprintf(stdout,"\n\t\t --==[ SPLAT! v%s Available Options... ]==--\n\n",splat_version);
- fprintf(stdout," -t txsite(s).qth (max of 4)\n");
+ fprintf(stdout," -t txsite(s).qth (max of 4 with -c, max of 30 with -L)\n");
fprintf(stdout," -r rxsite.qth\n");
fprintf(stdout," -c plot coverage of TX(s) with an RX antenna at X feet/meters AGL\n");
fprintf(stdout," -L plot path loss map of TX based on an RX at X feet/meters AGL\n");
fprintf(stdout," -s filename(s) of city/site file(s) to import (5 max)\n");
- fprintf(stdout," -b filename(s) of cartographic boundary file(s) to import (max of 5)\n");
+ fprintf(stdout," -b filename(s) of cartographic boundary file(s) to import (5 max)\n");
fprintf(stdout," -p filename of terrain profile graph to plot\n");
fprintf(stdout," -e filename of terrain elevation graph to plot\n");
fprintf(stdout," -h filename of terrain height graph to plot\n");
fprintf(stdout," -o filename of topographic map to generate (.ppm)\n");
fprintf(stdout," -u filename of user-defined terrain file to import\n");
fprintf(stdout," -d sdf file directory path (overrides path in ~/.splat_path file)\n");
- fprintf(stdout," -n no analysis, brief report\n");
- fprintf(stdout," -N no analysis, no report\n");
fprintf(stdout," -m earth radius multiplier\n");
+ fprintf(stdout," -n do not plot LOS paths in .ppm maps\n");
+ fprintf(stdout," -N do not produce unnecessary site or obstruction reports\n");
fprintf(stdout," -f frequency for Fresnel zone calculation (MHz)\n");
fprintf(stdout," -R modify default range for -c or -L (miles/kilometers)\n");
fprintf(stdout," -db maximum loss contour to display on path loss maps (80-230 dB)\n");
fprintf(stdout," -nf do not plot Fresnel zones in height plots\n");
- fprintf(stdout," -plo filename of path-loss output file\n");
+ fprintf(stdout," -fz Fresnel zone clearance percentage (default = 60)\n");
+ fprintf(stdout," -ngs display greyscale topography as white in .ppm files\n");
+ fprintf(stdout," -erp override ERP in .lrp file (Watts)\n");
fprintf(stdout," -pli filename of path-loss input file\n");
+ fprintf(stdout," -plo filename of path-loss output file\n");
fprintf(stdout," -udt filename of user defined terrain input file\n");
+ fprintf(stdout," -kml generate Google Earth (.kml) compatible output\n");
fprintf(stdout," -geo generate an Xastir .geo georeference file (with .ppm output)\n");
- fprintf(stdout," -kml generate a Google Earth .kml file (for point-to-point links)\n");
fprintf(stdout," -metric employ metric rather than imperial units for all user I/O\n\n");
-
fprintf(stdout,"If that flew by too fast, consider piping the output through 'less':\n");
fprintf(stdout,"\n\tsplat | less\n\n");
fprintf(stdout,"Type 'man splat', or see the documentation for more details.\n\n");
fflush(stdout);
+
return 1;
}
y=argc-1;
+ kml=0;
+ geo=0;
metric=0;
rxfile[0]=0;
txfile[0]=0;
sdf_path[0]=0;
udt_file[0]=0;
path.length=0;
+ max_txsites=30;
LR.frq_mhz=0.0;
+ fzone_clearance=0.6;
rx_site.lat=91.0;
rx_site.lon=361.0;
+ longley_file[0]=0;
plo_filename[0]=0;
pli_filename[0]=0;
earthradius=EARTHRADIUS;
tx_site[x].lon=361.0;
}
- for (x=0; x<MAXSLOTS; x++)
+ for (x=0; x<MAXPAGES; x++)
{
dem[x].min_el=32768;
dem[x].max_el=-32768;
}
}
+ if (strcmp(argv[x],"-fz")==0)
+ {
+ z=x+1;
+
+ if (z<=y && argv[z][0] && argv[z][0]!='-')
+ {
+ sscanf(argv[z],"%lf",&fzone_clearance);
+
+ if (fzone_clearance<0.0 || fzone_clearance>100.0)
+ fzone_clearance=60.0;
+
+ fzone_clearance/=100.0;
+ }
+ }
+
if (strcmp(argv[x],"-o")==0)
{
z=x+1;
if (z<=y && argv[z][0] && argv[z][0]!='-')
{
sscanf(argv[z],"%lf",&altitude);
+ map=1;
coverage=1;
+ area_mode=1;
+ max_txsites=4;
}
}
{
strncpy(terrain_file,argv[z],253);
terrain_plot=1;
+ pt2pt_mode=1;
}
}
{
strncpy(elevation_file,argv[z],253);
elevation_plot=1;
+ pt2pt_mode=1;
}
}
{
strncpy(height_file,argv[z],253);
height_plot=1;
+ pt2pt_mode=1;
}
if (strcmp(argv[x],"-H")==0)
norm=0;
}
- if (strcmp(argv[x],"-n")==0)
- {
- report='n';
- map=1;
- }
-
- if (strcmp(argv[x],"-N")==0)
- {
- report='N';
- map=1;
- }
-
if (strcmp(argv[x],"-metric")==0)
metric=1;
if (strcmp(argv[x],"-nf")==0)
nf=1;
+ if (strcmp(argv[x],"-ngs")==0)
+ ngs=1;
+
+ if (strcmp(argv[x],"-n")==0)
+ nolospath=1;
+
+ if (strcmp(argv[x],"-N")==0)
+ {
+ nolospath=1;
+ nositereports=1;
+ }
+
if (strcmp(argv[x],"-d")==0)
{
z=x+1;
z=x+1;
- while (z<=y && argv[z][0] && argv[z][0]!='-' && txsites<4)
+ while (z<=y && argv[z][0] && argv[z][0]!='-' && txsites<30)
{
strncpy(txfile,argv[z],253);
tx_site[txsites]=LoadQTH(txfile);
if (z<=y && argv[z][0] && argv[z][0]!='-')
{
sscanf(argv[z],"%lf",&altitudeLR);
+ map=1;
+ LRmap=1;
+ area_mode=1;
if (coverage)
fprintf(stdout,"c and L are exclusive options, ignoring L.\n");
- else
- {
- LRmap=1;
- ReadLRParm(txfile);
- }
}
}
{
strncpy(longley_file,argv[z],253);
longley_plot=1;
- /* Doing this twice is harmless */
- ReadLRParm(txfile);
+ pt2pt_mode=1;
}
}
strncpy(rxfile,argv[z],253);
rx_site=LoadQTH(rxfile);
rxsite=1;
+ pt2pt_mode=1;
}
}
}
}
+ if (strcmp(argv[x],"-erp")==0)
+ {
+ z=x+1;
+
+ if (z<=y && argv[z][0] && argv[z][0]!='-')
+ {
+ sscanf(argv[z],"%lf",&forced_erp);
+
+ if (forced_erp<0.0)
+ forced_erp=-1.0;
+ }
+ }
+
if (strcmp(argv[x],"-plo")==0)
{
z=x+1;
map=0;
topomap=1;
- report='N';
}
else
{
y=LoadPLI(pli_filename);
- for (x=0; x<txsites; x++)
+ for (x=0; x<txsites && x<max_txsites; x++)
PlaceMarker(tx_site[x]);
if (rxsite)
LoadCities(city_file[x]);
}
- WritePPMLR(mapfile,geo);
+ WritePPMLR(mapfile,geo,kml,ngs,tx_site,txsites);
exit(0);
}
min_lon=(int)floor(tx_site[0].lon);
max_lon=(int)floor(tx_site[0].lon);
- for (y=0, z=0; z<txsites; z++)
+ for (y=0, z=0; z<txsites && z<max_txsites; z++)
{
txlat=(int)floor(tx_site[z].lat);
txlon=(int)floor(tx_site[z].lon);
max_lon=rxlon;
}
+
/* Load the required SDF files */
LoadTopoData(max_lon, min_lon, max_lat, min_lat);
- if (coverage | LRmap | topomap)
+ if (area_mode || topomap)
{
- if (LRmap)
- txsites=1;
-
- for (z=0; z<txsites; z++)
+ for (z=0; z<txsites && z<max_txsites; z++)
{
/* "Ball park" estimates used to load any additional
SDF files required to conduct this analysis. */
deg_range=max_range/69.0;
/* Prevent the demand for a really wide coverage
- from allocating more slots than are available
+ from allocating more "pages" than are available
in memory. */
- switch (MAXSLOTS)
+ switch (MAXPAGES)
{
case 2: deg_limit=0.25;
break;
if (deg_range_lon>deg_limit)
deg_range_lon=deg_limit;
-
north_min=(int)floor(tx_site[z].lat-deg_range);
north_max=(int)floor(tx_site[z].lat+deg_range);
LoadTopoData(max_lon, min_lon, max_lat, min_lat);
}
-
+
if (udt_file[0])
LoadUDT(udt_file);
- if (mapfile[0] && topomap==0)
- map=1;
-
- if (freq==0.0 && nf==0)
- freq=LR.frq_mhz;
-
- else if (nf==1)
- freq=0.0;
-
- if (coverage | LRmap)
- {
- for (x=0; x<txsites; x++)
- {
- if (coverage)
- PlotCoverage(tx_site[x],altitude);
-
- if (LRmap)
- PlotLRMap(tx_site[x],altitudeLR,plo_filename);
-
- PlaceMarker(tx_site[x]);
-
- if (report!='N')
- SiteReport(tx_site[x]);
- }
- map=1;
- }
+ /***** Let the SPLATting begin! *****/
- if (coverage==0 && LRmap==0)
+ if (pt2pt_mode)
{
PlaceMarker(rx_site);
- for (x=0; x<txsites; x++)
+ if (longley_plot)
{
- PlaceMarker(tx_site[x]);
+ /* Grab extension to determine graphic file type */
- if (report=='y')
- {
- switch (x)
- {
- case 0:
- PlotPath(tx_site[x],rx_site,1);
- break;
-
- case 1:
- PlotPath(tx_site[x],rx_site,8);
- break;
+ for (x=0; longley_file[x]!='.' && longley_file[x]!=0 && x<80; x++);
- case 2:
- PlotPath(tx_site[x],rx_site,16);
- break;
+ if (longley_file[x]=='.')
+ {
+ ext[0]='.';
+ for (y=1, z=x, x++; longley_file[x]!=0 && x<253 && y<14; x++, y++)
+ ext[y]=longley_file[x];
- case 3:
- PlotPath(tx_site[x],rx_site,32);
- }
+ ext[y]=0;
+ longley_file[z]=0;
}
- if (report!='N')
- ObstructionReport(tx_site[x],rx_site,report,freq);
-
- if (kml)
- WriteKML(tx_site[x],rx_site);
- }
- }
-
- if (map | topomap)
- {
- if (bfs)
- {
- for (x=0; x<bfs; x++)
- LoadBoundaries(boundary_file[x]);
- }
-
- if (cities)
- {
- for (x=0; x<cities; x++)
- LoadCities(city_file[x]);
+ else
+ {
+ ext[0]=0; /* No extension */
+ longley_file[x]=0;
+ }
}
- if (LRmap)
- WritePPMLR(mapfile,geo);
- else
- WritePPM(mapfile,geo);
- }
-
- if (terrain_plot)
- {
- if (txsites>1)
+ if (terrain_plot)
{
for (x=0; terrain_file[x]!='.' && terrain_file[x]!=0 && x<80; x++);
ext[0]=0; /* No extension */
terrain_file[x]=0;
}
-
- for (count=0; count<txsites; count++)
- {
- sprintf(string,"%s-%c%s%c",terrain_file,'1'+count,ext,0);
- GraphTerrain(tx_site[count],rx_site,string);
- }
}
- else
- GraphTerrain(tx_site[0],rx_site,terrain_file);
- }
-
- if (elevation_plot)
- {
- if (txsites>1)
+ if (elevation_plot)
{
for (x=0; elevation_file[x]!='.' && elevation_file[x]!=0 && x<80; x++);
ext[0]=0; /* No extension */
elevation_file[x]=0;
}
-
- for (count=0; count<txsites; count++)
- {
- sprintf(string,"%s-%c%s%c",elevation_file,'1'+count,ext,0);
- GraphElevation(tx_site[count],rx_site,string);
- }
}
- else
- GraphElevation(tx_site[0],rx_site,elevation_file);
- }
-
- if (height_plot)
- {
- if (txsites>1)
+ if (height_plot)
{
for (x=0; height_file[x]!='.' && height_file[x]!=0 && x<80; x++);
ext[0]=0; /* No extension */
height_file[x]=0;
}
-
- for (count=0; count<txsites; count++)
- {
- sprintf(string,"%s-%c%s%c",height_file,'1'+count,ext,0);
- GraphHeight(tx_site[count],rx_site,string,freq,norm);
- }
}
- else
- GraphHeight(tx_site[0],rx_site,height_file,freq,norm);
- }
-
- if (longley_plot)
- {
- if (txsites>1)
+ for (x=0; x<txsites && x<4; x++)
{
- for (x=0; longley_file[x]!='.' && longley_file[x]!=0 && x<80; x++);
+ PlaceMarker(tx_site[x]);
- if (longley_file[x]=='.') /* extension */
+ if (nolospath==0)
{
- ext[0]='.';
- for (y=1, z=x, x++; longley_file[x]!=0 && x<253 && y<14; x++, y++)
- ext[y]=longley_file[x];
+ switch (x)
+ {
+ case 0:
+ PlotPath(tx_site[x],rx_site,1);
+ break;
- ext[y]=0;
- longley_file[z]=0;
+ case 1:
+ PlotPath(tx_site[x],rx_site,8);
+ break;
+
+ case 2:
+ PlotPath(tx_site[x],rx_site,16);
+ break;
+
+ case 3:
+ PlotPath(tx_site[x],rx_site,32);
+ }
}
+ if (nositereports==0)
+ SiteReport(tx_site[x]);
+
+ if (kml)
+ WriteKML(tx_site[x],rx_site);
+
+ if (txsites>1)
+ sprintf(string,"%s-%c%s%c",longley_file,'1'+x,ext,0);
else
+ sprintf(string,"%s%s%c",longley_file,ext,0);
+
+ if (nositereports==0)
{
- ext[0]=0; /* No extension */
- longley_file[x]=0;
+ if (longley_file[0]==0)
+ {
+ ReadLRParm(tx_site[x],0);
+ PathReport(tx_site[x],rx_site,string,0); }
+
+ else
+ {
+ ReadLRParm(tx_site[x],1);
+ PathReport(tx_site[x],rx_site,string,longley_file[0]);
+ }
+ }
+
+ if (terrain_plot)
+ {
+ if (txsites>1)
+ sprintf(string,"%s-%c%s%c",terrain_file,'1'+x,ext,0);
+ else
+ sprintf(string,"%s%s%c",terrain_file,ext,0);
+
+ GraphTerrain(tx_site[x],rx_site,string);
+ }
+
+ if (elevation_plot)
+ {
+ if (txsites>1)
+ sprintf(string,"%s-%c%s%c",elevation_file,'1'+x,ext,0);
+ else
+ sprintf(string,"%s%s%c",elevation_file,ext,0);
+
+ GraphElevation(tx_site[x],rx_site,string);
+ }
+
+ if (height_plot)
+ {
+ if (freq==0.0 && nf==0)
+ freq=LR.frq_mhz;
+
+ if (txsites>1)
+ sprintf(string,"%s-%c%s%c",height_file,'1'+x,ext,0);
+ else
+ sprintf(string,"%s%s%c",height_file,ext,0);
+
+ GraphHeight(tx_site[x],rx_site,string,freq,norm);
}
+ }
+ }
+
+ if (area_mode && topomap==0)
+ {
+ for (x=0; x<txsites && x<max_txsites; x++)
+ {
+ if (coverage)
+ PlotCoverage(tx_site[x],altitude);
+
+ else if (ReadLRParm(tx_site[x],1))
+ PlotLRMap(tx_site[x],altitudeLR,plo_filename);
+
+ SiteReport(tx_site[x]);
+ }
+ }
+
+ if (map || topomap)
+ {
+ /* Label the map */
- for (count=0; count<txsites; count++)
+ if (cities)
+ {
+ if (kml==0)
{
- sprintf(string,"%s-%c%s%c",longley_file,'1'+count,ext,0);
- GraphLongley(tx_site[count],rx_site,string);
+ for (x=0; x<txsites && x<max_txsites; x++)
+ PlaceMarker(tx_site[x]);
}
+
+ for (y=0; y<cities; y++)
+ LoadCities(city_file[y]);
+ }
+
+ /* Load city and county boundary data files */
+
+ if (bfs)
+ {
+ for (y=0; y<bfs; y++)
+ LoadBoundaries(boundary_file[y]);
}
+ /* Plot the map */
+
+ if (coverage || pt2pt_mode || topomap)
+ WritePPM(mapfile,geo,kml,ngs);
+
else
- GraphLongley(tx_site[0],rx_site,longley_file);
+ {
+ if (LR.erp==0.0)
+ WritePPMLR(mapfile,geo,kml,ngs,tx_site,txsites);
+ else
+ WritePPMSS(mapfile,geo,kml,ngs,tx_site,txsites);
+ }
}
+ printf("\n");
+
+ /* That's all, folks! */
+
return 0;
}
-
projects / debian / splat / blobdiff