Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1077 1736 (N39E21 L= 96)  2440    442        3 ( 935/1606-1172/1831) -19.2 E35-E11 N30-N44
      2  460 1605 (N13E55 L= 62)  3948    738       11 ( 275/1446- 599/1740)  65.4 E67-E46 N00-N24
      3  692 1478 (N11E32 L= 85)  1961    562        5 ( 577/1355- 834/1563)  12.5 E42-E24 N04-N16
      4  401 1472 (N04E54 L= 64)   419    112        1 ( 354/1447- 446/1488)  29.7 E57-E49 N02-N06
      5 1501 1250 (N21W22 L=139)   885    184        4 (1458/1196-1546/1315) -26.6 W17-W25 N17-N23
      6  265 1186 (S15E55 L= 63)   490    122        0 ( 231/1130- 289/1241)  74.1 E56-E53 S18-S11
      7 1604  443 (S19W59 L=176)   488     92        1 (1566/ 424-1649/ 462)  45.0 W56-W61 S21-S16
      8 1206  249 (S44W46 L=163)  1616    282        2 (1081/ 221-1330/ 284)  24.0 W35-W53 S48-S39
  

    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