IACS Compliant Field Numbers

[Robert Crossley]

Hi All,

 

I have a lot of references to the IACS numbers for fields in England, and wish to geolocate them.

 

I have found a reference to be able to do this using the 100K map sheet reference and the grid reference, and if I wanted to spend a few hours on it, I am sure I could be able to do this. 

 

However, I suspect I would be reinventing the wheel.

 

Does anyone one have code or logic to be able to do this.  Happy to provide what ever tool I create back.

 

Regards

Rob.

 

-------------------------------------------

Robert Crossley

Managing Director
Agtrix P/L Australia
 
NSW Office:
Unit 6, 2 Bonanza Drive
Billinudgel NSW 2483
AUSTRALIA
Postal:
PO Box 63
New Brighton 2483
 
P: 61 (0) 2 6680 1309
M: 61 (0)419 718 642

UK: 44 (0)2070975139

E: rob...@agtrix.com
W: www.agtrix.com
S: robertcrossley
 
28°30'14.81"S
153°31'41.79"E
 
Brisbane Office:
109 Milsom St
Cooparoo  4151
Queensland
P: 61 7 3843 3363

 

 

 

 

[KK KTG]

Are these numbers something related to soil sampling sort ...

Kalai Selvan.K

KTG Infotech

GIS (Mineral Exploration)
Geological Mapping (Historical Map) Data Integration, Analysis, Capturing and Data Conversion
for Geochemistry,IP and Drill logs
Skypes-ktgcad

--
--
You received this message because you are subscribed to the
Google Groups "MapInfo-L" group.To post a message to this group, send
email to mapi...@googlegroups.com
To unsubscribe from this group, go to:
http://groups.google.com/group/mapinfo-l/subscribe?hl=en
For more options, information and links to MapInfo resources (searching
archives, feature requests, to visit our Wiki, visit the Welcome page at
http://groups.google.com/group/mapinfo-l?hl=en

---
You received this message because you are subscribed to the Google Groups "MapInfo-L" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mapinfo-l+...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Robert Crossley]

Sorry all,

 

This was a really poorly written request.

 

What I am looking for is if anyone has got the code that converts the Integrated Administration and Control System (IACS) number that is assigned to all fields where crops are grown in the EU and receive subsidies.  The number is a structured number that can be logically deciphered to give a georeference.  Somewhere there is supposed to be a repository of the spatial extent of each field as well.

 

I would imagine that soil samples may also be related to these numbers.  The diagram below  shows how these are calculated for the IACS number of NY12433412  It seems that it would be relatively easy to decipher except I don’t know enough about the numbering system for the 100K sheets and the projection for the grid reference they use, without a lot of thought and research.  I would appreciate if someone could point me to a worked example, or even better, so mapbasic code.

 

BONUS if someone could tell me a source for the actual spatial information of the field!  I have requested to the Agriculture and Rural Development of the European Commission, but expect that that request could take a while….

 

Rob.

 

 

If you require specific IACS or Rural Land Register field nu

[Seb]

Hi Rob, 

It is relatively easy to write this formula. I've used a variety of formulas,  though usually within excel or MSAccess. You may or may not want to use a lookup table to cross check the letter code. 

How do you want your final reference to look like? e.g., bg 123456,123456

Here is an example of a formula from Exel which creates the easting from a number such as SP605490 

 =CONCATENATE((LEFT(LOOKUP(LEFT(A6,2),Nat_grid_letters!$A:$A,Nat_grid_letters!$B:$B),1)),MID(A6,3,3),"00")  where A6 is the ref to the cell containing  to your field reference and Nat_grid_letters!$A:$A,Nat_grid_letters!$B:$B is the range and separate sheet containing the lookup values for the brit grid.

You can modify the above to fit your field reference, your ouput and software environment.

"BONUS if someone could tell me a source for the actual spatial information of the field!" Do you mean the parcel topology etc...? You will need to speak  to the Rural Payments Agency in the UK to obtain that, https://www.gov.uk/government/organisations/rural-payments-agency.

If you need any further advice please feel free to contact me directly.

HTH 

Cheers

Seb Hudson
+44 (0) 1436 676580
+44 (0) 7811 268789

[Robert Crossley]

Thanks for the leads Seb.

 

I will give them a call.

 

Rob.

[Seb]

No problems Rob. If you also wanted data for other parts of the UK such as Wales. You would need to contacts them seperately.

Cheers

[Robert Crossley]

Tom,

 

I found that I needed to go down to 1 km sheet identifiers, and it all seems relatively easy except to get the lower left hand coordinates of the 100km Letter.  I can them construct the coordinates from the rest of the number. 

 

I merged the data that you gave me on the first 2 letters, set the projection to be the British National Grid and looked for the minimum x and min y of the objects, but they don’t look right.

 

I could round them possibly, but am not sure if this is the way to go.  Any suggestions?

 

Seb, the lookup table that you referred to would probably work?

 

R

 

HP          399999.99            1200081.08

HT           300038.76            1098657.08

HU          399900.2              1100075.32

HW         100537.93            993547.21

HX          200215.97            996117.52

HY           300036.07            998704.1

HZ           400000  1000070.02

NA          924.935 891321.93

NB          100501.07            893764.24

NC          200201.14            896275.97

ND          300033.56            898780.59

NF          866.989 791660.24

NG         100467.87            793925.75

NI           300025.94            598867.98

NJ           300052.17            798788.88

NK          399901.03            800060.79

NL           805.011 691939.52

NM        100440.04            694122.39

NN         200175.28            696496.41

NO         300029.34            698891.59

NR          100411.54            594309.39

NS          200163.76            596612.53

NU         399901.7              600050.81

NW        121332.3              491505.73

NX          200153.03            496658.18

NY          300025.58            498904.88

NZ          399999.98            500051.82

SC           200146.03            396831.47

SD           300023.94            398869.25

SE           400000.02            400046.28

SH           200133.98            296859.15

SJ            300023.88            298974.25

SK           400000.02            300043.35

SM         100309.13            195052.14

SN          200128.52            197029.06

SO          300021.02            199007.03

SP           399999.98            200040.62

SR           100294.02            95208.75

SS           200117.74            97048.62

ST           300019.72            99038.37

SU          400000.02            100038.02

SV           512.46   -6474.73

SW         100271.59            -4636.97

SX           200111.49            -2787.96

SY           300018.44            -931.084

SZ           399999.81            35.679

TA           497877.59            400185.3

TF           497945.27            300177.45

TG          595794.27            300389.03

TL            498009.95            200162.34

TM         596022.52            200364.44

TQ          498071.65            100148.5

TR           596241.33            100341.59

TV           498130.42            146.477

 

From: Thomas Bacon [mailto:Thomas...@mouchel.com]
Sent: Tuesday, 3 February 2015 9:54 PM
To: Robert Crossley
Subject: RE: [MI-L] IACS Compliant Field Numbers

 

Rob,

 

Yes that’s right, the cells are referenced by their lower left corner and are referenced to the 1K sheet (in 5k increments). I thought I had the code I used to generate this kicking around somewhere but unfortunately I can’t seem to find it. Anyway, I hope this helps.

 

BR,

 

Tom Bacon

GIS Engineer, Mouchel

T 01444 472380 │ E thomas...@mouchel.com │ W www.mouchel.com

Our values: innovation │ excellence │ integrity │ responsibility

 

From: Robert Crossley [mailto:rob...@agtrix.com]
Sent: 03 February 2015 10:50
To: Thomas Bacon
Subject: RE: [MI-L] IACS Compliant Field Numbers

 

Thanks Thomas,

 

This looks pretty close to what is needed.  It seems that the grid you have produced is the grid down to the 1K sheet of the code.

 

It seems that if I open your table, and project it in the British National Grid, the grid cell TF4015 represents the cell with the bottom left corner as 540000,  315000.

 

I could use the file that you have provided to find the coordinates of the lower left corner of the 1KM sheet, and then could pull the next two numbers from the ACIS number to get the rest of the x and y coordinates to find the 10m grid.

 

Great lead.

 

Rob.

 

 

 

From: Thomas Bacon [mailto:Thomas...@mouchel.com]
Sent: Tuesday, 3 February 2015 9:19 PM
To: Robert Crossley
Subject: RE: [MI-L] IACS Compliant Field Numbers

 

Hi Rob,

 

I have a 5k UK grid table with the National Grid references in both this format (NA09NE) and this (NA0595) which I made a while back. I’ve attached a copy if it helps. It’s not quite down to the 10m resolution you’re looking for but it might be a start?

 

Regards,

 

Tom Bacon

GIS Engineer, Mouchel

T 01444 472380 │ E thomas...@mouchel.com │ W www.mouchel.com

Our values: innovation │ excellence │ integrity │ responsibility

Mouchel Limited (Mouchel) is registered in England and Wales with registered number 01686040 at Export House, Cawsey Way, Woking, Surrey, UK, GU21 6QX.  The information in this e-mail is confidential and may be legally privileged. It is intended solely for the addressee. Access to this email by anyone else is unauthorised. Any views or opinions expressed in this e-mail may be solely those of the author and are not necessarily those of Mouchel. No contracts may be concluded on behalf of Mouchel by means of email communications. Mouchel reserves the right to monitor and intercept emails sent and received on our network.

 

[Seb]

The lookup table would work and to get you started here are the letter codes and matching numbers.

cheers

NA	09
NB	19
ND	39
NF	08
NG	18
NH	28
NJ	38
NJ	38
NK	48
NM	17
NN	27
NO	37
NP	47
NR	49
NR	16
NS	26
NT	36
NU	46
NV	29
NW	15
NX	25
NY	35
NZ	45
SC	24
SD	34
SE	44
SH	23
SJ	33
SK	43
SM	12
SN	22
SO	32
SP	42
SR	11
SS	21
ST	31
SU	41
SV	00
SW	10
SX	20
SY	30
SZ	40
TA	54
TF	53
TG	63
TL	52
TM	62
TQ	51
TR	61
TV	50

[Thomas Bacon]

Rob,

 

Yes, that’s definitely not quite right. I remember how I made this now – I used the grid maker tool that comes with MapInfo. I’ve just tried it again and it seems to not fit itself perfectly to the extent coordinates I feed into it. Strange.

 

Seb’s lookup table is probably the easiest way to go. I’m going to have to take a look at this at some point and just write some code to generate my own grid rather than using grid maker I guess.

[Tony Pilkington]

Part of an old and bold package.

Modify as you see fit

 

Tony

No virus found in this message.
Checked by AVG - www.avg.com
Version: 2014.0.4800 / Virus Database: 4257/9049 - Release Date: 02/03/15

[Robert Crossley]

Thanks,

 

The key function that I need now is to get the minx and miny values of the 100K tile referred by the 2 letter code.  Currently I am just getting it from the objects of the map sheets that I created from the file that Tom gave me, which is not quite accurate.

 

I can see the answer is in the supplied code, especially after reading the brochure of how it is constructed, but it just might be too late at night to reverse the code to get the coordinates for the Sheet Letters.

 

I will try tomorrow unless a straighter answer comes up. 

 

Code so far:

 

'Open the field reference data table up and write the 2 letter code from the IACS code

Close All Interactive

Open Table "D:\Data\FMapData\Agdat_fammap\ABSUGAR\ANalysis_2015_Season\V_fields.TAB" Interactive

Commit Table V_fields As "D:\Data\FMapData\Agdat_fammap\ABSUGAR\ANalysis_2015_Season\V_fields_Locations.TAB" TYPE NATIVE Charset "WindowsLatin1" Interactive

Open Table "D:\Data\FMapData\Agdat_fammap\ABSUGAR\ANalysis_2015_Season\V_fields_Locations.TAB" Interactive

Alter Table "V_fields_Locations" ( add Grid100KLetter Char(2),NA0595_X Integer,NA0595_Y Integer,Lat Float,Long Float ) Interactive

Update V_fields_Locations Set Grid100KLetter = Left$(IACS_of_Field,2)

Commit Table V_fields_Locations Interactive

 

Open Table "D:\Data\FMapData\Agdat_fammap\ABSUGAR\ANalysis_2015_Season\BNG_Grid_100K.TAB"

Browse * From V_fields_Locations

 

'Update the coordinates to be lower x and y of the 100K map sheet based on the first 2 leters of the code

'And then add the appropriate offsets to this point based on the numbers in the IACS code

Set CoordSys Earth Projection 8, 79, "m", -2, 49, 0.9996012717, 400000, -100000

Add Column "V_fields_Locations" (NA0595_X )From BNG_Grid_100K Set To ObjectGeography(Obj,OBJ_GEO_MINX)  Where Grid100KLetter  = Sheet100K

Update V_fields_Locations Set NA0595_X = NA0595_X + Val(Mid$(IACS_of_Field,3,2))*1000 + Val(Mid$(IACS_of_Field,7,2))*10

 

Add Column "V_fields_Locations" (NA0595_Y)From BNG_Grid_100K Set To ObjectGeography(Obj,OBJ_GEO_MINY)  Where Grid100KLetter  = Sheet100K

Update V_fields_Locations Set NA0595_Y = NA0595_Y + Val(Mid$(IACS_of_Field,5,2))*1000 + Val(Mid$(IACS_of_Field,9,2))*10

Commit Table V_fields_Locations Interactive

Browse * From V_fields_Locations

 

'Create the points from the x and y coodinates, and then create a polygon the size of the planted area

Create Map For V_fields_Locations CoordSys Earth Projection 8, 79, "m", -2, 49, 0.9996012717, 400000, -100000

Update V_fields_Locations Set Obj = CreatePoint(NA0595_X, NA0595_Y)

Update V_fields_Locations Set Obj = Buffer(Obj,12,SQR(10000*Drilled_Area/3.14159),"m")

Map From V_fields_Locations

Set Map Zoom Entire Layer 1

[Robert Crossley]

Seb,

 

Am I correct in understanding that I use the numbers in the following way?

 

If I have a code of NB, the min  x coordinate is 100,000 and min y is 900,000?

 

Rob

[Seb]

Hi Robert,

yes...

cheers