Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  536 1559 (N15E47 L= 18)   376     98        0 ( 508/1522- 576/1592) -12.5 E51-E42 N13-N16
      2  328 1486 (N06E62 L=  3)   466    104        0 ( 299/1442- 346/1532) -74.1 E64-E60 N02-N09
      3  379 1493 (N07E57 L=  8)   584     86        0 ( 356/1456- 413/1520)  53.1 E60-E54 N04-N09
      4 1500 1411 (N32W20 L= 85)   480    113        2 (1460/1363-1531/1458) -39.8 W16-W23 N29-N34
      5  170 1101 (S20E66 L=359)  1172    190        3 ( 129/1025- 227/1164)  45.0 E73-E59 S24-S14
      6  373  732 (S38E40 L= 26)   397    119        0 ( 323/ 710- 420/ 758)  15.3 E45-E34 S39-S35
      7  696  617 (S37E10 L= 56)  1925    350        5 ( 575/ 530- 834/ 734)  53.1 E18-E03 S45-S25
      8 1703  585 (S09W60 L=125)  1540    330        2 (1606/ 511-1804/ 657)  47.3 W54-W67 S15-S02
      9 1503  552 (S17W45 L=110)   280     74        0 (1478/ 521-1539/ 586)  90.0 W43-W45 S19-S14
     10 1352  389 (S31W44 L=110)  1499    321        4 (1231/ 330-1486/ 462)  45.0 W37-W50 S37-S23
  

    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