Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  742 1726 (N33E41 L= 94)   418     58        2 ( 719/1703- 764/1753) -68.2 E42-E39 N30-N35
      2 1191 1680 (N42E08 L=126)  1197    206        2 (1145/1601-1248/1752) -34.7 E14-E01 N37-N47
      3  787 1493 (N19E26 L=108)   707    213        2 ( 740/1401- 844/1589) -81.3 E28-E24 N12-N25
      4  508 1504 (N13E46 L= 88)   575    112        0 ( 477/1439- 544/1539)  81.9 E47-E44 N08-N15
      5  568 1508 (N15E42 L= 93)   272     54        1 ( 550/1479- 594/1526)  76.0 E43-E40 N12-N16
      6  930 1387 (N16E14 L=120)   242     83        0 ( 895/1349- 954/1427)  80.5 E15-E13 N12-N18
      7  174 1252 (S08E69 L= 65)  2277    307        7 ( 121/1198- 207/1307)  33.7 E77-E64 S11-S04
      8 1191 1162 (N09W06 L=141)   217     56        1 (1167/1139-1216/1194) -18.4 W04-W07 N07-N09
      9  623  649 (S35E17 L=117)   276     34        1 ( 606/ 621- 644/ 680) -45.0 E18-E15 S36-S33
     10 1276  319 (S39W41 L=176)  1590    376        2 (1140/ 292-1468/ 350)  24.4 W30-W53 S42-S31
  

    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