Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  859 1549 (N32E19 L=287)   501     80        1 ( 826/1504- 899/1575)  63.4 E21-E15 N29-N34
      2 1806 1360 (N30W62 L=  8)   782    122        1 (1780/1297-1835/1403) -66.8 W59-W65 N26-N33
      3 1065 1216 (N13W00 L=306)   399    131        1 (1047/1144-1084/1277) -82.9 E01-W01 N08-N16
      4 1475 1127 (N12W28 L=334)   546    206        1 (1439/1035-1504/1214) -65.6 W24-W30 N06-N17
      5 1725 1077 (N11W48 L=354)  1551    297        1 (1689/ 940-1751/1202)  76.8 W45-W50 N02-N19
      6  578  942 (S09E27 L=278)   599    125        0 ( 551/ 884- 601/1003)  81.9 E29-E25 S12-S05
      7  176  917 (S15E65 L=241)   518     68        2 ( 151/ 885- 202/ 944) -31.0 E69-E61 S17-S13
      8  183  816 (S22E66 L=239)  3089    130        4 ( 142/ 703- 227/ 911)  45.0 E75-E58 S29-S15
      9 1639  725 (S11W46 L=352)   456    136        0 (1610/ 667-1659/ 790)  90.0 W44-W47 S14-S06
     10 1376  606 (S21W29 L=335)   241    124        2 (1336/ 569-1411/ 661)  71.6 W26-W30 S23-S17
     11 1688  478 (S27W64 L= 10)  4113    337        4 (1549/ 401-1803/ 577)  36.4 W50-W75 S33-S18
  

    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