Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1272 1727 (N41W45 L=239)   354     55        0 (1240/1712-1303/1741) -21.8 W42-W48 N39-N41
      2  975 1645 (N46W18 L=211)   248     47        1 ( 949/1628-1006/1659) -21.8 W15-W20 N44-N47
      3 1535 1576 (N23W53 L=246)   502    108        2 (1495/1534-1573/1606) -63.4 W50-W54 N20-N26
      4 1468 1449 (N19W40 L=234)  1582    419        0 (1292/1340-1622/1529) -53.6 W31-W46 N08-N27
      5 1843 1078 (S13W58 L=251)   429     46        2 (1817/1052-1867/1105)  56.3 W54-W60 S14-S11
      6 1889  910 (S25W62 L=255)   996     89        1 (1856/ 870-1920/ 942)  12.5 W56-W68 S27-S23
      7 1262  845 (S09W09 L=203)   267    136        2 (1244/ 796-1287/ 919)  74.1 W08-W10 S12-S04
      8 1620  522 (S42W30 L=223)  1513    206        1 (1528/ 457-1688/ 595)  17.4 W20-W37 S44-S38
      9  424  440 (S15E60 L=134)   522    138        0 ( 403/ 374- 451/ 496)  56.3 E63-E56 S19-S10
  

    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