Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1075 1340 (N22W12 L=109)   450    174        1 (1055/1274-1116/1427)  77.5 W10-W14 N17-N27
      2  758 1207 (N23E10 L= 86)  2085    588        5 ( 506/1039- 941/1332) -51.8 E25-E02 N09-N37
      3 1695  992 (S15W39 L=136)  2389    704        0 (1514/ 750-1797/1246)  28.6 W21-W54 S23-S04
      4 1947 1146 (S16W74 L=171)   992     71        1 (1932/1107-1964/1183) -63.4 W71-W79 S19-S13
      5 1875  958 (S24W57 L=154)  1595    173        4 (1814/ 884-1923/1019)  10.6 W48-W65 S26-S21
      6 1354  749 (S18W11 L=108)   992    292        0 (1299/ 619-1411/ 870)  42.7 W04-W17 S24-S12
      7 1064  726 (S11E06 L= 91)   301    156        1 (1008/ 668-1090/ 765)  51.3 E08-E03 S14-S09
      8  615  512 (S12E39 L= 58)   755    195        1 ( 593/ 427- 634/ 601)  54.0 E44-E35 S17-S06
      9  301  534 (S05E65 L= 32)  1466    132        3 ( 264/ 482- 354/ 578)  40.6 E71-E59 S07-S01
     10  523  520 (S10E45 L= 52)   376     46        1 ( 495/ 499- 552/ 541) -26.6 E47-E43 S10-S08
     11  471  393 (S17E58 L= 39)  1433    177        2 ( 434/ 322- 523/ 480)  67.4 E61-E54 S22-S10
  

    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