Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  582 1642 (N38E45 L=343)  1217    196        2 ( 497/1614- 671/1675)   3.2 E53-E34 N36-N40
      2  743 1501 (N30E25 L=  3)   262     59        2 ( 726/1479- 767/1524)  45.0 E26-E23 N28-N31
      3 1107 1397 (N26W02 L= 30)  1017    329        1 (1027/1317-1162/1509) -59.7 E02-W06 N20-N33
      4 1750 1456 (N35W61 L= 88)  1755    260        2 (1634/1404-1846/1506)  -7.6 W46-W76 N31-N38
      5  336 1386 (N18E55 L=332)   409    146        0 ( 304/1320- 359/1425) -31.0 E60-E50 N14-N20
      6 1827 1147 (N15W60 L= 87)   595    129        0 (1797/1112-1873/1196) -29.1 W57-W66 N12-N17
      7  653 1135 (N06E24 L=  3)   275    141        0 ( 603/1089- 707/1188)  54.5 E27-E21 N02-N09
      8  514  890 (S11E32 L=355)   233     53        1 ( 487/ 871- 545/ 906) -11.3 E34-E29 S11-S09
      9 1069  739 (S15W05 L= 32)   427    165        1 (1013/ 692-1145/ 781)  40.6 W01-W09 S18-S12
     10 1550  593 (S21W42 L= 69)   345    142        0 (1499/ 551-1599/ 642)  45.0 W38-W44 S24-S17
     11 1690  481 (S28W63 L= 90)   663     98        1 (1644/ 450-1734/ 512)  26.6 W58-W67 S30-S25
  

    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