Greetings, R helpers, I have been using the Tps() function in the fields package to model response surfaces for some residuals data. But i donĀ“t know how i can get more resolution and fill all the shape area? Thank you,
Dominic outresi <- Tps(x=as.matrix(ana),Y=resi) preresi <- predict.surface(outresi) image(preresi) resi <- structure(c(1.18650257765518, 0.922870034874496, 1.40456497838008, -1.98869637540279, 0.652835483390163, 3.38161411979703, 1.5301074186747, -0.409663707199911, 1.34434770758852, -0.416787026362009, 4.18757416663512, 3.28025478124903, -5.17448326255257, 2.95828686941439, -0.782530617455976, 1.82121460192511, 0.203539209775968, 0.508322787328424, -1.94498277323241, -4.89815558191742, 3.36600081231205, -3.43524939015204, -3.05944133924843, 1.19512325854811, -3.04309148543589, 1.40717743356163, -3.35495846881043, 2.72864346101486, 1.37230147197744, -0.106841569047528, -0.119344066154028, -0.608350506308074, 1.44743280454379, 3.30716673292547, 0.874190165066352, -0.916075316748891, 1.37181996156702, -1.2930355920363, -4.54134505448977, 1.80012571483517, -0.721990604703457, -2.65811792871444, 2.14103809342353, -2.41809548989305, -1.81952424936354, 3.31004331724789, -3.30108275773548, 2.9757633493691, 0.295994951412489, 2.30145808678062, 0.422002402195169, 0.817785741134315, 3.49852828825309, 0.409779315387873, 2.30965751639276, -1.00806803466553, 2.84451533239655, -2.35467760469766, 4.42619863347384, 2.70778104139562, 2.16518126265409, -2.88134006056208, -1.09216809797566, -1.52993750262298, 2.65269273615996, 3.15979996437503, -0.363350001435549, -6.39582581459343, -1.04031709170907, 1.33328272103405, -5.49726774890521, 4.08987284624148, 2.50912004072798, 1.8274159447864, -2.22385808283135, -4.86917684419856, -0.564375776024271, 1.8708963744, -3.78120181090081, -0.558708890799885, 4.05539617977911, -0.948307806738612, -1.55729490567548, 1.96565015968813, 2.91601715036197, -1.25289442884952, 2.88976270172768, -1.9916070689498, -1.48771804146596, -2.1835149830146, -2.15786396484737, -7.13845827165779, 2.54353180775172, -3.42663298459986), .Names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "23", "25", "26", "27", "28", "29", "30", "31", "32", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96", "98", "99")) ana <- structure(list(Long = c(521328, 643449, 547859, 495486, 510314, 570154, 651957, 591244, 653608, 616292, 562796, 669916, 648619, 670351, 607180, 519179, 610567, 625025, 631194, 640477, 531377, 618769, 684214, 516223, 623222, 524431, 511201, 560578, 641809, 590324, 497771, 550783, 502146, 648423, 574209, 510859, 549546, 639190, 625827, 585343, 598945, 682612, 489814, 485023, 527749, 560019, 513379, 589733, 623277, 583167, 526634, 613230, 570876, 518427, 660057, 636959, 558736, 505372, 594942, 594370, 532925, 658781, 654995, 615988, 522575, 555930, 639736, 522594, 612930, 576464, 538116, 547280, 566252, 540986, 525267, 498997, 598421, 536101, 516677, 640275, 544158, 597892, 664123, 587674, 518653, 540339, 670165, 632785, 629480, 657933, 523762, 520553, 534026, 585393), Lat = c(4702005, 4734442, 4801512, 4773043, 4700787, 4823868, 4772302, 4724987, 4730789, 4684126, 4696309, 4743224, 4800328, 4703493, 4644830, 4706070, 4832008, 4723028, 4833932, 4680282, 4716165, 4761493, 4690230, 4704557, 4779463, 4701030, 4649389, 4815885, 4706108, 4765782, 4711638, 4801765, 4764251, 4717609, 4644255, 4757557, 4685177, 4824470, 4812691, 4669787, 4786917, 4694030, 4738281, 4747735, 4695489, 4788103, 4798361, 4799728, 4703623, 4751251, 4658771, 4676811, 4718574, 4732806, 4730358, 4664595, 4775147, 4692421, 4688775, 4688943, 4685679, 4782554, 4822868, 4835419, 4706352, 4720490, 4780523, 4676863, 4741204, 4839683, 4799702, 4675170, 4753162, 4702027, 4771652, 4701610, 4777923, 4747354, 4694770, 4695132, 4741500, 4827467, 4691376, 4714920, 4698461, 4797835, 4758256, 4648384, 4645955, 4669863, 4675121, 4675090, 4725092, 4639803)), .Names = c("Long", "Lat"), class = "data.frame", row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "23", "25", "26", "27", "28", "29", "30", "31", "32", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96", "98", "99")) [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo