Rotating fishnet made of minimum bounding boxes?

Question

Given a grid of points, I had to create something similar to a fishnet, and did so using the buffer -> minimum bounding geometry option as listed here (option A).

Since my grid comes from Climate Models, it is rotated (alpha=15.16 degrees counterclockwise). The fishnet obtained has no such rotation, however, and I need to rotate it by the same alpha angle.

My fishnet has 170090 individual objects (466 columns and 365 rows), and this is because each square had to inherit an order field from the point it contains, so I have not dissolved the squares into a unique object.

How can I rotate all my squares by alpha?

Each square has to rotate relative to its center.

Image:

Answer 1

Script below will create a new polygon feature class with features rotated counterclockwise around centroid.

import arcpy,math
arcpy.env.overwriteOutput=1
original_fc=r'C:\TEST.gdb\unrotatedPolygons'
copied_fc=r'C:\TEST.gdb\rotatedPolygons'
arcpy.CopyFeatures_management(in_features=original_fc, out_feature_class=copied_fc)
rotation=15.16 #counterclockwise rotation in degrees
with arcpy.da.UpdateCursor(copied_fc,'SHAPE@') as cursor:
    for row in cursor:
        array1=row[0].getPart()
        centroid=(row[0].centroid.X,row[0].centroid.Y)
        for vertice in range(row[0].pointCount):
            pnt=array1.getObject(0).getObject(vertice)
            pnt.X=math.cos(math.radians(rotation))*(pnt.X-centroid[0])-math.sin(math.radians(rotation))*(pnt.Y-centroid[1])+centroid[0]
            pnt.Y=math.sin(math.radians(rotation))*(pnt.X-centroid[0])+math.cos(math.radians(rotation))*(pnt.Y-centroid[1])+centroid[1]
        row[0]=arcpy.Polygon(array1)
        cursor.updateRow(row)

Math explained: How can I rotate a coordinate around a circle?