Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1498 1679 (N32W54 L=330)   520     81        0 (1469/1660-1528/1696)  -9.5 W51-W57 N30-N33
      2 1661 1410 (N13W53 L=329)   447     94        0 (1625/1384-1702/1434) -38.7 W50-W55 N10-N14
      3 1712 1271 (N04W51 L=327)   538     84        1 (1675/1249-1742/1289) -31.0 W48-W54 N01-N05
      4  418 1097 (N23E38 L=238)  2428    595        4 ( 235/ 987- 634/1226) -35.0 E54-E24 N12-N33
      5  175  962 (N18E64 L=212)   625    131        0 ( 136/ 940- 202/ 977)  -7.1 E68-E60 N16-N18
      6  342  892 (N12E48 L=228)   380     86        2 ( 315/ 866- 367/ 921) -56.3 E49-E45 N10-N13
      7  563  323 (S28E52 L=224)   500     88        1 ( 546/ 292- 581/ 345)  38.7 E54-E49 S30-S26
      8  638  280 (S32E49 L=227)   656    119        1 ( 613/ 246- 659/ 316)  56.3 E52-E45 S35-S29
  

    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