Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1044 1768 (N40E22 L=262)   306     28        0 (1020/1754-1070/1783)  33.7 E23-E19 N38-N41
      2 1136 1709 (N38E13 L=271)   743    128        1 (1089/1657-1186/1748) -18.4 E16-E07 N35-N40
      3  788 1666 (N25E33 L=251)   445    116        0 ( 766/1615- 806/1716) -60.3 E35-E30 N21-N28
      4  516 1591 (N15E48 L=236)  1548    263        2 ( 494/1489- 547/1676) -47.5 E53-E41 N08-N20
      5 1095 1317 (N11E03 L=281)  1409    400        2 (1057/1143-1127/1528) -66.5 E08-W02 N01-N24
      6  624 1412 (N06E33 L=251)   780    190        1 ( 548/1374- 697/1466)  -9.5 E39-E27 N04-N07
      7 1823 1041 (N16W52 L=335)  1185    157        1 (1787/ 968-1858/1114) -41.2 W47-W56 N12-N19
      8  322  979 (S24E45 L=239)  1168    283        1 ( 228/ 908- 411/1052)  35.5 E53-E38 S30-S19
      9 1820  831 (N04W58 L=342)   880    179        3 (1767/ 791-1858/ 879)  45.0 W53-W62 N01-N08
     10  136  839 (S35E63 L=221)   723    254        0 ( 105/ 803- 165/ 875)  35.5 E66-E59 S37-S32
     11 1448  592 (S20W38 L=322)   336     79        0 (1421/ 550-1468/ 633) -78.7 W36-W39 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