Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1859 1279 (S02W67 L=178)   548    143        0 (1821/1241-1892/1309) -35.5 W63-W70 S05-N00
      2 1592 1277 (N06W41 L=152)   267     47        1 (1572/1254-1608/1304)  33.7 W39-W43 N05-N07
      3  214 1207 (N34E56 L= 55)   458     69        1 ( 192/1183- 241/1231) -36.9 E58-E53 N32-N35
      4 1787  741 (S31W46 L=157)  2397    381        3 (1676/ 617-1896/ 878)  23.5 W34-W59 S36-S25
      5 1337  707 (S19W11 L=121)   216     81        0 (1294/ 691-1377/ 722) -31.0 W08-W13 S20-S17
      6  930  633 (S13E15 L= 95)   451    137        3 ( 886/ 588- 976/ 685)  23.2 E18-E11 S15-S10
      7 1584  637 (S31W27 L=137)   448     76        1 (1556/ 599-1604/ 679)  59.0 W24-W28 S33-S28
      8 1108  558 (S22E07 L=104)   275     99        0 (1075/ 523-1150/ 602) -74.1 E07-E04 S25-S18
      9  876  499 (S20E22 L= 88)   296     92        3 ( 844/ 470- 907/ 525) -76.0 E23-E20 S22-S18
  

    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