Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1444 1781 (N55W57 L=165)  1095     96        1 (1415/1740-1483/1815)  12.1 W48-W76 N51-N57
      2  535 1683 (N49E50 L= 57)   523     73        0 ( 507/1651- 564/1714)  80.5 E51-E47 N46-N52
      3 1482 1579 (N38W41 L=148)   327     36        0 (1462/1557-1507/1597)  11.3 W38-W43 N36-N39
      4  750 1438 (N31E19 L= 88)  2241    540        2 ( 547/1299- 948/1560) -27.9 E39-E04 N21-N39
      5 1297 1268 (N18W19 L=126)   483    210        0 (1259/1189-1363/1332)  56.3 W15-W23 N13-N22
      6  994  997 (N03E02 L=105)   332     61        1 ( 960/ 969-1020/1025)  45.0 E04-W00 N00-N04
      7 1693  938 (S04W46 L=153)   339     82        2 (1660/ 913-1721/ 973) -45.0 W42-W48 S05-S01
      8  845  660 (S18E13 L= 95)  2002    648        4 ( 644/ 431- 948/ 893)  50.0 E31-E05 S34-S03
      9 1507  733 (S16W32 L=139)   326    106        0 (1477/ 686-1536/ 775)  59.0 W29-W33 S18-S13
     10 1598  529 (S31W44 L=151)  4391    706        6 (1334/ 362-1782/ 736)  37.1 W24-W57 S42-S17
     11  588  626 (S21E31 L= 76)   249     71        1 ( 556/ 596- 610/ 655)  45.0 E33-E29 S22-S18
     12 1144  445 (S35W07 L=115)  1610    376        2 (1010/ 354-1276/ 538)  33.7 E03-W16 S41-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