topic Re: Polygon Concentration in Agriculture Questions

Polygon Concentration

Md_Abdullah_All_Sourav — Mon, 07 May 2018 16:58:52 GMT

I have yield data of a field. Almost quarter million of polygons with the size of 30ft by 1.2ft representing the yield of the region. I have been asked to concentrate (dissolve) every 5 polygons, thus the result will be 30 ft by 6ft (1.2*5). The new polygon will contain the average yield of the old 5 polygons. I need to contrate them from left to right. Is there any way to do that?

Thanks

Re: Polygon Concentration

DanPatterson_Retired — Mon, 07 May 2018 17:06:55 GMT

Are they in a sequential order within your attribute table? And if so, is the pattern left to right, left to right.

Are all 1/4 million fields in one table?

It would be quite simple to produce a sequential sum in groupings of 5 if that was the case.

There are several means, but clarification would help

Re: Polygon Concentration

Md_Abdullah_All_Sourav — Tue, 08 May 2018 16:12:15 GMT

Hi Dan,

They are in sequential order in the attribute table and their location is from left to right (I might be wrong though). Can you please tell me how can I do the sequential sum in the grouping of 5?

Re: Polygon Concentration

DanPatterson_Retired — Tue, 08 May 2018 16:25:54 GMT

The reason that I ask, is because each row may or may not be equally divisible by 5, so I wondered if there was a Row, Column indicator/key in the table to facilitate the groupings. If you had a sample of the table, even a hundred rows or so, it would be easier to show you

Re: Polygon Concentration

Md_Abdullah_All_Sourav — Tue, 08 May 2018 16:54:04 GMT

This is the problem we talked about in our meeting. There could be different polygons in each row. For the first pass of the harvester there could be 10000 points (later converted to polygons), may for the second pass, it could be 10003. They are inconsistency. I have attached the data here. I also provided details of the problem here.

Re: Polygon Concentration

JoshuaBixby — Tue, 08 May 2018 21:53:50 GMT

I tried downloading the data from Box, but it says the file isn't available to me, so it looks like you need to update your permissions.

Re: Polygon Concentration

DanPatterson_Retired — Tue, 08 May 2018 21:59:43 GMT

me too... not part of that University...

Just remove any sensitive fields and zip and attach, but not when you are responding to the message, you actually have to go back to your question directly

Re: Polygon Concentration

Md_Abdullah_All_Sourav — Wed, 09 May 2018 03:26:50 GMT

Can you please check now. It's here.

Re: Polygon Concentration

DanPatterson_Retired — Wed, 09 May 2018 03:34:52 GMT

Got it.... for others, make sure you delete the lock file before you try to use it

Re: Polygon Concentration

Md_Abdullah_All_Sourav — Wed, 09 May 2018 04:08:05 GMT

I have deleted the lock file from the folder.

Re: Polygon Concentration

DanPatterson_Retired — Wed, 09 May 2018 20:37:28 GMT

What field is it you want the summary on? I currently have it in array format for processing

a.dtype
dtype([('FID', '<i4'), ('Shape', '<f8', (2,)), ('Field', '<U254'), ('Dataset', '<U254'), ('Product', '<U254'), ('Obj__Id', '<f8'),

('Time', '<M8[us]'), ('Duration_s', '<f8'), ('Area_Count', '<U254'), ('Elevation_', '<f8'), ('Swth_Wdth_', '<f8'),

('Y_Offset_f', '<f8'), ('Hrv_Mass_W', '<f8'), ('Distance_f', '<f8'), ('Crop_Flw_M', '<f8'), ('Fuel_Con_T', '<f8'),

('Track_deg_', '<f8'), ('Moisture__', '<f8'), ('Pass_Num', '<f8'), ('Crop_Flw_V', '<f8'), ('Yld_Mass_W', '<f8'),

('Yld_Mass_D', '<f8'), ('Yld_Vol_We', '<f8'), ('Yld_Vol_Dr', '<f8'), ('One_Second', '<f8'), ('Speed_mph_', '<f8'),

('Prod_ac_h_', '<f8'), ('Date', '<M8[us]')])

a[0] first record

list(zip(nms, a[0]))

[('FID', 0),
('Shape', array([-93.69, 41.98])),
('Field', 'Been'),
('Dataset', '3/20/2018'),
('Product', 'NO Product'),
('Obj__Id', 1.0),
('Time', numpy.datetime64('2017-10-30T00:00:00.000000')),
('Duration_s', 0.20499999999999999),
('Area_Count', 'On'),
('Elevation_', 1040.4000000000001),
('Swth_Wdth_', 30.0),
('Y_Offset_f', 0.0),
('Hrv_Mass_W', 1.4137),
('Distance_f', 0.53149999999999997),
('Crop_Flw_M', 3.7258),
('Fuel_Con_T', 10.48),
('Track_deg_', 359.62400000000002),
('Moisture__', 16.899999999999999),
('Pass_Num', 1.0),
('Crop_Flw_V', 239.52000000000001),
('Yld_Mass_W', 2086.6444000000001),
('Yld_Mass_D', 1993.0999999999999),
('Yld_Vol_We', 37.261499999999998),
('Yld_Vol_Dr', 35.590000000000003),
('One_Second', 1.0),
('Speed_mph_', 1.7677),
('Prod_ac_h_', 6.4028),
('Date', numpy.datetime64('2017-10-30T00:00:00.000000'))]

and you want a 5 plot running sum?

Any further clarification you want to provide?

Re: Polygon Concentration

Md_Abdullah_All_Sourav — Fri, 11 May 2018 03:17:45 GMT

Thanks Dan. The objects of interest are Yld_Vol_Dr'(bu/ac) and Moisture content. Thus, I just need to merge every 5 polygons, make one finally which will contain the average of the Yld_Vol_Dr' and moisture of the previous 5 polygons.

Re: Polygon Concentration

DanPatterson_Retired — Sat, 12 May 2018 01:37:41 GMT

Some results, each pair represents the Yield and Moisture average in 5 unit blocks

The first block before the ... represents 8 blocks of 5. The last block is the last 8 blocks of 5

Check manually to make sure it is alright.

The code is simple, but I am still concerned about 'end-wrapping'

It seems that you have 40 blocks per row, so if the number of blocks is not divisible by 5,then you can't wrap the end or one rows values to the beginning of the next, so you have to return a truncated average.

Advise on any anomolies

Sample results

array([( 35.59, 16.9 ), ( 61.14, 16.9 ), (122.4 , 16.9 ), (188.36, 16.9 ),
(260.38, 16.9 ), (289.45, 16.08), (259.13, 16.08), (242.38, 16.08),

...,
       (208.47, 15.86), (230.03, 15.86), (221.9 , 15.86), (229.2 , 15.86),
       (224.42, 15.86), (233.92, 15.86), (206.63, 15.86), (218.36, 15.86)],
      dtype=[('Yield', '<f8'), ('Moisture', '<f8')])