Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1004 1783 (N41E22 L= 98)   305     58        1 ( 985/1757-1028/1800)  90.0 E23-E20 N38-N43
      2  862 1415 (N12E17 L=102)   437    136        1 ( 832/1362- 910/1483) -39.8 E20-E13 N09-N16
      3  494 1338 (N02E39 L= 80)   243     54        0 ( 480/1314- 513/1366) -45.0 E40-E37 N00-N03
      4  405 1289 (S02E44 L= 75)   934    165        2 ( 353/1214- 469/1363)  60.9 E47-E41 S07-N02
      5 1500 1189 (N13W25 L=144)  1123    311        2 (1450/1034-1529/1342) -58.0 W19-W29 N04-N20
      6  156 1134 (S14E65 L= 54)  3326    540        6 (  84/ 977- 217/1247)  59.0 E72-E60 S23-S06
      7 1730 1040 (N10W44 L=163)   473    103        0 (1711/ 983-1745/1094) -81.9 W42-W45 N06-N13
      8 1877  961 (N10W60 L=180)  1726    269        5 (1832/ 872-1905/1049) -84.8 W57-W64 N04-N15
      9 1267  481 (S33W30 L=149)   583    181        0 (1186/ 467-1361/ 506)  13.0 W23-W36 S35-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