Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  694 1591 (N20E34 L= 73)   467    182        0 ( 655/1524- 731/1658) -38.7 E38-E28 N16-N24
      2  132 1259 (S06E75 L= 32)  1123    134        0 ( 107/1181- 156/1319) -77.5 E78-E72 S11-S02
      3  556 1228 (S04E32 L= 75)  1030    232        2 ( 464/1128- 641/1319)  56.3 E37-E28 S10-N02
      4  283 1160 (S11E52 L= 55)   822    126        2 ( 254/1106- 323/1209)  63.4 E54-E48 S14-S08
      5 1668 1080 (N11W38 L=145)  1011    291        2 (1636/ 979-1689/1186) -57.5 W34-W41 N05-N16
      6  178 1039 (S20E60 L= 46)   710    128        1 ( 157/ 984- 206/1097) -50.2 E64-E56 S22-S16
      7 1341  550 (S27W32 L=138)   252     43        0 (1313/ 535-1374/ 565)  45.0 W29-W33 S28-S25
  

    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