Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1822 (N43E35 L=184)   460     52        0 ( 964/1795-1002/1846) -33.7 E38-E32 N40-N45
      2 1366 1686 (N46W01 L=220)  1391    171        3 (1322/1608-1421/1745) -81.9 E01-W04 N39-N53
      3 1810 1361 (N40W50 L=269)   849    105        2 (1766/1322-1849/1393)  -4.4 W43-W56 N38-N42
      4 1703 1336 (N34W36 L=255)   418     78        0 (1687/1297-1717/1376) -78.7 W34-W38 N31-N36
      5 1491 1329 (N25W19 L=238)   276     73        1 (1454/1307-1524/1358)   0.0 W15-W21 N24-N26
      6  384 1183 (S13E44 L=175)   299     32        1 ( 361/1165- 410/1199)  36.9 E46-E42 S14-S11
      7  927  671 (S28W05 L=224)   454    114        0 ( 880/ 635- 992/ 703)  54.5 W02-W08 S31-S24
      8 1450  610 (S15W38 L=258)   383    146        0 (1406/ 563-1483/ 659)  81.9 W36-W40 S18-S11
      9  694  564 (S42E07 L=212)   240     32        1 ( 675/ 541- 717/ 580)  71.6 E08-E05 S43-S40
     10 1114  393 (S39W30 L=249)   749    121        1 (1064/ 354-1161/ 442)  82.9 W26-W31 S42-S34
     11 1259  309 (S38W46 L=265)  1135    272        2 (1167/ 243-1363/ 407) -86.4 W44-W49 S45-S29
     12  849  284 (S57W20 L=239)   332     54        0 ( 835/ 251- 864/ 314) -45.0 W17-W24 S59-S54
  

    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