Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  801 1858 (N48E48 L=340)  1039    186        2 ( 766/1819- 856/1886)  -9.5 E53-E41 N43-N50
      2  579 1459 (N14E38 L=349)  1675    430        7 ( 486/1345- 691/1580) -23.5 E50-E27 N09-N19
      3 1546 1325 (N28W27 L= 55)   658    211        0 (1478/1240-1604/1428) -39.3 W21-W32 N23-N32
      4 1197 1296 (N17W04 L= 32)   432     92        2 (1161/1263-1242/1336) -31.0 W02-W07 N15-N19
      5  565 1149 (S05E30 L=357)   434    114        1 ( 509/1110- 610/1169)  39.8 E33-E27 S07-S02
      6 1219 1138 (N08W09 L= 36)   283     23        0 (1196/1117-1241/1164)   0.0 W07-W10 N07-N09
      7  828  969 (S10E10 L= 17)   530    201        1 ( 732/ 937- 936/1006)   8.7 E16-E03 S11-S07
      8  204  806 (S33E57 L=331)   965    263        0 ( 160/ 753- 255/ 861)  39.3 E62-E51 S36-S27
      9  907  528 (S35W05 L= 32)   219     52        0 ( 879/ 518- 943/ 542)  26.6 W02-W07 S36-S33
  

    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