Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1090 1340 (N21W13 L= 79)   301    137        0 (1036/1309-1149/1358) -15.9 W08-W15 N18-N22
      2  470 1224 (N31E30 L= 36)  1753    344        3 ( 337/1177- 623/1265) -14.6 E42-E19 N25-N34
      3  397 1010 (N20E40 L= 27)  1159    302        0 ( 288/ 956- 522/1084) -35.2 E48-E31 N14-N26
      4  679  926 (N09E22 L= 45)  1074    346        0 ( 626/ 865- 753/1027)  77.5 E25-E16 N05-N14
      5  178  936 (N20E64 L=  3)  1371    213        2 ( 140/ 840- 246/1017) -47.5 E69-E57 N13-N25
      6 1810  891 (S25W49 L=116)   371     54        0 (1786/ 866-1835/ 917)   9.5 W46-W52 S26-S24
      7 1680  564 (S43W34 L=101)  1275    243        4 (1573/ 504-1778/ 624)  -5.0 W23-W46 S45-S40
      8  429  385 (S17E64 L=  3)  3175    392        7 ( 313/ 254- 536/ 497)  81.5 E71-E58 S28-S08
      9  908  273 (S37E34 L= 33)  1701    231        4 ( 821/ 198- 990/ 356)  10.8 E45-E23 S41-S32
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y