Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  947 1593 (N41E07 L=342)   488     88        2 ( 918/1554- 982/1629)  90.0 E09-E04 N37-N43
      2  324 1416 (N25E55 L=294)   422    138        1 ( 297/1379- 352/1457) -33.7 E58-E51 N23-N27
      3 1687 1382 (N26W49 L= 38)  3100    722        9 (1455/1262-1890/1455) -16.7 W30-W70 N17-N30
      4  975 1311 (N21E04 L=345)   526    143        0 ( 945/1241-1005/1380) -66.0 E05-E01 N16-N25
      5 1688 1203 (N14W46 L= 35)  1071    236        0 (1591/1115-1751/1291) -42.5 W39-W51 N08-N20
      6 1310 1072 (N06W18 L=  7)   479    155        1 (1281/1004-1349/1155) -71.6 W16-W20 N02-N11
      7 1439  982 (N01W26 L= 16)   869    194        2 (1354/ 946-1506/1034)  21.8 W20-W31 S01-N04
      8 1238  908 (S04W13 L=  3)   225     74        2 (1201/ 887-1272/ 925)   0.0 W11-W15 S04-S02
  

    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