Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1365 1278 (N13W27 L=  7)  1636    517        3 (1247/1157-1530/1475) -76.0 W23-W34 N03-N27
      2  705 1288 (N30E14 L=326)  1082    385        4 ( 570/1222- 865/1379) -41.6 E23-E04 N22-N38
      3 1814 1232 (S02W60 L= 40)   462     31        1 (1793/1212-1836/1252)  -9.5 W57-W63 S03-N00
      4  331  986 (N18E48 L=293)  1511    408        2 ( 193/ 949- 512/1039)  -3.6 E63-E31 N15-N21
      5 1835  975 (S17W55 L= 36)   439    142        0 (1804/ 910-1866/1038)  29.1 W51-W60 S19-S14
      6  650  905 (N08E25 L=315)   487    120        1 ( 602/ 875- 687/ 951)  26.6 E28-E21 N06-N10
      7  589  924 (N10E28 L=312)   263     50        2 ( 564/ 902- 610/ 950) -33.7 E30-E27 N08-N11
      8 1207  460 (S32E02 L=338)   428     93        0 (1183/ 417-1227/ 510)  56.3 E04-W00 S34-S28
      9  465  405 (S20E57 L=283)   529    116        1 ( 427/ 352- 511/ 445)  82.9 E59-E55 S24-S16
  

    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