Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1044 1582 (N34W09 L=193)   235     53        0 (1016/1564-1074/1599)  14.0 W07-W11 N32-N34
      2 1416 1364 (N15W29 L=213)   240    120        0 (1370/1326-1465/1384)  15.9 W25-W32 N13-N15
      3 1806 1257 (N04W58 L=241)   649    142        1 (1772/1185-1848/1322) -69.4 W55-W61 N00-N08
      4  191 1010 (N09E59 L=125)   895    168        3 ( 158/ 943- 234/1083) -60.9 E62-E55 N04-N13
      5  302  905 (N01E50 L=134)   882    173        0 ( 269/ 833- 323/ 997)  53.1 E53-E46 S02-N06
      6 1346  923 (S10W18 L=202)   216     81        0 (1299/ 900-1376/ 941) -26.6 W15-W19 S11-S08
      7 1161  797 (S15W05 L=189)   212     66        2 (1118/ 777-1196/ 817)   0.0 W02-W07 S16-S14
      8  311  640 (S14E55 L=128)  1630    291        6 ( 252/ 551- 372/ 735) -90.0 E60-E50 S19-S08
  

    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