Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  458 1605 (N36E55 L=130)   670    110        2 ( 429/1557- 494/1652) -26.6 E63-E48 N32-N39
      2 1371 1461 (N32W23 L=208)  1673    478        1 (1232/1319-1549/1559) -33.7 W12-W35 N23-N38
      3  683 1395 (N23E26 L=159)   678    219        1 ( 640/1310- 753/1475) -42.3 E30-E19 N18-N28
      4 1672 1264 (N20W46 L=230)   577    156        2 (1615/1240-1715/1318)  -6.3 W40-W49 N18-N23
      5 1729  824 (S07W51 L=236)  1244    182        1 (1657/ 758-1798/ 874)  36.0 W45-W58 S11-S03
      6 1325  726 (S15W21 L=206)   213     80        1 (1306/ 694-1342/ 766)  90.0 W19-W22 S16-S11
      7 1480  628 (S20W34 L=219)   518    110        2 (1444/ 589-1516/ 682)  66.8 W31-W36 S23-S16
      8  653  645 (S24E23 L=161)   436    116        2 ( 595/ 615- 707/ 676)  20.6 E27-E19 S26-S21
      9  470  547 (S33E40 L=144)   938    182        4 ( 380/ 532- 553/ 565)   3.2 E50-E32 S34-S30
  

    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