Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1485 1530 (N24W45 L=106)   621    132        0 (1423/1489-1536/1561) -54.5 W41-W46 N20-N27
      2 1124 1488 (N32W19 L= 80)   491     64        1 (1093/1453-1154/1518) -63.4 W16-W20 N29-N34
      3 1493 1357 (N14W37 L= 98)   409    110        0 (1477/1300-1505/1414)  50.2 W35-W40 N11-N17
      4 1695 1317 (N07W51 L=112)  1445    186        4 (1623/1259-1776/1380) -36.0 W46-W58 N02-N11
      5  409 1325 (N38E38 L= 23)   300     33        0 ( 386/1304- 434/1343) -31.0 E40-E35 N35-N39
      6  640 1225 (N27E20 L= 41)   241     60        3 ( 619/1195- 666/1250)  33.7 E21-E17 N25-N28
      7 1518 1091 (S01W31 L= 92)   654    185        0 (1502/1009-1534/1191)  73.3 W28-W32 S05-N05
      8  973 1120 (N14E01 L= 60)   644    279        0 ( 897/1052-1041/1174) -54.5 E05-W01 N08-N16
      9 1914 1048 (S15W67 L=128)  4017    243        3 (1884/ 936-1945/1152)  72.9 W61-W75 S22-S09
     10  199 1106 (N27E61 L=360)   509     93        1 ( 176/1069- 227/1148) -35.5 E64-E57 N24-N29
     11 1848 1035 (S14W57 L=118)   619     93        1 (1832/ 985-1869/1083)  71.6 W54-W60 S16-S10
     12  472 1021 (N18E35 L= 26)  1663    416        2 ( 322/ 964- 642/1067) -24.6 E47-E23 N11-N22
     13 1819  924 (S19W52 L=113)  1841    240        6 (1745/ 865-1888/ 986) -33.7 W44-W61 S23-S13
     14 1339  760 (S15W13 L= 74)   613    169        0 (1293/ 690-1379/ 831)  41.2 W08-W16 S18-S11
     15 1118  642 (S17E03 L= 58)   392    115        0 (1080/ 614-1172/ 684) -56.3 E04-W00 S19-S13
     16 1574  513 (S39W26 L= 88)  5901    813        4 (1209/ 311-1841/ 683)   5.4 E06-W58 S44-S34
     17  574  531 (S13E40 L= 21)   613    141        1 ( 522/ 477- 616/ 580) -74.1 E42-E37 S16-S09
     18 1410  212 (S64W10 L= 71)  1753    127        2 (1368/ 175-1460/ 252)   0.0 W02-W19 S70-S57
  

    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