Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  544 1606 (N15E48 L=281)   358     76        1 ( 526/1570- 561/1647) -59.0 E49-E46 N12-N17
      2 1282 1545 (N31E00 L=329)   344     90        0 (1235/1510-1328/1566)   0.0 E03-W03 N29-N32
      3 1690 1516 (N46W34 L=  3)   715    112        1 (1667/1464-1707/1556)  76.0 W31-W35 N41-N49
      4 1516 1436 (N32W19 L=347)   325     79        1 (1495/1404-1541/1480) -45.0 W16-W20 N30-N34
      5  495 1422 (N02E43 L=286)   432     94        1 ( 477/1380- 523/1462) -53.1 E45-E40 N00-N04
      6  703 1382 (N05E28 L=301)   237     78        0 ( 666/1347- 728/1408) -18.4 E31-E25 N03-N06
      7  250 1209 (S15E57 L=272)  2172    462        3 ( 117/1170- 368/1240)  10.9 E71-E45 S18-S12
      8 1817  990 (N16W52 L= 21)   565     95        1 (1796/ 937-1833/1036) -51.3 W49-W54 N12-N18
      9  815  868 (S21E08 L=321)   192     70        0 ( 804/ 833- 832/ 910) -53.1 E09-E06 S22-S18
     10 1749  841 (N05W51 L= 20)   761    196        2 (1692/ 799-1806/ 890)   4.4 W44-W57 N03-N06
     11 1340  798 (S10W24 L=353)   215     99        0 (1321/ 744-1361/ 849) -74.1 W23-W25 S13-S06
     12 1891  799 (N09W68 L= 37)  1147    270        3 (1849/ 765-1941/ 843)  38.7 W62-W75 N05-N13
     13 1499  694 (S10W38 L=  7)   292     69        2 (1468/ 671-1526/ 713)  45.0 W35-W39 S12-S08
     14  986  575 (S33W12 L=341)  1039    248        1 ( 850/ 539-1090/ 632)  47.5 W03-W16 S38-S26
     15 1587  595 (S13W48 L= 17)   339     71        1 (1549/ 575-1624/ 620)  45.0 W46-W50 S14-S10
     16 1429  418 (S27W47 L= 16)  1761    388        3 (1303/ 341-1529/ 495)   0.0 W34-W58 S28-S24
  

    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