Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  586 1732 (N30E50 L=286)  2484    373        8 ( 501/1664- 650/1811)  85.6 E57-E42 N23-N36
      2  564 1482 (N12E38 L=298)  1546    282        4 ( 488/1383- 651/1601)  76.0 E41-E36 N05-N21
      3  712 1462 (N13E28 L=309)   451    111        2 ( 683/1410- 741/1505) -80.5 E29-E26 N10-N16
      4 1531 1426 (N28W25 L=  1)   250     72        0 (1512/1388-1558/1457) -36.9 W23-W27 N26-N30
      5  879 1310 (N07E14 L=323)   333     99        1 ( 846/1270- 915/1349)  78.7 E15-E12 N04-N09
      6 1828 1226 (N24W53 L= 30)   541     73        0 (1804/1191-1854/1264) -31.0 W50-W56 N22-N26
      7 1528 1226 (N15W27 L=  3)   427    111        0 (1480/1182-1559/1256) -21.8 W23-W29 N13-N16
      8 1348 1215 (N11W15 L=352)   322     72        1 (1325/1173-1365/1253) -53.1 W13-W16 N08-N12
      9  192 1067 (S17E62 L=275)  2080    242        2 ( 156/ 962- 249/1150)  59.7 E66-E55 S23-S11
     10 1765 1035 (N10W48 L= 25)   387     69        0 (1750/ 989-1782/1068)  90.0 W46-W49 N07-N12
     11  560  815 (S27E26 L=310)  1598    541        4 ( 362/ 678- 668/1002)  46.2 E44-E20 S39-S14
     12 1331  943 (S05W19 L=356)   225     72        0 (1318/ 905-1344/ 981) -63.4 W17-W20 S07-S03
     13 1808  892 (N03W55 L= 32)   366     35        2 (1788/ 866-1827/ 917) -33.7 W53-W57 N02-N04
     14 1771  672 (S09W61 L= 37)   547     86        1 (1731/ 648-1806/ 699)  20.6 W55-W64 S10-S06
     15 1583  515 (S23W52 L= 28)  1626    290        2 (1525/ 392-1634/ 643) -63.4 W47-W56 S31-S15
     16 1777  530 (S16W72 L= 48)  1286    121        4 (1732/ 496-1815/ 568)  63.4 W67-W79 S19-S13
  

    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