Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  947 1431 (N28W05 L=335)   295    112        0 ( 903/1396- 998/1464) -56.3 W02-W07 N24-N30
      2 1824 1259 (S03W59 L= 29)  4232    363        3 (1727/1068-1896/1361) -80.0 W52-W69 S13-N04
      3 1682 1302 (N03W48 L= 17)   435    103        0 (1636/1268-1719/1330)   6.3 W42-W51 N02-N04
      4  902 1313 (N22E01 L=329)   232     68        0 ( 861/1302- 942/1327) -26.6 E03-W01 N20-N22
      5 1313 1288 (N11W23 L=352)   351     96        0 (1264/1266-1363/1306) -26.6 W19-W25 N08-N13
      6 1930 1143 (S12W70 L= 40)   680     79        1 (1919/1109-1948/1187)  63.4 W67-W72 S14-S09
      7 1538  955 (S12W29 L=358)   223     20        0 (1518/ 935-1558/ 977) -63.4 W27-W30 S13-S11
      8 1284  664 (S24W07 L=337)  3241   1002        3 ( 914/ 387-1614/ 887)  15.3 E23-W32 S32-S17
      9 1054  791 (S11E03 L=326)   443    170        0 ( 973/ 757-1127/ 838)   0.0 E08-W01 S11-S09
     10  205  779 (N05E63 L=267)  2230    233        4 ( 121/ 730- 322/ 833)  -8.4 E79-E51 N00-N08
     11  532  304 (S30E59 L=271)  1599    160        4 ( 464/ 251- 603/ 344)   7.1 E73-E49 S32-S27
  

    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