Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  416 1571 (N10E57 L=276)   770    155        3 ( 389/1516- 448/1624) -74.1 E59-E54 N06-N13
      2  647 1508 (N11E37 L=297)   307     69        0 ( 611/1491- 683/1523)  11.3 E38-E33 N10-N12
      3  224 1288 (S10E63 L=271)   612     95        0 ( 201/1239- 255/1332)  90.0 E64-E60 S13-S07
      4  183 1103 (S22E61 L=273)   610    121        0 ( 162/1044- 215/1150)  66.8 E62-E58 S25-S18
      5 1865 1026 (N20W56 L= 30)   393     97        0 (1838/1007-1889/1040)   9.5 W53-W59 N18-N21
      6  499  955 (S24E30 L=304)   728    143        1 ( 439/ 915- 562/1002) -11.3 E35-E25 S24-S22
      7 1142  483 (S33W26 L=359)   398     95        0 (1113/ 458-1184/ 519)  -9.5 W22-W28 S33-S31
      8 1611  467 (S18W58 L= 31)   405     96        0 (1575/ 450-1653/ 489)  45.0 W54-W60 S19-S14
      9 1344  397 (S31W43 L= 17)   901    181        2 (1271/ 366-1413/ 439)   9.5 W35-W48 S32-S28
  

    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