Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  944 1866 (N62W39 L= 68)   971    153        1 ( 914/1811- 985/1922)  14.0 W28-W56 N56-N66
      2  758 1756 (N61W09 L= 37)   505    116        1 ( 722/1714- 807/1788)  12.1 W02-W16 N58-N62
      3 1550 1607 (N22W55 L= 83)  1694    300        6 (1407/1569-1654/1648) -38.7 W45-W63 N17-N29
      4 1230 1317 (N17W20 L= 48)   570    262        2 (1151/1268-1292/1365) -20.6 W15-W23 N14-N20
      5 1323 1248 (N11W23 L= 51)   340    101        2 (1280/1198-1356/1286)  15.9 W18-W25 N09-N12
      6 1826 1265 (S04W60 L= 89)   453      7        0 (1805/1245-1846/1286) -56.3 W57-W62 S05-S02
      7 1041 1240 (N18W07 L= 35)   419    157        1 ( 960/1224-1113/1254) -32.0 W02-W10 N15-N20
      8 1372 1201 (N07W24 L= 53)   273     80        2 (1346/1166-1403/1229)   0.0 W22-W26 N05-N08
      9 1342  990 (S04W17 L= 45)   540    122        3 (1296/ 929-1385/1032)  29.7 W12-W20 S06-S01
     10  329  909 (N16E48 L=341)   360    138        0 ( 311/ 843- 342/ 971)  81.9 E49-E46 N12-N19
     11 1104  826 (S06E01 L= 28)   895    338        2 (1049/ 686-1192/ 918)  50.7 E04-W06 S14-S02
     12  628  543 (S11E37 L=352)   364     78        1 ( 614/ 503- 648/ 589)  59.0 E38-E34 S13-S08
     13 1660  556 (S41W31 L= 60)   734     81        3 (1620/ 526-1698/ 582)  -7.1 W27-W36 S44-S38
     14  902  518 (S19E21 L=  7)   363    148        0 ( 840/ 497- 981/ 540) -35.5 E24-E16 S22-S17
     15  579  367 (S21E50 L=338)  1549    174        3 ( 550/ 282- 613/ 449)  38.2 E59-E43 S26-S15
     16 1262  310 (S45E08 L= 21)  1909    433        3 (1107/ 219-1437/ 417)   3.7 E23-W09 S49-S41
  

    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