Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  273 1434 (N10E61 L= 22)  1344    215        2 ( 236/1391- 315/1502) -37.9 E65-E54 N07-N14
      2  186 1382 (N06E70 L= 13)  1133    198        2 ( 159/1345- 221/1428)  90.0 E74-E65 N04-N09
      3  536 1318 (N05E35 L= 48)   322     91        0 ( 494/1301- 585/1332)  18.4 E37-E31 N03-N05
      4  213 1107 (S10E58 L= 25)   941    256        0 ( 137/1085- 298/1134)   9.5 E67-E49 S11-S08
      5 1498 1092 (N05W28 L=111)   244    125        0 (1443/1060-1550/1114)  39.8 W24-W30 N02-N07
      6  546 1016 (S13E29 L= 54)   377    138        0 ( 498/ 961- 599/1063)  54.5 E32-E26 S16-S09
      7  369  898 (S22E42 L= 41)   683    257        0 ( 276/ 864- 450/ 925)   7.6 E50-E35 S24-S20
      8 1283  874 (S11W18 L=101)   543    237        0 (1183/ 831-1401/ 905)  30.3 W12-W24 S14-S07
      9  703  833 (S22E17 L= 66)   337    126        0 ( 639/ 815- 763/ 865)  26.6 E21-E13 S23-S19
     10  917  798 (S21E03 L= 80)   239    112        0 ( 903/ 745- 930/ 865) -81.9 E04-E02 S24-S17
     11  174  776 (S31E61 L= 22)   541    150        2 ( 146/ 748- 198/ 808)  38.7 E64-E58 S32-S28
     12 1492  773 (S13W33 L=116)   444    145        2 (1426/ 748-1548/ 800)  14.0 W29-W37 S15-S11
     13  526  680 (S34E29 L= 54)   484    128        0 ( 484/ 629- 575/ 714)  40.6 E33-E25 S37-S31
     14  789  535 (S40E08 L= 75)   268     83        0 ( 742/ 523- 835/ 551)  23.2 E12-E04 S41-S38
     15 1620  396 (S32W59 L=142)   457    104        0 (1586/ 382-1654/ 409)  18.4 W55-W61 S33-S30
     16 1559  371 (S35W55 L=138)   385     51        0 (1539/ 352-1582/ 391) -11.3 W52-W57 S35-S33
     17 1287  311 (S46W38 L=121)   331     87        0 (1272/ 280-1301/ 344) -38.7 W35-W41 S47-S43
  

    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