Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  377 1646 (N19E64 L= 29)  1162    104        4 ( 323/1610- 418/1687)  51.3 E69-E59 N16-N22
      2  463 1622 (N18E53 L= 40)   391     72        1 ( 433/1590- 494/1650)  63.4 E54-E51 N15-N20
      3 1719 1451 (N36W40 L=134)   946    164        0 (1650/1413-1788/1483)   3.8 W33-W48 N33-N38
      4 1415 1446 (N26W15 L=108)   276     85        0 (1383/1411-1451/1479)  63.4 W13-W16 N23-N29
      5  328 1227 (S08E49 L= 44)   266     35        1 ( 307/1208- 350/1250)  63.4 E50-E47 S08-S06
      6  156 1127 (S15E66 L= 28)  4637    415        9 ( 113/ 951- 226/1239)  70.6 E74-E58 S25-S08
      7 1246 1183 (N06W09 L=103)   213     61        0 (1212/1163-1275/1199)  14.0 W06-W10 N05-N07
      8  470 1105 (S13E36 L= 58)   442    135        2 ( 405/1064- 524/1129)  -6.3 E40-E31 S14-S11
      9 1900 1049 (N16W61 L=155)  1062    171        2 (1848/1001-1947/1088)  -9.5 W55-W70 N13-N18
     10 1081  899 (S13W06 L= 99)  1042    535        0 ( 863/ 822-1291/ 979)  34.3 E05-W17 S21-S05
     11  483  927 (S23E32 L= 61)   441    122        1 ( 413/ 905- 542/ 954)  24.0 E37-E27 S24-S20
     12 1207  645 (S24W20 L=114)   329     84        0 (1164/ 630-1256/ 663)   9.5 W17-W23 S25-S22
     13 1341  526 (S28W33 L=127)   379     51        1 (1303/ 503-1368/ 549)  11.3 W30-W35 S28-S25
  

    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