Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1538 1687 (N47W35 L=260)  1693    198        1 (1463/1628-1611/1727) -17.1 W26-W42 N43-N50
      2  688 1576 (N25E30 L=195)   596    107        1 ( 664/1519- 708/1631) -49.4 E33-E27 N20-N28
      3  420 1522 (N19E50 L=176)   823    127        1 ( 393/1466- 445/1582) -63.4 E51-E46 N15-N23
      4 1391 1397 (N22W19 L=244)  1506    395        1 (1226/1273-1527/1481) -23.2 W07-W28 N15-N25
      5  749 1424 (N15E22 L=203)   294     60        1 ( 715/1398- 780/1446) -14.0 E24-E20 N13-N16
      6 1770 1363 (N26W51 L=276)   343     58        1 (1753/1330-1788/1394) -63.4 W49-W52 N24-N28
      7 1800 1302 (N23W53 L=278)   536     53        1 (1775/1279-1828/1324)  26.6 W50-W56 N21-N24
      8 1879 1258 (N22W63 L=288)   671     80        1 (1859/1228-1901/1298) -76.0 W60-W67 N19-N24
      9  209 1253 (N01E62 L=163)   726     89        3 ( 174/1231- 249/1285)  -5.7 E66-E56 N00-N03
     10  259 1021 (S13E53 L=172)   582    126        3 ( 239/ 956- 276/1078) -74.1 E55-E51 S16-S09
     11  517 1022 (S11E32 L=193)   301     98        0 ( 491/ 974- 551/1054)  53.1 E33-E29 S13-S08
     12  329  838 (S23E47 L=178)   521     89        0 ( 280/ 819- 379/ 856)   0.0 E52-E42 S24-S22
     13  887  835 (S18E06 L=219)   244    104        0 ( 842/ 815- 939/ 849)  18.4 E08-E02 S20-S17
     14 1733  807 (S08W51 L=276)   291     59        2 (1708/ 782-1756/ 831)  45.0 W49-W53 S09-S06
     15  672  677 (S31E20 L=205)   605    121        1 ( 632/ 626- 719/ 726)  53.1 E22-E16 S34-S26
     16 1524  538 (S28W43 L=269)  2685    480        2 (1329/ 425-1710/ 671)  46.4 W32-W54 S37-S16
     17 1448  362 (S40W48 L=273)   441    127        0 (1411/ 325-1499/ 416)  90.0 W46-W49 S43-S35
  

    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