Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  667 1704 (N61E06 L=245)   479     92        2 ( 643/1663- 702/1741)  24.0 E09-W00 N57-N63
      2 1285 1650 (N35W40 L=290)  1468    268        4 (1175/1592-1396/1694) -53.1 W33-W45 N28-N40
      3  452 1623 (N62E38 L=212)   405     44        1 ( 427/1601- 473/1644)   0.0 E43-E34 N59-N63
      4  983 1584 (N41W16 L=266)   242     61        1 ( 963/1554-1003/1618)  68.2 W13-W17 N38-N43
      5 1014 1411 (N29W10 L=261)   224     82        0 ( 979/1396-1055/1420) -36.9 W08-W13 N27-N30
      6  329 1157 (N31E46 L=205)   878    160        4 ( 252/1131- 414/1182) -19.4 E54-E37 N27-N33
      7 1785  861 (S25W47 L=297)  3067    475        7 (1652/ 699-1903/1022)  19.0 W32-W64 S32-S18
      8 1716  649 (S36W39 L=289)   344     41        1 (1696/ 627-1731/ 676)  45.0 W36-W40 S38-S34
      9  935  592 (S15E17 L=234)   238     66        0 ( 900/ 575- 968/ 617) -36.9 E18-E14 S17-S13
     10  909  442 (S24E24 L=227)   760    278        0 ( 852/ 401- 967/ 510) -56.3 E27-E19 S27-S20
     11 1565  440 (S47W23 L=273)   521    106        0 (1531/ 397-1608/ 474)  26.6 W19-W27 S49-S44
     12  541  352 (S22E55 L=195)  2318    212        7 ( 462/ 299- 638/ 412) -51.3 E63-E48 S26-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