Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1206 1746 (N54W07 L= 84)   387     84        0 (1178/1726-1250/1763)  -7.1 W03-W11 N51-N55
      2  482 1511 (N26E47 L= 30)  1304    249        0 ( 437/1422- 519/1634) -65.0 E54-E44 N19-N34
      3 1347 1418 (N29W18 L= 96)  2211    530        1 (1140/1245-1466/1552) -35.3 W02-W27 N19-N37
      4  507 1403 (N19E40 L= 37)   367     59        1 ( 480/1372- 537/1428) -31.0 E43-E37 N17-N20
      5 1167 1227 (N15W07 L= 84)   368    145        1 (1137/1175-1212/1292) -54.5 W04-W10 N12-N19
      6 1349 1258 (N19W19 L= 96)   255    128        1 (1314/1201-1376/1288) -59.0 W16-W20 N15-N20
      7  465  721 (S23E36 L= 41)   258     48        0 ( 436/ 696- 494/ 744)  38.7 E38-E33 S25-S21
      8 1287  672 (S18W21 L= 98)   635    176        0 (1220/ 622-1362/ 728)  45.0 W16-W25 S21-S13
      9 1708  654 (S15W55 L=132)   580     46        0 (1683/ 630-1734/ 677) -15.9 W51-W58 S16-S13
     10 1644  582 (S21W51 L=129)   524     75        0 (1621/ 548-1674/ 625)  78.7 W49-W53 S22-S17
     11 1703  580 (S20W58 L=135)   738    127        0 (1673/ 537-1754/ 615)  68.2 W53-W62 S23-S17
     12  415  488 (S41E48 L= 30)   645     52        1 ( 392/ 463- 448/ 517) -26.6 E51-E43 S43-S38
  

    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