Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1288 1683 (N45W34 L=127)   293     32        0 (1270/1656-1303/1704)  71.6 W32-W36 N43-N46
      2 1546 1468 (N26W45 L=137)   379    121        0 (1502/1435-1605/1498) -39.8 W42-W48 N23-N28
      3 1787 1421 (N19W68 L=160)  1436    137        1 (1749/1364-1816/1486)  20.6 W61-W78 N15-N21
      4  554 1301 (N27E30 L= 62)   595    220        0 ( 463/1250- 612/1330)  -9.5 E37-E25 N23-N29
      5 1447 1299 (N17W31 L=124)   290    115        1 (1400/1270-1496/1329) -39.8 W28-W34 N14-N19
      6 1270 1164 (N11W17 L=109)   201     54        0 (1240/1146-1291/1184)   0.0 W14-W18 N10-N12
      7  176 1058 (N13E66 L= 26)  1115    120        1 ( 140/1003- 195/1114)  74.1 E72-E63 N09-N16
      8  987  872 (S03E04 L= 88)  1651    335        4 ( 934/ 724-1037/1035) -87.0 E05-E01 S12-N07
      9 1753  853 (S15W50 L=142)   498     19        0 (1729/ 831-1776/ 879) -36.9 W47-W52 S16-S12
     10  192  813 (S02E67 L= 25)   508     70        0 ( 173/ 783- 213/ 841) -36.9 E69-E64 S03-N00
     11  293  570 (S18E64 L= 28)   510     64        0 ( 263/ 551- 324/ 584)   5.7 E69-E59 S19-S16
     12 1524  524 (S33W32 L=124)   318     76        0 (1483/ 513-1566/ 539)   0.0 W28-W35 S34-S32
  

    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