Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  494 1717 (N21E62 L=312)  1142     78        2 ( 465/1677- 528/1749) -71.6 E66-E57 N18-N24
      2  445 1266 (S07E42 L=332)   283     87        0 ( 435/1221- 455/1306) -68.2 E43-E40 S09-S04
      3  871 1250 (N02E14 L=359)   347    114        2 ( 826/1197- 916/1290) -18.4 E17-E10 N00-N03
      4 1768  954 (N11W49 L= 62)  1919    470        2 (1679/ 778-1827/1139) -35.5 W36-W58 N02-N18
      5  523  992 (S21E30 L=344)   231     60        0 ( 502/ 965- 551/1023)  76.0 E30-E28 S22-S18
      6  835  882 (S20E07 L=  6)   219     56        0 ( 804/ 861- 869/ 903)   0.0 E09-E04 S20-S18
      7 1857  852 (N09W62 L= 75)   538     44        1 (1838/ 828-1878/ 886)  90.0 W59-W64 N07-N11
      8  666  702 (S35E14 L=360)   620    147        1 ( 596/ 658- 727/ 749)  -5.2 E19-E08 S36-S33
      9 1128  445 (S36W27 L= 40)   248     76        0 (1093/ 428-1165/ 467)   9.5 W23-W29 S37-S35
     10 1456  400 (S28W52 L= 65)   460     83        1 (1412/ 379-1498/ 426)   6.3 W47-W56 S28-S26
     11  463  379 (S66E37 L=337)   816    103        1 ( 422/ 359- 511/ 411)   8.5 E46-E26 S69-S62
     12  614  336 (S64E10 L=  3)   363     35        1 ( 590/ 314- 635/ 354)  -9.5 E13-E07 S65-S62
  

    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