Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  714 1830 (N35E56 L=262)  1779    169        2 ( 680/1762- 742/1894) -32.3 E71-E48 N28-N40
      2 1650 1537 (N47W31 L=349)   838     75        0 (1624/1492-1674/1579) -74.1 W27-W33 N43-N50
      3 1084 1401 (N17E07 L=311)   444    172        2 (1049/1322-1110/1480) -48.4 E11-E03 N12-N21
      4 1563 1052 (N12W30 L=349)   400     92        0 (1528/1011-1601/1088)  63.4 W28-W32 N08-N14
      5  809  999 (S13E11 L=307)   607    266        0 ( 714/ 942- 929/1034)  12.1 E16-E02 S16-S11
      6 1844  932 (N15W58 L= 16)   887    201        3 (1797/ 859-1893/1023) -20.6 W50-W67 N12-N18
      7  662  849 (S25E17 L=302)   267     67        1 ( 639/ 818- 693/ 881) -14.0 E18-E14 S26-S23
      8  370  796 (S37E39 L=279)   450     97        0 ( 324/ 770- 409/ 817)  24.0 E44-E35 S39-S34
      9 1670  751 (S01W48 L=  6)   676    109        0 (1639/ 728-1722/ 792)   0.0 W44-W53 S02-N00
     10  780  638 (S35E03 L=315)   815    174        0 ( 695/ 619- 870/ 659)  20.0 E09-W02 S37-S31
     11 1434  360 (S29W54 L= 12)  3141    476        8 (1250/ 301-1670/ 432)  45.0 W43-W70 S39-S16
  

    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