Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  484 1334 (N35E33 L= 90)   322     86        0 ( 445/1306- 519/1359) -29.7 E36-E29 N33-N37
      2 1727 1018 (S09W46 L=170)  1050    289        2 (1644/ 961-1815/1062)   3.4 W38-W55 S12-S06
      3 1776  900 (S17W49 L=173)   783    222        0 (1730/ 866-1822/ 956)  18.4 W44-W54 S19-S13
      4 1897  915 (S21W65 L=189)   552     77        0 (1869/ 884-1921/ 954)  12.5 W60-W70 S22-S19
      5 1433  927 (S07W23 L=146)   211     57        1 (1400/ 902-1458/ 950) -45.0 W20-W24 S09-S05
      6 1524  822 (S16W27 L=151)   272     51        0 (1508/ 788-1543/ 852)  76.0 W25-W28 S17-S13
      7  913  722 (S09E12 L=111)   220     44        0 ( 898/ 694- 927/ 748)  56.3 E13-E11 S10-S07
      8 1084  654 (S16E04 L=120)   628    269        1 (1052/ 551-1172/ 717)  37.9 E07-W02 S22-S12
      9  437  544 (S13E51 L= 73)   406     77        0 ( 417/ 505- 459/ 582)  35.5 E54-E47 S15-S10
     10 1521  448 (S42W25 L=149)   391     40        1 (1500/ 430-1542/ 472)  14.0 W23-W27 S44-S39
     11 1416  376 (S45W15 L=139)   369     61        1 (1400/ 351-1442/ 408) -80.5 W13-W17 S47-S41
  

    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