Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  621 1722 (N63E09 L=298)   458     71        0 ( 592/1689- 650/1754)  24.0 E13-E03 N60-N65
      2 1477 1682 (N28W55 L=  2)   395    121        0 (1415/1663-1525/1708)  -8.7 W48-W61 N27-N29
      3 1118 1677 (N41W30 L=337)   415     45        0 (1091/1651-1145/1700)   9.5 W27-W33 N39-N41
      4  969 1620 (N42W17 L=324)   270     70        0 ( 944/1591-1004/1645)  18.4 W14-W20 N40-N42
      5  316 1574 (N60E59 L=248)   788     93        2 ( 266/1549- 359/1596)   2.6 E72-E49 N58-N61
      6 1346 1516 (N24W35 L=342)   382     78        1 (1303/1487-1388/1533)   0.0 W31-W37 N22-N24
      7 1419 1444 (N18W36 L=343)   415    137        0 (1384/1388-1441/1505)  45.0 W31-W38 N15-N21
      8 1871 1174 (S11W62 L=  9)  1365    190        3 (1832/1100-1919/1258)  23.2 W55-W70 S13-S07
      9  330 1122 (N28E44 L=263)  1046    229        2 ( 242/1053- 433/1179) -31.4 E53-E35 N22-N33
     10 1643  844 (S21W33 L=340)  1595    306        4 (1541/ 719-1713/ 968)  29.4 W23-W40 S26-S15
     11  196  832 (N14E63 L=244)   541    110        1 ( 153/ 784- 223/ 871) -50.2 E68-E59 N10-N16
     12 1150  783 (S10W01 L=308)   640    155        2 (1122/ 708-1174/ 860)  54.5 E02-W03 S14-S06
     13 1103  548 (S23E09 L=298)   304     91        0 (1085/ 516-1133/ 598)  80.5 E09-E07 S25-S19
     14 1267  483 (S32E01 L=306)   894    169        0 (1196/ 449-1355/ 532) -10.3 E05-W06 S33-S29
     15 1028  329 (S35E22 L=285)  1646    273        2 ( 908/ 276-1152/ 396)  -2.9 E32-E11 S37-S32
     16  577  335 (S23E53 L=254)   469     92        1 ( 542/ 302- 613/ 359) -68.2 E55-E50 S26-S21
     17  753  157 (S41E59 L=248)   942     80        2 ( 721/ 131- 776/ 186)  13.2 E69-E52 S43-S39
  

    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