Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1153 1554 (N31W13 L=132)  1256    318        0 (1034/1464-1270/1645) -45.0 W05-W19 N23-N38
      2 1406 1584 (N31W32 L=151)   603    133        0 (1342/1547-1467/1629) -45.0 W27-W35 N28-N34
      3  404 1349 (N22E41 L= 78)   906    134        0 ( 373/1300- 456/1407) -41.2 E44-E36 N18-N26
      4  161 1343 (N24E72 L= 47)  1436     97        1 ( 133/1301- 186/1391)  90.0 E78-E68 N21-N27
      5 1837 1132 (N01W59 L=178)   972    168        1 (1762/1081-1893/1167) -13.0 W51-W66 S02-N02
      6 1870  873 (S15W64 L=183)  1253    164        1 (1806/ 809-1914/ 934)  30.3 W56-W69 S18-S11
      7  782  850 (S11E16 L=103)   383    124        0 ( 739/ 806- 833/ 896)  35.5 E19-E12 S13-S08
      8 1677  606 (S31W48 L=167)  3611    571        1 (1495/ 444-1883/ 758)  28.6 W34-W70 S42-S21
      9  823  623 (S25E16 L=103)   626    197        1 ( 798/ 558- 884/ 711) -78.7 E17-E12 S29-S19
     10  505  511 (S30E43 L= 76)   409    137        0 ( 452/ 496- 583/ 527)  -5.2 E47-E36 S31-S28
     11 1373  471 (S39W23 L=143)   666    113        1 (1339/ 415-1401/ 521) -49.4 W20-W27 S43-S35
     12 1387  253 (S59W35 L=154)   500     48        2 (1359/ 238-1418/ 273)   0.0 W31-W41 S60-S56
  

    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