[R] Inserting color into an irregular grid comprised of polygons
Jim Lemon
jim at bitwrit.com.au
Thu Jan 16 21:34:00 CET 2014
On 01/17/2014 05:12 AM, Morway, Eric wrote:
> As a follow up to this thread started nearly a month ago, I'm in need of
> help sorting out the R code that will create a log-scale legend in the call
> to "image.plot" below. As the last line of the code provided below shows,
> I attempted to force the labeling through the argument "legend.lab", to no
> avail. The -4 on the legend label should actually be 0.0001, -3 = 0.001,
> etc etc according to the inverse of the log base 10.
>
> If possible, it would be nice to create legend similar to the one created
> here<http://r.789695.n4.nabble.com/sppolot-fill-below-minimum-legend-value-td902841.html>,
> only that uses spplot and I need to stick with the base graphics since I'm
> using irregular polygons to draw the figure. However, if someone is able
> to show how to replace the -4 with 0.0001 etc, that would be a great place
> to start. Here's the reproducible R code (note that log10 is taken of the
> matrix "vals" and is what guides the color fill in the plot):
>
> library(gsubfn) #uses paste0 func
> library(colorRamps) #uses blue2green2red()
> library(fields) #uses image.plot(..., legend.only=TRUE, ...)
>
>
> z.space<-
> c(0.6790521,0.3454826,0.1872356,0.0891079,0.1525315,0.1088516,0.0950484,0.1128700,0.1247511,0.1188105,0.1143682,0.1232529,0.0930168,0.0751814,0.0511553,0.0244765,0.0424162,0.0435835,0.0577441,0.0471291,0.0974984,0.0303579,0.0234230,0.0378371,0.0396388,0.0278040,0.0427108,0.0450803,0.0735903,0.1499654,0.0235646,0.0309285,0.0770295,0.0687763,0.1007385,0.0666026,0.1083643,0.1092819,0.1372624,0.2248670,0.2620903,0.4606435,0.6262846,1.7111480,1.7111480,1.7111480,1.7111480,1.7111480,1.7111480,1.7111480,0.8662780,1.1220410,0.5368302)
> x.space<-
> c(0.4477580,0.4058683,0.5047908,0.3488354,0.3170296,0.2280360,0.2371574,0.1658813,0.2098874,0.2441864,0.3050745,0.4087275,0.4448988,0.4416195,0.4020654,0.0862620,0.0332546,0.0871109,0.3531576,0.3037825,0.2396926,0.2351304,0.2144404,0.0733572,0.0338528,0.2016122,0.1533454,0.1265044,0.0932833,0.0481462,0.0662010,0.0150457,0.0481462,0.0270822,0.0318521,0.0995603,0.0583223,0.0371142,0.0854215,0.0577332,0.0883671,0.0786467,0.0787786,0.1135672,0.0897309,0.1659446,0.8536263,0.8536263,0.8536263,0.8536263,0.5794602,0.2741660)
>
> #Obs<- read.csv("Obs_Loc_for_R.txt",header=T)
>
> x.range<- c(0,sum(x.space))
> z.range<- c(0,sum(z.space))
> z.sum<- sum(z.space)
>
> ###################
> # Read HK
> ###################
>
> vals<-
> matrix(c(0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,
>
> 0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.1498912E-03,0.5414670E-02,
>
> 0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.9835267E-03,0.6889354E-01,0.1004814E-01,
>
> 0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.6619056E-03,0.4637940E-02,0.6664155E-01,0.1453067E-01,
>
> 0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.6710035E-03,0.1815741E-01,0.9034086E-01,0.2995279E-01,
>
> 0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.5446172E-02,0.3912347E-01,0.6302311E-01,0.3142162E-01,
>
> 0.1624253E-03,0.8807862E-04,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.2207096E-02,0.2737333E-01,0.5499471E-01,0.4478347E-01,0.3109830E-01,
>
> 0.1459764E-03,0.1272388E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.3896962E-02,0.4446387E-01,0.6582496E-01,0.4850795E-01,0.6290437E-01,
>
> 0.1106772E-03,0.2563490E-03,0.2326099E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.1038738E-01,0.4728762E-01,0.6748933E-01,0.5443041E-01,0.7390492E-01,
>
> 0.1491018E-03,0.4064940E-03,0.1790473E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.3760740E-03,0.3576420E-01,0.3302301E-01,0.6788430E-01,0.6265076E-01,0.7570447E-01,
>
> 0.2646171E-03,0.6089240E-03,0.1563286E-03,0.2083332E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.2977450E-04,0.1590487E-03,0.1240888E+00,0.2019374E-01,0.6820355E-01,0.7159009E-01,0.7410751E-01,
>
> 0.1238532E-02,0.2645979E-02,0.2909128E-03,0.3345401E-03,0.3404909E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.4212634E-04,0.3471802E-04,0.3374332E-04,0.2903843E-04,0.1301713E-03,0.2447790E-01,0.2479411E-01,0.6293514E-01,0.7443890E-01,0.7670832E-01,
>
> 0.1309739E-01,0.2563637E-01,0.7939482E-03,0.4279433E-03,0.4211852E-03,0.4679168E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.4513439E-04,0.5510390E-04,0.5106100E-04,0.3568510E-04,0.1932924E-04,0.2309740E-04,0.3484642E-04,0.2303861E-03,0.7116658E-03,0.4851235E-01,0.5422235E-01,0.7393300E-01,0.5992470E-01,
>
> 0.5216245E-01,0.7901551E-01,0.1998644E-02,0.5713268E-03,0.4577809E-03,0.3623582E-03,0.4273419E-03,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.3327895E-04,0.1661508E-04,0.4661848E-04,0.9217941E-04,0.9985013E-04,0.1027270E-03,0.5621401E-04,0.1546029E-04,0.2530272E-04,0.5851532E-04,0.4824676E-03,0.1373458E-03,0.7548517E-01,0.5205899E-01,0.7774446E-01,0.6158688E-01,
>
> 0.5384069E-01,0.9219270E-01,0.4106641E-02,0.7405762E-03,0.4043071E-03,0.3408503E-03,0.4168981E-03,0.1437158E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.1232444E-02,0.4622970E-03,0.2436128E-03,0.1953873E-03,0.2284356E-03,0.2535794E-03,0.2194460E-03,0.1364029E-03,0.9066684E-04,0.3719590E-04,0.4699193E-04,0.1076661E-03,0.1080930E-02,0.3043486E-03,0.9257647E-01,0.5773572E-01,0.1080207E+00,0.6257250E-01,
>
> 0.5281573E-01,0.9269738E-01,0.6195772E-02,0.8786800E-03,0.3837434E-03,0.3435221E-03,0.4230360E-03,0.1425388E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.2614737E-02,0.2557613E-02,0.2429492E-02,0.2279371E-02,0.1199601E-02,0.9997486E-03,0.8834830E-03,0.7776493E-03,0.5830090E-03,0.3870902E-03,0.2149853E-03,0.1332767E-03,0.7279420E-04,0.7398164E-04,0.1629017E-03,0.1867165E-02,0.6806691E-03,0.1005101E+00,0.6195525E-01,0.1334367E+00,0.7303004E-01,
>
> 0.5949255E-01,0.8381032E-01,0.8720866E-02,0.1105375E-02,0.3968406E-03,0.3551831E-03,0.4371457E-03,0.1429820E-02,0.2397559E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.7961092E-02,0.3563717E-02,0.3750839E-02,0.3960508E-02,0.4107561E-02,0.2951662E-02,0.3639369E-02,0.3454647E-02,0.2436310E-02,0.1279606E-02,0.6589286E-03,0.3147616E-03,0.1889157E-03,0.1128230E-03,0.1231117E-03,0.2348349E-03,0.3053442E-02,0.1429933E-02,0.1094067E+00,0.6660549E-01,0.1604875E+00,0.7308186E-01,
>
> 0.6010304E-01,0.8025564E-01,0.1286649E-01,0.1385014E-02,0.4028917E-03,0.3922943E-03,0.4770908E-03,0.1454291E-02,0.2336916E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.1647675E-01,0.1011602E-01,0.1095920E-01,0.6187991E-02,0.7484856E-02,0.8808350E-02,0.1342270E-01,0.1996092E-01,0.2151685E-01,0.1059982E-01,0.3474249E-02,0.1293193E-02,0.5032435E-03,0.2808538E-03,0.1582102E-03,0.1883053E-03,0.3602254E-03,0.5583650E-02,0.3567805E-02,0.1237909E+00,0.7336210E-01,0.2113802E+00,0.7904304E-01,
>
> 0.6252230E-01,0.6468046E-01,0.2035429E-01,0.1870555E-02,0.4749717E-03,0.5111303E-03,0.6577721E-03,0.1507024E-02,0.2246557E-02,0.2849760E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.6675010E-02,0.1722321E-01,0.2048179E-01,0.1607310E-01,0.1081760E-01,0.1509175E-01,0.2016206E-01,0.4317668E-01,0.1272580E+00,0.2093490E+00,0.4788971E-01,0.9771262E-02,0.2621271E-02,0.7990651E-03,0.3928475E-03,0.2674300E-03,0.2513305E-03,0.5292210E-03,0.1045705E-01,0.7563207E-02,0.1316331E+00,0.8177958E-01,0.2974619E+00,0.7785880E-01,
>
> 0.7063367E-01,0.5191958E-01,0.2540953E-01,0.3055553E-02,0.7248933E-03,0.8903009E-03,0.1012931E-02,0.1588296E-02,0.1805916E-02,0.3083439E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.1098240E-01,0.1308467E-01,0.1692970E-01,0.2387417E-01,0.3202713E-01,0.1601916E-01,0.2529038E-01,0.4072065E-01,0.1058370E+00,0.2987763E+00,0.3754694E+00,0.1221606E+00,0.2176033E-01,0.4616003E-02,0.1101954E-02,0.4218264E-03,0.2537318E-03,0.2518645E-03,0.7206248E-03,0.1972273E-01,0.9191978E-02,0.1231698E+00,0.8643236E-01,0.2099450E+00,0.7579815E-01,
>
> 0.9006456E-01,0.4065609E-01,0.3351370E-01,0.5199741E-02,0.1720484E-02,0.2490027E-02,0.2412121E-02,0.1822400E-02,0.1183448E-02,0.2931238E-02,0.4199585E-02,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.0000000E+00,0.5264653E-01,0.7212964E-02,0.8577601E-02,0.1486374E-01,0.2690323E-01,0.2518052E-01,0.3804494E-01,0.6426261E-01,0.7835973E-01,0.2499310E+00,0.6706516E+00,0.5886371E+00,0.2415061E+00,0.4384018E-01,0.7361790E-02,0.1282174E-02,0.3500221E-03,0.1682523E-03,0.1549343E-03,0.7264481E-03,0.2955076E-01,0.5826306E-02,0.1042731E+00,0.9299631E-01,0.9913468E-01,0.7644407E-01,
>
> 0.1364559E+00,0.3198620E-01,0.3564783E-01,0.8317512E-02,0.3968001E-02,0.6579916E-02,0.6097231E-02,0.2168262E-02,0.7698590E-03,0.3656941E-02,0.4060849E-02,0.2150453E-02,0.4767156E-02,0.2768307E-01,0.1354281E+00,0.1058603E+00,0.1143571E+00,0.1094348E+00,0.5325599E-02,0.2128505E-03,0.1202116E-02,0.2333544E-02,0.4377179E-01,0.3853912E+00,0.7224898E+00,0.1993333E+00,0.4044349E-01,0.6103219E-02,0.5073220E-02,0.1318864E-01,0.1894912E-01,0.3331478E-01,0.5236714E-01,0.9241241E-01,0.1486214E+00,0.3895471E+00,0.1581072E+01,0.1518637E+01,0.4936490E+00,0.6101139E-01,0.8288936E-02,0.1151226E-02,0.2353228E-03,0.6690447E-04,0.8802950E-04,0.5986074E-03,0.3520550E-01,0.3772795E-02,0.9010564E-01,0.9541733E-01,0.5234505E-01,0.6758159E-01,
>
> 0.1561256E+00,0.2962665E-01,0.3461261E-01,0.9982125E-02,0.5716534E-02,0.1020576E-01,0.8432660E-02,0.2423745E-02,0.6104935E-03,0.3819768E-02,0.4032974E-02,0.2105578E-02,0.4989981E-02,0.3189802E-01,0.1498627E+00,0.1180727E+00,0.1422254E+00,0.1287981E+00,0.5247850E-02,0.1696853E-03,0.1279081E-02,0.2291563E-02,0.4118333E-01,0.4988456E+00,0.1068528E+01,0.1951103E+00,0.3736538E-01,0.4848815E-02,0.4966058E-02,0.1334049E-01,0.2114189E-01,0.3716179E-01,0.5811247E-01,0.1018696E+00,0.1635320E+00,0.4505116E+00,0.2271184E+01,0.2674734E+01,0.5933481E+00,0.6399865E-01,0.8009510E-02,0.1046734E-02,0.2136105E-03,0.3783259E-04,0.6570102E-04,0.5359800E-03,0.3639964E-01,0.3200932E-02,0.8560315E-01,0.9563935E-01,0.3937848E-01,0.6603998E-01,
>
> 0.1770920E+00,0.2815360E-01,0.3558350E-01,0.1314717E-01,0.7972147E-02,0.1458678E-01,0.1163480E-01,0.2854676E-02,0.4650055E-03,0.4173236E-02,0.4049313E-02,0.2245366E-02,0.5771731E-02,0.3960153E-01,0.1692946E+00,0.1246485E+00,0.1705126E+00,0.1496826E+00,0.5232939E-02,0.1510017E-03,0.1469552E-02,0.2433426E-02,0.3593167E-01,0.5412888E+00,0.1745517E+01,0.1674432E+00,0.3545550E-01,0.4399809E-02,0.4365363E-02,0.1490021E-01,0.2468639E-01,0.4213111E-01,0.6413683E-01,0.1088669E+00,0.1709121E+00,0.4688794E+00,0.3171674E+01,0.5556163E+01,0.6425149E+00,0.5869946E-01,0.7152567E-02,0.1006172E-02,0.1717671E-03,0.2206736E-04,0.4950469E-04,0.5135077E-03,0.3835706E-01,0.2940349E-02,0.8324928E-01,0.9540258E-01,0.2839582E-01,0.6439322E-01,
>
> 0.1617234E+00,0.2846193E-01,0.3669326E-01,0.1763696E-01,0.1130720E-01,0.2036462E-01,0.1592378E-01,0.4274657E-02,0.7020108E-03,0.4963929E-02,0.4933829E-02,0.2861772E-02,0.8078864E-02,0.5661243E-01,0.1991557E+00,0.1146786E+00,0.1467128E+00,0.1219753E+00,0.5339475E-02,0.1856591E-03,0.1964900E-02,0.3162345E-02,0.2701869E-01,0.2862444E+00,0.7358219E+00,0.1079985E+00,0.3152878E-01,0.6321066E-02,0.6576352E-02,0.2238160E-01,0.3182091E-01,0.4969640E-01,0.7057704E-01,0.1102752E+00,0.1622570E+00,0.4555364E+00,0.2129294E+01,0.2940004E+01,0.4349839E+00,0.4321371E-01,0.6097398E-02,0.8030692E-03,0.1375305E-03,0.1775944E-04,0.4023215E-04,0.5216482E-03,0.4872911E-01,0.3371748E-02,0.8572043E-01,0.9949564E-01,0.2557490E-01,0.6253364E-01,
>
> 0.1260119E+00,0.3063453E-01,0.4038322E-01,0.2054567E-01,0.1437562E-01,0.2457226E-01,0.1929946E-01,0.5728094E-02,0.1552052E-02,0.6059005E-02,0.6596521E-02,0.4066243E-02,0.1144022E-01,0.9162235E-01,0.2302253E+00,0.9336674E-01,0.9732815E-01,0.7804549E-01,0.5186473E-02,0.2981081E-03,0.2791004E-02,0.4666605E-02,0.1915244E-01,0.9123682E-01,0.1511566E+00,0.5955763E-01,0.3277183E-01,0.1322382E-01,0.1484807E-01,0.3359206E-01,0.4177701E-01,0.5769788E-01,0.7475993E-01,0.1047389E+00,0.1410312E+00,0.3144056E+00,0.8223500E+00,0.8705581E+00,0.2073945E+00,0.2937436E-01,0.4418871E-02,0.6247617E-03,0.1188419E-03,0.2145310E-04,0.3793062E-04,0.4487192E-03,0.4651812E-01,0.4650101E-02,0.9288612E-01,0.9784722E-01,0.2775293E-01,0.6115219E-01,
>
> 0.9568319E-01,0.3376210E-01,0.4138662E-01,0.2358972E-01,0.1969589E-01,0.2843634E-01,0.2342540E-01,0.8103615E-02,0.3448099E-02,0.7816843E-02,0.9035400E-02,0.6497264E-02,0.2002942E-01,0.1395259E+00,0.2419057E+00,0.6914363E-01,0.6008841E-01,0.4781502E-01,0.5576929E-02,0.5630654E-03,0.4329984E-02,0.7708156E-02,0.9577785E-02,0.2238534E-01,0.2953239E-01,0.2761710E-01,0.3779950E-01,0.3350705E-01,0.4099815E-01,0.5825405E-01,0.5754150E-01,0.6740542E-01,0.7738841E-01,0.9395976E-01,0.1125712E+00,0.1853697E+00,0.2867987E+00,0.2497954E+00,0.8689786E-01,0.1633998E-01,0.2972100E-02,0.4797695E-03,0.9907543E-04,0.2630055E-04,0.4326147E-04,0.3915109E-03,0.4403545E-01,0.7822213E-02,0.1077334E+00,0.9601451E-01,0.2970043E-01,0.6044319E-01,
>
> 0.6726907E-01,0.3973379E-01,0.4293448E-01,0.2719793E-01,0.2464437E-01,0.3161275E-01,0.2721916E-01,0.1279970E-01,0.8913528E-02,0.1111176E-01,0.1399803E-01,0.1374758E-01,0.4348021E-01,0.2415760E+00,0.3361519E+00,0.4495556E-01,0.3283257E-01,0.2604829E-01,0.6277758E-02,0.1306557E-02,0.6591999E-02,0.1515555E-01,0.5407024E-02,0.3577704E-02,0.3391261E-02,0.9091800E-02,0.3680739E-01,0.1161744E+00,0.1314589E+00,0.1119177E+00,0.9781276E-01,0.9123342E-01,0.8549694E-01,0.8449271E-01,0.8546127E-01,0.9224446E-01,0.7830644E-01,0.6116871E-01,0.2562409E-01,0.7241758E-02,0.1713816E-02,0.3386083E-03,0.8170580E-04,0.2688931E-04,0.4080362E-04,0.3391886E-03,0.4114522E-01,0.1512450E-01,0.1287728E+00,0.9372337E-01,0.3146907E-01,0.5930399E-01,
>
> 0.4348354E-01,0.5008505E-01,0.4569119E-01,0.3135827E-01,0.3470170E-01,0.3984557E-01,0.3394159E-01,0.2304775E-01,0.3091791E-01,0.1793484E-01,0.2416584E-01,0.3275779E-01,0.1272687E+00,0.4726028E+00,0.3376644E+00,0.2599599E-01,0.1546338E-01,0.1059649E-01,0.7716056E-02,0.3790854E-02,0.1462767E-01,0.3513828E-01,0.2741873E-02,0.3088193E-03,0.1875218E-03,0.2618570E-02,0.5382948E-01,0.5677117E+00,0.7196847E+00,0.2989115E+00,0.1548267E+00,0.1045880E+00,0.7731146E-01,0.5816138E-01,0.4743601E-01,0.3096705E-01,0.1780926E-01,0.1099431E-01,0.5715238E-02,0.2362370E-02,0.7833299E-03,0.2072879E-03,0.5806183E-04,0.1971219E-04,0.3364425E-04,0.2923626E-03,0.3849743E-01,0.3147306E-01,0.1735495E+00,0.9093874E-01,0.3217609E-01,0.5213468E-01,
>
> 0.3434681E-01,0.5781734E-01,0.5847079E-01,0.4302881E-01,0.3907984E-01,0.5918604E-01,0.5853324E-01,0.5683731E-01,0.8518384E-01,0.3769261E-01,0.5476855E-01,0.1028494E+00,0.5784550E+00,0.8782803E+00,0.1842485E+00,0.1272144E-01,0.7322852E-02,0.6657565E-02,0.1269100E-01,0.2443688E-01,0.2280098E-01,0.3138442E-01,0.1557003E-02,0.1445880E-03,0.1064606E-03,0.1557382E-02,0.9419733E-01,0.2483386E+01,0.2825258E+01,0.6148053E+00,0.1762852E+00,0.7715444E-01,0.4212902E-01,0.2130515E-01,0.1281251E-01,0.4857580E-02,0.1488796E-02,0.9052465E-03,0.5064251E-03,0.3255549E-03,0.1811494E-03,0.8082773E-04,0.3742215E-04,0.1403465E-04,0.3321273E-04,0.2598714E-03,0.3770767E-01,0.5739681E-01,0.3059424E+00,0.8195814E-01,0.3618274E-01,0.5168135E-01,
>
> 0.4512117E-01,0.5308073E-01,0.7272963E-01,0.4508510E-01,0.3880418E-01,0.1067216E+00,0.1105220E+00,0.7333525E-01,0.7108819E-01,0.5104649E-01,0.7576904E-01,0.1618577E+00,0.1010573E+01,0.7501726E+00,0.6456181E-01,0.1009602E-01,0.7550230E-02,0.7154825E-02,0.1866851E-01,0.8297640E-01,0.1359374E-01,0.6905828E-02,0.1773696E-02,0.1453961E-02,0.1664439E-02,0.4455402E-02,0.1107461E+00,0.1685297E+01,0.2000099E+01,0.5046504E+00,0.1229368E+00,0.4560262E-01,0.2151676E-01,0.9012380E-02,0.4638202E-02,0.1131346E-02,0.2825367E-03,0.1558241E-03,0.1162573E-03,0.9299316E-04,0.7714150E-04,0.4996025E-04,0.2825607E-04,0.1616293E-04,0.3807035E-04,0.2929824E-03,0.4176042E-01,0.4949051E-01,0.4732706E+00,0.7777966E-01,0.3910970E-01,0.5760659E-01,
>
> 0.4790651E-01,0.5215655E-01,0.7854529E-01,0.4591308E-01,0.3749941E-01,0.1251934E+00,0.1302131E+00,0.7672322E-01,0.6730324E-01,0.5428094E-01,0.8278674E-01,0.1791058E+00,0.1099262E+01,0.6625134E+00,0.4602054E-01,0.9400134E-02,0.7094789E-02,0.6855966E-02,0.2054928E-01,0.1278550E+00,0.1072797E-01,0.3905681E-02,0.1861920E-02,0.2373821E-02,0.3431623E-02,0.6137554E-02,0.1084819E+00,0.1793933E+01,0.2067821E+01,0.4455505E+00,0.1021708E+00,0.3657212E-01,0.1675937E-01,0.6748084E-02,0.3349477E-02,0.7342977E-03,0.1548510E-03,0.8225830E-04,0.6979348E-04,0.6679614E-04,0.6250224E-04,0.4505874E-04,0.2710121E-04,0.1526862E-04,0.3366982E-04,0.3120460E-03,0.4289762E-01,0.4583776E-01,0.5531438E+00,0.7605083E-01,0.4017489E-01,0.5781028E-01,
>
> 0.5279086E-01,0.5032045E-01,0.9581324E-01,0.4837249E-01,0.3563856E-01,0.1632324E+00,0.1840193E+00,0.8057278E-01,0.6066804E-01,0.5927094E-01,0.9734113E-01,0.2154147E+00,0.1165791E+01,0.4914498E+00,0.2073040E-01,0.8098976E-02,0.5923700E-02,0.5595430E-02,0.2332763E-01,0.2208645E+00,0.8225261E-02,0.1758009E-02,0.2555679E-02,0.6282730E-02,0.1250278E-01,0.1153100E-01,0.7383288E-01,0.1936265E+01,0.2125543E+01,0.2634855E+00,0.5367346E-01,0.2204685E-01,0.9842224E-02,0.3805011E-02,0.1883615E-02,0.3511793E-03,0.4528087E-04,0.1905631E-04,0.2738672E-04,0.4089296E-04,0.4579933E-04,0.2804448E-04,0.2116800E-04,0.1238481E-04,0.3598009E-04,0.3693456E-03,0.4379154E-01,0.3980652E-01,0.7051615E+00,0.7173526E-01,0.4188387E-01,0.5828505E-01,
>
> 0.5191508E-01,0.5968262E-01,0.1211241E+00,0.5550964E-01,0.3769149E-01,0.1756077E+00,0.1964576E+00,0.9167825E-01,0.6101080E-01,0.7878678E-01,0.1071014E+00,0.2325657E+00,0.1161060E+01,0.2543801E+00,0.1150904E-01,0.7423799E-02,0.5517107E-02,0.6101956E-02,0.2248316E-01,0.1662242E+00,0.4561628E-02,0.7951080E-03,0.1915489E-02,0.4712472E-02,0.6888872E-02,0.9463635E-02,0.5949265E-01,0.5709903E+00,0.5442112E+00,0.9465075E-01,0.2232800E-01,0.8378362E-02,0.3942512E-02,0.1611552E-02,0.7982667E-03,0.1764814E-03,0.2541264E-04,0.1421835E-04,0.1722847E-04,0.2776491E-04,0.2751415E-04,0.3112745E-04,0.3019421E-04,0.2355444E-04,0.6594906E-04,0.4409241E-03,0.3780599E-01,0.4199751E-01,0.7580482E+00,0.6412683E-01,0.4505653E-01,0.5444700E-01,
>
> 0.4355716E-01,0.8733416E-01,0.1536599E+00,0.6560311E-01,0.4475639E-01,0.1441608E+00,0.1640185E+00,0.9005611E-01,0.6951861E-01,0.8253337E-01,0.1129339E+00,0.2554535E+00,0.7104740E+00,0.1192783E+00,0.7677205E-02,0.5773198E-02,0.6082294E-02,0.6760821E-02,0.1740481E-01,0.5968277E-01,0.2647498E-02,0.5596312E-03,0.1288748E-02,0.7751010E-03,0.8822999E-03,0.3176315E-02,0.2331444E-01,0.7107395E-01,0.4164753E-01,0.1642068E-01,0.6303825E-02,0.2976200E-02,0.1612777E-02,0.7591693E-03,0.4160732E-03,0.1108386E-03,0.2556193E-04,0.1785732E-04,0.1403372E-04,0.1937110E-04,0.2752077E-04,0.4076820E-04,0.5958297E-04,0.1021522E-03,0.2122076E-03,0.6466907E-03,0.3016595E-01,0.5712829E-01,0.6557645E+00,0.5549022E-01,0.5029169E-01,0.5510478E-01,
>
> 0.3504792E-01,0.1205022E+00,0.1855510E+00,0.7573815E-01,0.5111053E-01,0.1125733E+00,0.1269330E+00,0.8580401E-01,0.7008798E-01,0.8702227E-01,0.1132957E+00,0.2337266E+00,0.3684450E+00,0.4996030E-01,0.4962245E-02,0.4465351E-02,0.5276125E-02,0.6395389E-02,0.1176227E-01,0.1897196E-01,0.2114873E-02,0.5134127E-03,0.1242211E-02,0.1699282E-03,0.9218296E-04,0.1476873E-02,0.1088367E-01,0.9746083E-02,0.5051154E-02,0.3829346E-02,0.2159273E-02,0.1257526E-02,0.7798878E-03,0.4205173E-03,0.2524269E-03,0.7814084E-04,0.2024509E-04,0.1311006E-04,0.1400313E-04,0.1909558E-04,0.2873475E-04,0.5047805E-04,0.1034147E-03,0.3474390E-03,0.4601162E-03,0.7914868E-03,0.2259713E-01,0.6977478E-01,0.5253417E+00,0.4968414E-01,0.5419997E-01,0.4578392E-01,
>
> 0.3644789E-01,0.1094320E+00,0.1812871E+00,0.8090766E-01,0.5458792E-01,0.8453539E-01,0.9459070E-01,0.7997739E-01,0.7218789E-01,0.9252239E-01,0.1235988E+00,0.1916332E+00,0.1600187E+00,0.1689259E-01,0.2238496E-02,0.3474775E-02,0.4246865E-02,0.4639725E-02,0.7562945E-02,0.8655204E-02,0.3218143E-02,0.1432361E-02,0.2987746E-02,0.6836831E-03,0.4578503E-03,0.2736592E-02,0.9189470E-02,0.5100553E-02,0.2354053E-02,0.1784556E-02,0.1184593E-02,0.7540732E-03,0.5015396E-03,0.2932982E-03,0.1871755E-03,0.6034159E-04,0.2081635E-04,0.1371099E-04,0.1420733E-04,0.1977560E-04,0.2733073E-04,0.4674405E-04,0.8729282E-04,0.2633706E-03,0.3779292E-03,0.6857174E-03,0.1703467E-01,0.6216488E-01,0.2697471E+00,0.4844839E-01,0.5147192E-01,0.4584366E-01,
>
> 0.4538495E-01,0.6428195E-01,0.1408200E+00,0.7812387E-01,0.5196887E-01,0.6079464E-01,0.6857989E-01,0.7737246E-01,0.7458069E-01,0.9937596E-01,0.1602969E+00,0.1203187E+00,0.4573522E-01,0.4332148E-02,0.5447117E-03,0.2177756E-02,0.2675605E-02,0.3269132E-02,0.4752257E-02,0.4392496E-02,0.1062331E-01,0.1656332E-01,0.2009005E-01,0.6549678E-01,0.8628400E-01,0.4396579E-01,0.1700950E-01,0.7778500E-02,0.3640412E-02,0.1845144E-02,0.1257515E-02,0.6835750E-03,0.4559327E-03,0.2127205E-03,0.1413101E-03,0.5491232E-04,0.2119532E-04,0.1395551E-04,0.1418517E-04,0.1809416E-04,0.2645285E-04,0.2941790E-04,0.3840854E-04,0.4570941E-04,0.9958091E-04,0.3501749E-03,0.1081307E-01,0.3869956E-01,0.7878832E-01,0.5268212E-01,0.4303591E-01,0.4563945E-01,
>
> 0.5201690E-01,0.4414320E-01,0.1032115E+00,0.7373682E-01,0.5161830E-01,0.4897009E-01,0.5589062E-01,0.7558857E-01,0.7586090E-01,0.9935160E-01,0.1867155E+00,0.8613827E-01,0.1779336E-01,0.1280139E-02,0.2480400E-03,0.1470505E-02,0.1857580E-02,0.2378422E-02,0.3442274E-02,0.3341513E-02,0.2902278E-01,0.1017186E+00,0.1262739E+00,0.1368409E+01,0.3036598E+01,0.3310379E+00,0.2551310E-01,0.1309108E-01,0.5509888E-02,0.2353525E-02,0.1199714E-02,0.7416096E-03,0.5203202E-03,0.2366006E-03,0.1637158E-03,0.7003936E-04,0.2485285E-04,0.1414354E-04,0.1572256E-04,0.1952128E-04,0.2294974E-04,0.2407856E-04,0.2091599E-04,0.1400297E-04,0.3944428E-04,0.2032251E-03,0.8302215E-02,0.2844989E-01,0.3576908E-01,0.5440925E-01,0.3841882E-01,0.4544563E-01,
>
> 0.5348080E-01,0.5381636E-01,0.7901072E-01,0.6067793E-01,0.6508222E-01,0.5235143E-01,0.6034862E-01,0.8425818E-01,0.8353563E-01,0.9121266E-01,0.1426220E+00,0.7975989E-01,0.1563352E-01,0.1676089E-02,0.5479662E-03,0.1221801E-02,0.1447088E-02,0.2353543E-02,0.3522641E-02,0.6034007E-02,0.2974472E-01,0.6227752E-01,0.1050294E+00,0.1745723E+00,0.3152606E+00,0.1936336E+00,0.3980809E-01,0.1048688E-01,0.3411023E-02,0.2170874E-02,0.1468369E-02,0.1162395E-02,0.9268029E-03,0.7087942E-03,0.2933519E-03,0.1954615E-03,0.1320647E-03,0.4318345E-04,0.3817508E-04,0.3576300E-04,0.4651525E-04,0.4589578E-04,0.3474174E-04,0.3873098E-04,0.8491570E-04,0.2805419E-03,0.6853279E-02,0.3381941E-01,0.5052855E-01,0.5274747E-01,0.4520723E-01,0.4574984E-01,
>
> 0.6701137E-01,0.6752905E-01,0.7728650E-01,0.5609375E-01,0.6649279E-01,0.7695681E-01,0.8341227E-01,0.7343599E-01,0.8764342E-01,0.8050204E-01,0.1281606E+00,0.8062495E-01,0.2190935E-01,0.4653459E-02,0.1494571E-02,0.1400683E-02,0.1561912E-02,0.2936219E-02,0.4495304E-02,0.2470658E-01,0.1016565E-01,0.1789034E-01,0.2605223E-01,0.2764470E-01,0.8395130E-01,0.1121967E-01,0.1370545E-01,0.9507278E-02,0.2042823E-02,0.1674927E-02,0.1423714E-02,0.1752843E-02,0.1606748E-02,0.1452867E-02,0.6601716E-03,0.5640681E-03,0.3646661E-03,0.7374791E-04,0.1665530E-03,0.1576760E-03,0.9556996E-04,0.1555420E-03,0.1590785E-03,0.6956725E-04,0.1766699E-03,0.2260075E-03,0.6065029E-02,0.2681332E-01,0.5102418E-01,0.4827190E-01,0.5395776E-01,0.5384292E-01,
>
> 0.6821958E-01,0.6653726E-01,0.7542274E-01,0.7230803E-01,0.6761272E-01,0.7240158E-01,0.7374620E-01,0.8704310E-01,0.8197942E-01,0.8302580E-01,0.7716634E-01,0.3827731E-01,0.9503171E-02,0.7745612E-02,0.3412459E-02,0.3226514E-02,0.3284226E-02,0.2882712E-02,0.7848977E-02,0.7891713E-02,0.3232360E-01,0.1524918E-01,0.1142893E-01,0.1129006E-01,0.1630594E-01,0.1556832E-01,0.4080241E-02,0.3704535E-02,0.4824599E-02,0.4622310E-02,0.4481132E-02,0.4393321E-02,0.1195263E-02,0.3901052E-03,0.3871015E-03,0.3799590E-03,0.3708666E-03,0.4182171E-03,0.4100097E-03,0.4008761E-03,0.3923923E-03,0.6484714E-03,0.6547080E-03,0.6695485E-03,0.1108473E-02,0.1177582E-02,0.4044055E-02,0.3956076E-01,0.6226875E-01,0.4705076E-01,0.5098740E-01,0.6245828E-01,
>
> 0.8092227E-01,0.7079308E-01,0.5342720E-01,0.7015040E-01,0.6996244E-01,0.6826767E-01,0.6489854E-01,0.6764144E-01,0.7897275E-01,0.7333872E-01,0.5892679E-01,0.3420132E-01,0.2084583E-01,0.8810038E-02,0.6224935E-02,0.6325575E-02,0.5697807E-02,0.5838716E-02,0.1256433E-01,0.1533775E-01,0.1439642E-01,0.2953991E-01,0.3538589E-01,0.3428081E-01,0.3253277E-01,0.1729425E-01,0.9545859E-02,0.6373140E-02,0.3294692E-02,0.2472626E-02,0.2706315E-02,0.2923991E-02,0.3116356E-02,0.3363619E-02,0.3564003E-02,0.3989297E-02,0.2756725E-02,0.2779364E-02,0.2624145E-02,0.2345233E-02,0.2014759E-02,0.1687530E-02,0.1476065E-02,0.1200352E-02,0.2155286E-02,0.2989681E-02,0.2095317E-01,0.3906499E-01,0.5602789E-01,0.5702719E-01,0.5946530E-01,0.5269162E-01,
>
> 0.8286110E-01,0.8168069E-01,0.8009119E-01,0.6620687E-01,0.5965812E-01,0.5684869E-01,0.5182595E-01,0.5181318E-01,0.5585941E-01,0.4079785E-01,0.4666164E-01,0.3397221E-01,0.2914809E-01,0.2380620E-01,0.2234047E-01,0.2262664E-01,0.2265226E-01,0.2268667E-01,0.1174864E-01,0.1307145E-01,0.1285912E-01,0.9885800E-02,0.7340377E-02,0.1148781E-01,0.1142617E-01,0.9100934E-02,0.8887787E-02,0.1076426E-01,0.1067490E-01,0.1006696E-01,0.1005100E-01,0.1004431E-01,0.1004119E-01,0.1003923E-01,0.1003862E-01,0.1003851E-01,0.1184939E-01,0.1184676E-01,0.1184078E-01,0.9852745E-02,0.9845543E-02,0.9843818E-02,0.9854952E-02,0.1435728E-01,0.1446773E-01,0.1388169E-01,0.1355964E-01,0.2942844E-01,0.5460724E-01,0.6521013E-01,0.7582597E-01,0.6257769E-01,
>
> 0.8771370E-01,0.8181528E-01,0.7944667E-01,0.7697722E-01,0.7199347E-01,0.7784528E-01,0.8212147E-01,0.8462927E-01,0.8066515E-01,0.8133087E-01,0.7760409E-01,0.7819823E-01,0.7600749E-01,0.5320767E-01,0.5417101E-01,0.5467827E-01,0.5436387E-01,0.4824388E-01,0.4374537E-01,0.3483792E-01,0.3009322E-01,0.3343575E-01,0.4591080E-01,0.5441335E-01,0.5488610E-01,0.5580195E-01,0.5175355E-01,0.4840761E-01,0.4244217E-01,0.4241133E-01,0.4272235E-01,0.4309691E-01,0.4346220E-01,0.4396723E-01,0.4440640E-01,0.4548217E-01,0.4686106E-01,0.5463197E-01,0.5586066E-01,0.5728871E-01,0.5872028E-01,0.6020900E-01,0.6119667E-01,0.6130693E-01,0.7987627E-01,0.7459150E-01,0.5077481E-01,0.8016669E-01,0.6848454E-01,0.7942270E-01,0.8315815E-01,0.8296445E-01,
>
> 0.7679867E-01,0.7674027E-01,0.6930419E-01,0.5388183E-01,0.5383458E-01,0.5734111E-01,0.5738815E-01,0.5235089E-01,0.5239711E-01,0.5243716E-01,0.6887815E-01,0.7278056E-01,0.7093403E-01,0.7079363E-01,0.6956196E-01,0.4564178E-01,0.4561833E-01,0.4559211E-01,0.4547448E-01,0.4837488E-01,0.4495171E-01,0.4502418E-01,0.4065792E-01,0.4074185E-01,0.4076913E-01,0.4081822E-01,0.6037227E-01,0.6041060E-01,0.6043662E-01,0.6045376E-01,0.6046853E-01,0.6047970E-01,0.5194132E-01,0.7593730E-01,0.7594074E-01,0.7595118E-01,0.7596734E-01,0.7597808E-01,0.7599201E-01,0.7600705E-01,0.7601944E-01,0.5757346E-01,0.5756919E-01,0.5755191E-01,0.5751974E-01,0.5746212E-01,0.6995818E-01,0.6608821E-01,0.8221844E-01,0.7311027E-01,0.7306487E-01,0.7300587E-01,
>
> 0.1488105E+00,0.8080313E-01,0.4571154E-01,0.4635561E-01,0.4475472E-01,0.2700748E-01,0.2637547E-01,0.3611866E-01,0.4028207E-01,0.4128201E-01,0.4139032E-01,0.4522617E-01,0.7106137E-01,0.6954663E-01,0.6308198E-01,0.9954325E-01,0.1133250E+00,0.1178772E+00,0.9240813E-01,0.8103993E-01,0.4973169E-01,0.4093944E-01,0.6562766E-01,0.8338632E-01,0.8787832E-01,0.7978603E-01,0.6623685E-01,0.5602208E-01,0.5580108E-01,0.6166632E-01,0.6674261E-01,0.7047142E-01,0.7348517E-01,0.7726602E-01,0.8042598E-01,0.8832898E-01,0.9922590E-01,0.1021007E+00,0.9744413E-01,0.9089941E-01,0.8604769E-01,0.8247492E-01,0.8056927E-01,0.7934239E-01,0.7672020E-01,0.7089683E-01,0.7225678E-01,0.7808220E-01,0.7029370E-01,0.6940968E-01,0.5973309E-01,0.6364381E-01,
>
> 0.5795809E-01,0.5495343E-01,0.4706498E-01,0.4043636E-01,0.3340418E-01,0.4968254E-01,0.4446836E-01,0.4453667E-01,0.4226248E-01,0.5675332E-01,0.5687277E-01,0.5951097E-01,0.7391368E-01,0.8420131E-01,0.7043754E-01,0.8099187E-01,0.8097482E-01,0.8095984E-01,0.8691131E-01,0.6951261E-01,0.6299660E-01,0.6670299E-01,0.6703420E-01,0.6357251E-01,0.6357595E-01,0.6960068E-01,0.6960717E-01,0.7648194E-01,0.7649686E-01,0.6949229E-01,0.6949931E-01,0.6950378E-01,0.6950686E-01,0.6950989E-01,0.6951174E-01,0.6955834E-01,0.6955715E-01,0.6955381E-01,0.6954642E-01,0.6976538E-01,0.6974361E-01,0.6971327E-01,0.6968047E-01,0.6601531E-01,0.6596813E-01,0.6120099E-01,0.6473158E-01,0.6225431E-01,0.6083302E-01,0.5828087E-01,0.5924804E-01,0.5926056E-01,
>
> 0.4313182E-01,0.3932900E-01,0.4598860E-01,0.4888564E-01,0.4748531E-01,0.4642788E-01,0.4995710E-01,0.5042158E-01,0.5906870E-01,0.6690776E-01,0.7476615E-01,0.6681214E-01,0.6951392E-01,0.7723163E-01,0.8187254E-01,0.8067562E-01,0.8044949E-01,0.6721825E-01,0.6651182E-01,0.7611940E-01,0.7749233E-01,0.7050008E-01,0.7061531E-01,0.7043605E-01,0.6766590E-01,0.6739765E-01,0.6707212E-01,0.6718514E-01,0.6485243E-01,0.6530887E-01,0.6574852E-01,0.6608203E-01,0.6634366E-01,0.6664646E-01,0.6686886E-01,0.6728150E-01,0.6754559E-01,0.6815954E-01,0.6797798E-01,0.6752418E-01,0.6685620E-01,0.6594741E-01,0.6504641E-01,0.6398232E-01,0.6165998E-01,0.6054787E-01,0.5572481E-01,0.5542183E-01,0.5286020E-01,0.5496611E-01,0.5387216E-01,0.5257342E-01,
>
> 0.4266687E-01,0.3810273E-01,0.3690985E-01,0.3896502E-01,0.3859067E-01,0.3862647E-01,0.3433147E-01,0.3431797E-01,0.3826348E-01,0.3831679E-01,0.3952323E-01,0.3967597E-01,0.4049171E-01,0.5028000E-01,0.5332739E-01,0.5319454E-01,0.5323796E-01,0.5383471E-01,0.5444790E-01,0.6002122E-01,0.6141461E-01,0.6374808E-01,0.6341618E-01,0.6295031E-01,0.6453129E-01,0.6409577E-01,0.6352402E-01,0.6320708E-01,0.7022721E-01,0.7017468E-01,0.7016072E-01,0.7016493E-01,0.7017563E-01,0.7019608E-01,0.7021748E-01,0.7027970E-01,0.7037533E-01,0.7679827E-01,0.7689483E-01,0.7700915E-01,0.7711989E-01,0.7722896E-01,0.7730565E-01,0.7735583E-01,0.6871950E-01,0.6865267E-01,0.7006264E-01,0.6946852E-01,0.6585055E-01,0.6492389E-01,0.6023150E-01,0.6041276E-01,
>
> 0.4844484E-01,0.3642741E-01,0.3376257E-01,0.3720178E-01,0.3909988E-01,0.3859439E-01,0.3467487E-01,0.3462980E-01,0.3797071E-01,0.3823806E-01,0.3942306E-01,0.3967370E-01,0.4195464E-01,0.5007371E-01,0.4961523E-01,0.4843083E-01,0.4845073E-01,0.4904005E-01,0.5135238E-01,0.6128201E-01,0.6851908E-01,0.7265307E-01,0.6970225E-01,0.6703834E-01,0.6760444E-01,0.6559991E-01,0.6312431E-01,0.6185145E-01,0.6704390E-01,0.6698888E-01,0.6707250E-01,0.6718721E-01,0.6730217E-01,0.6746334E-01,0.6760543E-01,0.6796423E-01,0.6846146E-01,0.7383548E-01,0.7432497E-01,0.7492594E-01,0.7553371E-01,0.7615437E-01,0.7659570E-01,0.7685904E-01,0.6966449E-01,0.6909129E-01,0.6832164E-01,0.7162107E-01,0.6921861E-01,0.6137763E-01,0.6072355E-01,0.6161803E-01,
>
> 0.3875395E-01,0.3867135E-01,0.3861343E-01,0.3863410E-01,0.3866171E-01,0.3867369E-01,0.3414714E-01,0.3414922E-01,0.3834129E-01,0.3834565E-01,0.3956075E-01,0.3957593E-01,0.4008044E-01,0.5057959E-01,0.5459315E-01,0.5460132E-01,0.5460802E-01,0.5518256E-01,0.5523199E-01,0.5954568E-01,0.5962998E-01,0.6156283E-01,0.6155794E-01,0.6153870E-01,0.6338990E-01,0.6336696E-01,0.6333258E-01,0.6331021E-01,0.7101443E-01,0.7100755E-01,0.7100359E-01,0.7100160E-01,0.7100052E-01,0.7099973E-01,0.7099944E-01,0.7099982E-01,0.7100176E-01,0.7787250E-01,0.7787547E-01,0.7787941E-01,0.7788356E-01,0.7788800E-01,0.7789153E-01,0.7789456E-01,0.6864183E-01,0.6864185E-01,0.7050855E-01,0.6883876E-01,0.6499121E-01,0.6486494E-01,0.6486809E-01,0.6489217E-01,
>
> 0.3869122E-01,0.3868914E-01,0.3868745E-01,0.3868744E-01,0.3868803E-01,0.3868850E-01,0.3412816E-01,0.3412841E-01,0.3834546E-01,0.3834570E-01,0.3956086E-01,0.3956131E-01,0.4005100E-01,0.5062821E-01,0.5469843E-01,0.5469926E-01,0.5469957E-01,0.5527038E-01,0.5527203E-01,0.5950757E-01,0.5950988E-01,0.6141567E-01,0.6141596E-01,0.6141577E-01,0.6328529E-01,0.6328486E-01,0.6328411E-01,0.6328352E-01,0.7105376E-01,0.7105348E-01,0.7105327E-01,0.7105315E-01,0.7105307E-01,0.7105298E-01,0.7105292E-01,0.7105280E-01,0.7105269E-01,0.7797483E-01,0.7797478E-01,0.7797474E-01,0.7797477E-01,0.7797477E-01,0.7797478E-01,0.7797483E-01,0.6865252E-01,0.6865261E-01,0.7053513E-01,0.6878572E-01,0.6492829E-01,0.6492439E-01,0.6492390E-01,0.6492462E-01),
> nrow = length(z.space), ncol = length(x.space), byrow =
> TRUE,
> dimnames = list(paste0("R",seq(1,length(z.space),by=1)),
> paste0("C",seq(1,length(x.space),by=1))))
>
>
> b2<- vals
> for(k in 1:length(z.space)){
> for(j in 1: length(x.space)){
> if (b2[k,j]==0) {
> b2[k,j]<- NA
> }
> }
> }
>
> vals2<- log10(b2)
> cv<- matrix(as.integer(cut(vals2, breaks=100)), dim(vals2))
> pal<- blue2green2red(100)
>
> plot.new()
> par(mar=c(4,4,1,4))
> plot(1,1, col="white" ,xlim=x.range, ylim=z.range, xaxs="i", yaxs="i",
> las=1, xlab="X", ylab="Z")
>
> #plot the irregular grid cells polygon-by-polygon
> for (k in (1:length(z.space))){
> for (j in (1:length(x.space))){
>
> if (j==1){
> x = c(0, 0, x.space[j], x.space[j])
> } else {
> x = c(sum(x.space[1:(j-1)]), sum(x.space[1:(j-1)]),
> sum(x.space[1:j]), sum(x.space[1:j]))
> }
>
> if (k==1){
> y = c(z.sum, z.sum - sum(z.space[1]), z.sum - sum(z.space[1]), z.sum)
> } else {
> y = c(z.sum - sum(z.space[1:(k-1)]), z.sum - sum(z.space[1:k]),
> z.sum - sum(z.space[1:k]), z.sum - sum(z.space[1:(k-1)]))
> }
>
> polygon(x=x, y=y, border=NA, col=pal[cv[k, j]])
> }
> }
>
> image.plot(log10(b2), legend.only=TRUE, legend.mar=3,legend.lab=c(0.0001,
> 0.001, 0.01, 0.1, 1.0))
>
Hi Eric,
Have a look at sampcolorlegend in the shape package.
Jim
More information about the R-help
mailing list