Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  611 1735 (N26E50 L=192)   750    174        0 ( 558/1658- 640/1796) -61.4 E56-E47 N19-N30
      2  474 1722 (N22E62 L=180)   823    102        2 ( 454/1672- 499/1768) -45.0 E68-E56 N18-N26
      3  380 1565 (N10E59 L=184)   645     62        2 ( 348/1532- 414/1601)  36.9 E60-E55 N07-N12
      4 1148 1453 (N20E04 L=239)   231     75        0 (1119/1419-1181/1481)  68.2 E05-E02 N17-N22
      5  204 1308 (S08E65 L=178)  2222    291        6 ( 173/1182- 260/1428)  86.2 E71-E59 S15-N00
      6  264 1326 (S06E59 L=184)   984    139        2 ( 231/1275- 295/1397)  82.9 E61-E56 S09-S01
      7  343 1269 (S08E49 L=193)   587    139        1 ( 281/1230- 408/1296)  16.7 E54-E44 S10-S05
      8  649 1113 (S11E24 L=219)   332    100        0 ( 624/1062- 674/1167) -45.0 E26-E21 S13-S08
      9 1120 1082 (S02W04 L=246)   188    115        0 (1069/1047-1163/1117)  -9.5 W00-W06 S02-N00
     10 1351  586 (S22W32 L=275)   226    109        0 (1319/ 539-1386/ 625) -14.0 W28-W36 S23-S20
     11 1628  535 (S16W54 L=297)   753    107        2 (1576/ 500-1674/ 573)  33.7 W51-W58 S19-S14
     12  951  388 (S46W17 L=259)  1961    433        4 ( 781/ 296-1136/ 474)   1.6 W00-W36 S49-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