Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  900 1826 (N44E41 L=322)  1260    154        2 ( 836/1777- 954/1860)  -7.1 E49-E32 N41-N46
      2  506 1564 (N15E50 L=313)   864    229        2 ( 439/1482- 553/1627) -22.4 E59-E42 N11-N18
      3  311 1514 (N09E67 L=296)   544     91        0 ( 294/1472- 334/1554) -63.4 E68-E64 N05-N11
      4  569  908 (S21E25 L=338)   516    224        0 ( 518/ 833- 630/ 996) -38.7 E29-E19 S24-S16
      5  801  876 (S17E09 L=354)   214     75        2 ( 777/ 838- 823/ 911) -45.0 E11-E07 S19-S15
      6  969  644 (S27W07 L= 10)   506    168        2 ( 884/ 625-1036/ 662)  21.8 W01-W11 S29-S24
      7 1537  363 (S27W62 L= 65)   707     75        1 (1504/ 337-1580/ 384)   5.2 W56-W67 S28-S25
      8 1389  274 (S37W59 L= 61)  2364    198        5 (1315/ 228-1483/ 325)  56.3 W48-W67 S42-S32
      9 1264  269 (S42W49 L= 52)   349     37        0 (1251/ 246-1281/ 301) -26.6 W44-W52 S43-S39
  

    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