Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1041 1771 (N40E22 L=263)   317     38        0 (1014/1756-1071/1784)  45.0 E23-E20 N38-N41
      2 1127 1713 (N38E13 L=271)   672     95        2 (1084/1678-1175/1749)  -7.1 E16-E08 N35-N40
      3  778 1671 (N26E33 L=251)   344    101        0 ( 762/1627- 796/1710) -63.4 E35-E31 N22-N28
      4  548 1520 (N11E42 L=242)  2091    487        2 ( 484/1376- 667/1670) -32.0 E54-E29 N04-N20
      5 1083 1327 (N11E04 L=281)  1351    393        2 (1044/1150-1116/1529) -67.8 E08-W01 N01-N23
      6  112 1115 (S19E70 L=214)  3723    432        5 (  75/1038- 152/1197)  77.5 E77-E65 S23-S13
      7 1815 1046 (N16W51 L=335)  1013    149        1 (1776/ 983-1851/1114) -26.6 W45-W55 N12-N18
      8  315  983 (S24E46 L=239)  1146    273        1 ( 230/ 915- 400/1057)  37.6 E53-E39 S29-S19
      9 1817  830 (N04W58 L=342)   820    183        2 (1776/ 798-1866/ 868)  16.7 W54-W64 N01-N07
     10 1436  595 (S20W37 L=322)   306     72        0 (1413/ 554-1454/ 636) -78.7 W35-W38 S22-S17
  

    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