Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  917 1757 (N40E31 L=126)   811     97        4 ( 873/1733- 967/1782)  33.7 E35-E27 N37-N42
      2  546 1714 (N27E57 L=100)   512     79        0 ( 524/1673- 574/1749) -18.4 E62-E50 N25-N29
      3 1007 1661 (N35E19 L=138)  1465    227        2 ( 960/1566-1066/1743) -31.0 E26-E11 N29-N40
      4 1610 1598 (N50W33 L=190)  1389    164        2 (1544/1555-1689/1636)   0.0 W24-W43 N47-N52
      5  431 1303 (S00E44 L=114)  1630    281        2 ( 388/1170- 484/1405) -42.5 E48-E36 S06-N05
      6 1298 1324 (N21W09 L=166)   473    165        0 (1233/1262-1346/1375) -24.0 W04-W13 N17-N22
      7  930 1088 (S02E07 L=150)   192     35        0 ( 913/1065- 950/1114) -33.7 E08-E05 S02-N00
      8  547  961 (S17E28 L=130)   248     41        0 ( 530/ 934- 568/ 984)  90.0 E28-E26 S18-S15
      9 1044  826 (S14W06 L=163)  1023    317        0 ( 909/ 755-1182/ 874)  45.0 W00-W12 S20-S08
     10 1499  794 (S05W34 L=191)   334    108        1 (1464/ 747-1520/ 836) -26.6 W30-W36 S07-S02
     11 1175  599 (S24W20 L=177)  1178    262        2 (1058/ 578-1307/ 639)  20.6 W12-W28 S26-S20
     12  460  566 (S45E32 L=126)   486     43        1 ( 429/ 543- 494/ 592)  20.6 E35-E27 S47-S42
     13 1537  365 (S28W60 L=217)   408     68        0 (1522/ 341-1563/ 397) -23.2 W56-W63 S29-S25
  

    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