Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  454 1582 (N31E54 L=  0)   834    149        2 ( 436/1512- 479/1641) -45.0 E60-E50 N26-N35
      2  752 1473 (N27E24 L= 31)   873    229        1 ( 718/1389- 784/1576) -58.4 E29-E20 N21-N34
      3 1622 1375 (N29W42 L= 96)  3264    654        3 (1389/1171-1800/1527) -29.2 W22-W57 N16-N38
      4  669 1388 (N21E28 L= 27)   321     86        0 ( 627/1376- 724/1402)   8.1 E30-E23 N19-N21
      5 1668 1181 (N17W43 L= 98)   459    124        0 (1642/1114-1683/1238) -76.0 W41-W44 N12-N20
      6 1519 1166 (N14W31 L= 86)   426    140        0 (1483/1130-1561/1213) -39.8 W28-W34 N12-N17
      7  241  910 (S13E57 L=357)  1063    246        2 ( 161/ 817- 319/ 987) -27.9 E67-E48 S18-S09
      8  613  805 (S16E25 L= 30)   298     88        0 ( 579/ 782- 661/ 822)  26.6 E27-E21 S17-S14
      9  578  732 (S21E27 L= 27)   255     65        0 ( 564/ 694- 592/ 763)  78.7 E28-E26 S23-S18
     10 1595  642 (S17W44 L= 98)  1015    218        1 (1512/ 575-1671/ 718)  48.0 W39-W48 S22-S11
     11  674  643 (S26E20 L= 34)   265     64        0 ( 643/ 622- 703/ 663) -26.6 E22-E18 S26-S24
     12 1414  491 (S29W35 L= 89)   471    101        2 (1369/ 468-1464/ 519)  38.7 W31-W38 S31-S26
  

    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