Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1219 1652 (N45W20 L=213)   326     80        0 (1181/1636-1262/1671) -23.2 W16-W23 N43-N46
      2 1387 1546 (N36W31 L=224)   448     62        3 (1350/1518-1420/1570) -45.0 W28-W33 N34-N38
      3  380 1302 (N22E47 L=147)   477     96        0 ( 341/1270- 408/1326)   9.5 E50-E43 N19-N23
      4 1025 1262 (N19W01 L=194)   242    130        0 ( 989/1214-1066/1313) -60.3 E01-W03 N15-N22
      5  258 1192 (N15E57 L=137)  1258    191        2 ( 233/1104- 283/1271) -84.3 E60-E54 N09-N19
      6 1621 1030 (N02W39 L=233)   277     52        1 (1592/1009-1645/1054)  14.0 W37-W41 N00-N03
      7  842  877 (S05E12 L=182)   996    290        5 ( 729/ 783- 940/ 977) -45.0 E18-E05 S10-N01
      8 1516  867 (S08W31 L=225)   621    125        1 (1479/ 814-1556/ 910)  45.0 W28-W34 S10-S05
      9 1353  778 (S12W20 L=214)   561    184        2 (1296/ 717-1379/ 837)  54.5 W16-W22 S16-S08
     10  528  699 (S15E34 L=159)  1604    301        4 ( 419/ 631- 645/ 774) -13.2 E42-E25 S19-S10
     11  923  588 (S23E08 L=186)  1092    279        0 ( 844/ 485-1021/ 674)  45.0 E14-E01 S30-S17
     12 1512  587 (S26W34 L=227)   988    248        3 (1426/ 526-1628/ 643)  21.8 W27-W42 S30-S22
     13  665  588 (S23E26 L=168)   598    162        2 ( 616/ 526- 719/ 629) -54.5 E29-E22 S27-S20
     14  578  560 (S25E33 L=160)   313     58        1 ( 562/ 531- 596/ 590) -76.0 E34-E32 S26-S22
     15 1511  497 (S32W36 L=230)   780    154        0 (1456/ 432-1549/ 555)  60.9 W32-W38 S37-S28
     16  800  375 (S39E20 L=173)   471     37        1 ( 772/ 355- 829/ 396)  11.3 E22-E17 S40-S37
     17 1241  319 (S46W16 L=210)   987    111        2 (1195/ 284-1295/ 351) -12.5 W12-W21 S49-S42
     18  858  223 (S54E22 L=172)   576     63        2 ( 824/ 199- 885/ 249)   7.1 E25-E17 S57-S51
     19 1032  197 (S58E03 L=190)   469     42        1 (1011/ 177-1060/ 219) -31.0 E05-W00 S60-S55
  

    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