Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1190 1489 (N32W18 L=319)  1207    340        1 (1046/1433-1330/1533) -26.6 W08-W26 N26-N36
      2 1720 1473 (N22W62 L=  3)   493    105        1 (1673/1445-1764/1498) -35.5 W58-W66 N19-N24
      3  621 1292 (N26E24 L=277)   340    126        0 ( 563/1280- 684/1304)  -6.3 E28-E19 N25-N27
      4  415 1146 (N19E40 L=261)   504     54        1 ( 392/1112- 439/1179)  63.4 E42-E38 N16-N20
      5 1538 1002 (S02W32 L=333)   268     76        1 (1501/ 970-1563/1032) -53.1 W30-W33 S04-N00
      6 1924  972 (S12W72 L= 14)  1162    107        1 (1905/ 925-1949/1023) -66.8 W68-W81 S15-S08
      7 1793  803 (S19W54 L=355)  3956    523        4 (1702/ 613-1883/ 991)  52.7 W46-W65 S30-S08
      8 1405  841 (S10W22 L=323)   703    186        2 (1353/ 790-1449/ 898)  81.9 W18-W24 S13-S06
      9 1557  861 (S11W32 L=333)   287     69        0 (1520/ 834-1587/ 883) -38.7 W29-W34 S13-S09
     10  188  753 (S05E69 L=232)  3599    320        9 ( 143/ 603- 227/ 878)  63.4 E82-E61 S15-N03
     11 1482  713 (S19W26 L=327)   582    206        0 (1414/ 670-1572/ 762)  14.0 W21-W33 S21-S16
     12  886  698 (S13E12 L=289)   276     86        0 ( 841/ 683- 926/ 716) -18.4 E15-E09 S13-S11
     13  834  540 (S23E19 L=282)  1465    373        4 ( 767/ 454- 955/ 678) -63.4 E22-E11 S29-S13
     14  635  442 (S28E36 L=265)   729    140        1 ( 606/ 381- 667/ 499)  52.1 E41-E31 S32-S23
     15 1609  390 (S48W46 L=348)  1165    103        2 (1562/ 366-1651/ 427)  16.7 W40-W53 S50-S44
     16 1080  316 (S43E06 L=295)   638     59        2 (1046/ 287-1116/ 345) -11.3 E08-E03 S45-S40
  

    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