Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1276 1728 (N39W41 L=202)  1111    237        1 (1169/1677-1383/1760)  -2.5 W29-W52 N36-N40
      2 1433 1539 (N23W40 L=201)   259     59        0 (1403/1516-1467/1562) -63.4 W38-W41 N20-N24
      3  493 1339 (N33E29 L=132)   778    189        1 ( 407/1319- 592/1357) -10.6 E36-E20 N31-N35
      4 1109 1271 (N15W11 L=171)   289     86        0 (1060/1253-1153/1286) -31.0 W07-W12 N13-N16
      5  695 1234 (N22E15 L=145)   310     81        0 ( 658/1213- 735/1251)   0.0 E17-E12 N21-N23
      6 1812 1218 (S05W58 L=218)   344     58        0 (1786/1188-1829/1246) -63.4 W55-W59 S06-S02
      7 1743 1054 (S12W47 L=207)  1374    286        2 (1674/ 975-1805/1136) -59.7 W42-W53 S17-S05
      8  274 1004 (N18E51 L=109)  1698    303        2 ( 138/ 942- 372/1094)   4.1 E68-E40 N14-N21
      9 1842 1023 (S17W56 L=217)   499     54        0 (1822/ 996-1862/1059)  36.9 W53-W59 S18-S14
     10 1439  947 (S11W22 L=183)   257     69        1 (1423/ 910-1455/ 982)  45.0 W20-W24 S12-S08
     11 1568  800 (S22W29 L=189)   316     49        0 (1547/ 770-1583/ 835)  63.4 W26-W30 S24-S20
     12 1395  730 (S22W16 L=176)   720    184        0 (1330/ 649-1443/ 807)  37.9 W09-W19 S25-S18
     13 1185  645 (S22W01 L=161)  1012    252        2 (1144/ 545-1250/ 777)  50.2 E04-W07 S28-S15
     14 1419  599 (S31W15 L=176)   491    129        0 (1363/ 572-1470/ 643)   6.3 W10-W19 S32-S29
     15  363  585 (S08E55 L=105)   330     23        0 ( 348/ 561- 380/ 608)  33.7 E57-E53 S09-S06
     16  841  493 (S24E25 L=136)   571    165        2 ( 807/ 412- 857/ 575)  48.0 E29-E20 S28-S18
     17 1635  451 (S50W36 L=196)  1088    121        1 (1568/ 427-1695/ 482)  -3.4 W27-W44 S52-S47
  

    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