Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  207 1206 (N34E57 L= 55)   502     59        1 ( 188/1181- 231/1229) -45.0 E58-E54 N32-N35
      2 1685  958 (S14W40 L=151)   424    134        2 (1628/ 921-1717/1012)  23.2 W34-W42 S16-S11
      3  221  978 (N20E57 L= 54)   356     71        1 ( 196/ 959- 247/ 998) -36.9 E59-E54 N18-N21
      4 1781  740 (S31W45 L=157)  2482    385        3 (1666/ 618-1890/ 872)  21.4 W34-W58 S35-S25
      5  918  704 (S09E14 L= 98)   358    148        1 ( 877/ 637- 978/ 768)  29.1 E17-E08 S12-S05
      6 1320  701 (S19W10 L=121)   216     81        0 (1278/ 686-1361/ 716) -36.9 W07-W12 S20-S17
      7 1576  645 (S30W26 L=138)   525    126        2 (1544/ 586-1608/ 711)  54.5 W23-W29 S33-S26
      8  918  634 (S13E16 L= 96)   487    136        2 ( 869/ 592- 963/ 683)  14.0 E19-E11 S15-S10
      9 1095  550 (S22E08 L=104)   344    109        0 (1058/ 516-1147/ 595) -66.8 E09-E05 S25-S18
     10  658  441 (S19E39 L= 73)   311     86        2 ( 626/ 422- 697/ 461) -11.3 E41-E36 S21-S18
     11  633  303 (S29E49 L= 63)   380    111        0 ( 615/ 273- 659/ 346)  63.4 E50-E46 S31-S25
  

    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