Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  492 1533 (N51E33 L=246)   290     56        1 ( 475/1503- 512/1561)  53.1 E34-E30 N48-N52
      2 1236 1516 (N31W27 L=306)  2313    765        2 ( 947/1417-1564/1552) -30.1 W07-W45 N16-N39
      3 1844 1217 (S03W63 L=342)   572     86        2 (1823/1176-1874/1256) -63.4 W59-W65 S06-N00
      4 1189 1210 (N14W14 L=293)   424    271        0 (1081/1140-1276/1249) -42.0 W07-W17 N08-N17
      5  336 1088 (N23E46 L=233)  1004    245        2 ( 298/ 998- 384/1200) -72.9 E49-E42 N17-N30
      6  834 1092 (N15E10 L=269)   777    241        3 ( 808/ 958- 866/1186)  85.9 E11-E07 N06-N20
      7 1775  785 (S26W47 L=326)  1283    234        2 (1706/ 687-1825/ 881)  31.6 W40-W53 S31-S21
      8  427  703 (S01E45 L=234)   841    199        2 ( 399/ 612- 489/ 765)  37.9 E50-E39 S05-N02
      9  364  527 (S11E59 L=220)  6849    543       13 ( 264/ 384- 468/ 663) -81.5 E71-E50 S21-S01
     10 1311  404 (S39W04 L=283)   300     48        0 (1285/ 379-1334/ 431)  26.6 W02-W06 S41-S37
  

    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