Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1415 1584 (N25W42 L=212)  3824    715        9 (1116/1525-1684/1618) -22.5 W22-W63 N15-N33
      2  942 1438 (N30W07 L=177)   203     40        0 ( 912/1420- 969/1459) -12.8 W05-W08 N28-N32
      3 1896 1260 (S07W70 L=240)   733    188        3 (1870/1238-1930/1286) -15.3 W65-W76 S09-S05
      4 1918 1149 (S14W67 L=237)  3476    260        5 (1875/1062-1956/1225)  78.7 W62-W74 S19-S09
      5 1724  781 (S29W39 L=209)  2698    515        7 (1576/ 602-1785/1005)  39.3 W24-W47 S37-S17
      6  676  910 (N07E23 L=148)   256     66        0 ( 647/ 883- 713/ 928)  14.0 E24-E20 N05-N07
      7 1586  774 (S25W28 L=198)   723    217        1 (1520/ 666-1648/ 878)  34.7 W21-W34 S29-S20
      8 1289  610 (S26W05 L=175)   361    146        1 (1241/ 538-1317/ 675)  45.0 W00-W07 S29-S22
      9  933  529 (S21E19 L=151)   630    254        0 ( 844/ 447-1036/ 623)  17.4 E26-E10 S24-S18
     10  371  415 (S16E65 L=105)  1058    212        3 ( 342/ 379- 390/ 452)  24.0 E70-E61 S18-S13
     11  838  356 (S30E32 L=138)   382     87        1 ( 802/ 326- 878/ 383)   8.1 E35-E28 S31-S28
  

    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