Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1208 1569 (N35W27 L=171)   409     90        1 (1166/1541-1259/1589)  -8.1 W23-W30 N33-N36
      2 1010 1528 (N38W12 L=155)   221     69        0 ( 991/1502-1027/1560)  31.0 W09-W14 N36-N39
      3  566 1379 (N38E25 L=119)   832    147        2 ( 518/1320- 622/1431) -45.0 E28-E20 N33-N41
      4  330 1184 (N29E46 L= 98)  1066    180        3 ( 285/1109- 377/1254) -59.0 E49-E41 N23-N33
      5 1768 1175 (S02W53 L=197)   379     72        1 (1738/1147-1790/1193) -21.8 W50-W55 S03-N00
      6  847 1147 (N18E08 L=135)   497    166        1 ( 773/1113- 936/1192) -31.0 E12-E02 N15-N21
      7 1709 1153 (S02W47 L=191)   349     64        1 (1690/1117-1728/1187)  45.0 W44-W49 S03-N00
      8 1093  925 (N00W02 L=145)  1081    426        2 (1014/ 804-1144/1026) -90.0 W00-W04 S07-N07
      9 1706  934 (S14W42 L=185)   761    112        2 (1679/ 868-1734/ 986)  56.3 W39-W44 S18-S11
     10 1833  912 (S19W54 L=197)   441    180        0 (1766/ 868-1856/ 960)  16.7 W46-W56 S21-S17
     11 1058  636 (S17E06 L=137)   753    144        2 ( 986/ 577-1097/ 670) -60.3 E09-E03 S20-S13
     12  848  314 (S35E31 L=113)  1047    174        3 ( 806/ 261- 886/ 366)  35.0 E36-E25 S38-S31
  

    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