Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1658 (N23E48 L=176)   440     91        0 ( 508/1636- 594/1680)  38.7 E51-E45 N21-N25
      2  898 1528 (N20E19 L=206)   694    207        0 ( 879/1434- 919/1633) -65.2 E22-E14 N14-N27
      3 1507 1515 (N34W22 L=246)  1434    338        1 (1388/1386-1589/1626) -35.2 W11-W29 N27-N39
      4  340  910 (S25E45 L=180)  1316    306        0 ( 287/ 829- 459/1027) -35.5 E50-E33 S29-S18
      5 1301  913 (S07W18 L=243)   248    110        0 (1276/ 860-1314/ 961) -51.3 W15-W19 S10-S05
      6  147  820 (S33E71 L=153)   921    112        1 ( 120/ 791- 192/ 871) -32.0 E77-E63 S34-S29
      7  454  476 (S53E39 L=186)   466     51        1 ( 425/ 460- 486/ 502)   0.0 E42-E34 S54-S51
      8 1292  302 (S43W45 L=269)   526     58        1 (1256/ 276-1321/ 326)  36.9 W41-W48 S44-S41
      9 1544  283 (S36W72 L=296)   938     88        1 (1505/ 266-1587/ 307)  10.6 W65-W82 S37-S33
  

    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