Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  478 1544 (N23E43 L= 87)   401     92        1 ( 455/1507- 508/1579) -33.7 E46-E40 N20-N25
      2  399 1439 (N16E47 L= 83)   582    133        0 ( 363/1393- 423/1474) -38.7 E50-E44 N12-N18
      3 1213 1185 (N06W10 L=140)   445    108        0 (1192/1135-1246/1248) -71.6 W08-W12 N03-N09
      4 1378 1199 (N09W20 L=150)   200     65        0 (1364/1167-1391/1232) -76.0 W19-W20 N06-N10
      5  682 1042 (S07E21 L=109)   217     55        1 ( 665/1014- 704/1081)  63.4 E21-E19 S08-S04
      6  545 1028 (S08E30 L=100)   696    178        1 ( 482/ 993- 636/1080) -10.3 E34-E23 S10-S05
      7   92 1017 (S10E75 L= 55)   873    100        2 (  75/ 980- 104/1057)  90.0 E78-E72 S11-S07
      8  550  747 (S25E29 L=101)   425    107        1 ( 532/ 692- 571/ 814) -90.0 E30-E27 S28-S21
      9 1838  675 (S16W65 L=195)   695    117        0 (1795/ 616-1870/ 726)  54.5 W62-W68 S19-S12
     10 1380  557 (S30W30 L=160)   720    141        2 (1331/ 511-1446/ 612)  58.0 W26-W33 S33-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