Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  686 1472 (N30E26 L= 15)  1107    300        2 ( 616/1376- 742/1565) -38.2 E34-E20 N23-N35
      2 1199 1327 (N22W10 L= 52)  1582    360        3 (1029/1226-1374/1435) -28.6 E01-W22 N16-N28
      3  304 1350 (N19E55 L=346)   888    173        0 ( 282/1275- 326/1425) -73.3 E57-E51 N14-N24
      4  867 1300 (N19E11 L= 30)   757    238        0 ( 820/1210- 948/1403) -56.3 E14-E05 N13-N25
      5  722 1270 (N17E21 L= 21)   261     59        2 ( 693/1242- 748/1299) -26.6 E22-E18 N15-N18
      6  486 1131 (N07E36 L=  6)   834    187        1 ( 454/1047- 542/1211) -56.3 E38-E31 N02-N11
      7  397  977 (S03E42 L=359)   388    104        1 ( 360/ 956- 441/1006)  -8.1 E45-E38 S04-S01
      8 1002  501 (S31W01 L= 42)   301     76        1 ( 975/ 473-1028/ 531) -63.4 E01-W02 S33-S29
      9 1189  261 (S51W20 L= 62)   598    137        0 (1137/ 227-1225/ 289) -26.6 W14-W26 S54-S48
  

    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